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

([7327426.165155034,
  7799377.79122286,
  7799390.667111297,
  7939804.684034693,
  7956131.510062348,
  8109033.885006448,
  8209116.604035509,
  8216208.976214902,
  8335283.104312408,
  8486787.5,
  8556695.013555547,
  8608280.537665417,
  8789835.9375,
  8830321.875,
  8893587.5,
  8893589.957035514,
  8909754.64975903,
  8921381.17891598,
  8922395.911433162,
  8968225.095157664,
  8968268.254809383,
  8969201.160707876,
  8970016.615093034,
  8979272.431575812,
  8979939.62047793,
  8993147.271797227,
  9007339.0625,
  9009068.049013903,
  9012516.867536742,
  9033400.48309831,
  9055802.4091866,
  9095473.4375,
  9095495.205135133,
  9106713.365973506,
  9131542.94382511,
  9143060.9375,
  9154457.631342614,
  9154467.82935047,
  9158136.056151805,
  9166148.13453587,
  9200683.784165535,
  9217205.438665004,
  9217976.397194602,
  9220791.782393917,
  9230026.246870363,
  9246042.1875,
  9251360.641888585,
  9278027.71716614,
  9285183.01007829,
  9296638.925568322,
  9299884.049148375,
  9345634.234690037,
  9346292.1875,
  9346492.1875,
  9347976.003853073,
  9348831.25,
  9348848.4375,
  9349648.16009287,
  9352925.251617005,
  9353805.306362716,
  9400261.237052396,
  9401669.840147445,
  9422966.644377237,
  9423806.690948995,
  9432250.182095895,
  9432695.216365414,
  9433114.0625,
  9443974.70048084,
  9446769.984196296,
  9450565.625,
  9453778.125,
  9455848.4375,
  9460913.985221656,
  9463617.348407395,
  9466659.647269633,
  9483931.777406503,
  9508966.540260786,
  9514589.957339026,
  9514901.62147496,
  9517695.3125,
  9521261.893517284,
  9522652.604017746,
  9523066.422492987,
  9523903.125,
  9524551.5625,
  9524792.721884431,
  9525302.290748097,
  9528440.391628487,
  9552768.36184259,
  9559687.178782301,
  9559999.8436646,
  9560656.25,
  9561341.903074494,
  9579026.212965334,
  9589400.645778496,
  9595899.504652567,
  9599269.238749076,
  9605294.50687347,
  9609730.7235725,
  9610199.167961268,
  9613146.628884837,
  9614226.275596865,
  9618464.663514046,
  9664529.894618401,
  9664663.465508314,
  9665144.442758262,
  9668211.700415751,
  9700166.790179495,
  9723083.913256293,
  9730400.555775732,
  9736848.83526096,
  9741232.42610989,
  9742377.981699638,
  9742823.4375,
  9765887.035624983,
  9766693.75,
  9766719.67725108,
  9770410.868924772,
  9780678.125,
  9781144.204739803,
  9785720.431165475,
  9788358.655751918,
  9808085.9375,
  9808583.588169497,
  9813212.39343075,
  9815598.193608623,
  9824971.875,
  9826076.701176478,
  9831420.93589777,
  9835396.19625881,
  9835639.0625,
  9835970.3125,
  9836261.35081455,
  9836515.450608224,
  9837967.475851655,
  9838809.375,
  9856243.75,
  9856481.048540115,
  9885837.5,
  9887371.517387621,
  9895981.25,
  9898386.723926332,
  9900897.107390271,
  9901403.346210882,
  9901406.651192874,
  9907061.564476078,
  9912917.60361787,
  9916457.8125,
  9916729.719260572,
  9944321.875,
  9945217.727706514,
  9945235.571232207,
  9950724.472448848,
  9953787.713318588,
  9964235.9375,
  9964544.201671857,
  9965747.700305002,
  9966795.3125,
  9968779.6875,
  9969680.755824512,
  9970300.0,
  9974410.9375,
  9974676.19605752,
  9986223.151097342,
  9990026.5625,
  9990045.3125,
  9999723.843523532,
  10003210.970740348,
  10007898.4375,
  10007951.926766513,
  10012539.194053937,
  10017620.1805305,
  10022336.710041035,
  10025506.25,
  10026787.009956984,
  10027405.22051793,
  10030017.1875,
  10034018.71479335,
  10043229.35362895,
  10064962.6623484,
  10064971.687430203,
  10071823.652972111,
  10078284.375,
  10080052.385105353,
  10081281.503864676,
  10086465.625,
  10089349.386021024,
  10090369.992063241,
  10095773.4375,
  10096840.183617594,
  10097425.0,
  10107094.355760887,
  10126921.832971035,
  10127321.875,
  10142575.362841457,
  10149932.273067508,
  10160179.6875,
  10180140.60061183,
  10180251.216514101,
  10183453.125,
  10190345.236790597,
  10201112.5,
  10223676.932226567,
  10224796.885978809,
  10236476.5625,
  10248871.140190983,
  10252125.59280453,
  10255750.0,
  10257955.108714499,
  10270148.309936251,
  10279271.600184523,
  10292621.210173285,
  10294989.041992733,
  10295976.5625,
  10296107.141611105,
  10296147.18283717,
  10308536.8888846,
  10312792.1875,
  10322008.61804829,
  10325643.714425078,
  10326206.25,
  10358231.25,
  10360787.5,
  10364907.8125,
  10370717.1875,
  10380779.60903569,
  10387392.1875,
  10387644.272717165,
  10388113.212121796,
  10389802.23951119,
  10410451.5625,
  10416700.0,
  10422747.914414411,
  10424831.648875559,
  10425022.380816001,
  10426523.468860503,
  10427165.625,
  10431043.770637976,
  10432017.376520744,
  10440187.61989328,
  10445104.6875,
  10445145.3125,
  10448706.25,
  10454005.239069661,
  10455601.5625,
  10455634.375,
  10458771.670926899,
  10458947.556441935,
  10458978.939359512,
  10461587.5,
  10461907.43368539,
  10462978.125,
  10462995.742463369,
  10465705.108299486,
  10465730.334377792,
  10474415.157231178,
  10477700.0,
  10478309.519962726,
  10479475.971809018,
  10481060.883694548,
  10487465.136796782,
  10495787.160922995,
  10501945.3125,
  10503381.006701568,
  10508401.175898485,
  10508667.1875,
  10509104.6875,
  10513143.679788072,
  10518677.883292628,
  10529090.275011314,
  10529220.572978143,
  10529242.73658636,
  10531850.246684035,
  10531909.937488614,
  10534096.875,
  10538176.38933034,
  10539913.423931781,
  10542124.748207744,
  10544145.3125,
  10550090.476358475,
  10565461.025291791,
  10566777.551967675,
  10567358.185831482,
  10568051.5625,
  10580235.826049905,
  10585429.060469655,
  10591216.259630535,
  10595089.806242248,
  10596014.607546404,
  10598147.886318216,
  10603913.213103363,
  10607828.399749015,
  10620724.530082503,
  10621599.355544081,
  10621700.723934634,
  10625354.53769265,
  10626814.546452133,
  10642616.993597535,
  10642817.068049252,
  10650596.99334149,
  10655201.5625,
  10656549.276483875,
  10662374.567019334,
  10663439.084789084,
  10665101.091205105,
  10667427.230182495,
  10667959.943337955,
  10668053.887993434,
  10672405.717701185,
  10683301.5625,
  10686336.977740964,
  10696645.3125,
  10698866.11348568,
  10708080.354843566,
  10708385.9375,
  10711884.375,
  10731785.761690062,
  10734933.554308563,
  10736203.125,
  10742700.892635493,
  10745604.6875,
  10746926.881369725,
  10748553.125,
  10751181.103338206,
  10751209.941040728,
  10760750.66745447,
  10767169.851767238,
  10767230.363151476,
  10785190.625,
  10787125.523350062,
  10787492.1875,
  10805005.962528912,
  10805239.0625,
  10806127.018236192,
  10807394.87982792,
  10810114.637627069,
  10812429.6875,
  10815804.6875,
  10822979.66390588,
  10857645.227232529,
  10863868.75,
  10918338.336707894,
  10919619.884979524,
  10945764.0625,
  10964353.57477003,
  10994955.989037808,
  10998174.385874284,
  10998653.125,
  11002853.593632428,
  11006450.996597279,
  11016062.26404281,
  11018686.699659476,
  11025283.719425179,
  11032134.266904904,
  11053821.941399861,
  11095335.763671527,
  11125766.250932002,
  11132193.025735555,
  11150030.413336134,
  11155254.404516628,
  11164391.960484443,
  11166353.799916288,
  11176426.433260964,
  11185793.75,
  11232250.489999108,
  11235793.198009068,
  11236318.40628443,
  11237661.520932302,
  11255545.455682302,
  11255893.971820354,
  11263565.344589077,
  11264400.054062905,
  11274832.30043398,
  11278684.95155125,
  11280390.625,
  11285168.75,
  11287679.6875,
  11312735.47120625,
  11324629.247025544,
  11324632.065573448,
  11328610.003886038,
  11329234.145130951,
  11337037.200278776,
  11357041.246543927,
  11375579.13056528,
  11375700.0,
  11403106.99226797,
  11405376.890927367,
  11413892.410228908,
  11415845.069434233,
  11416109.375,
  11442310.9375,
  11473715.625,
  11478972.423351046,
  11484207.64891236,
  11516124.755877858,
  11522331.25,
  11536606.746160336,
  11560892.797188083,
  11583493.392771447,
  11589855.767051289,
  11617951.5625,
  11620610.044228334,
  11620613.713904215,
  11623857.8125,
  11623964.0625,
  11625526.664737193,
  11625709.210187174,
  11656000.560084011,
  11657281.627008902,
  11657748.305358363,
  11668675.143138427,
  11669969.509165164,
  11739986.777641645,
  11745700.0,
  11759813.869591098,
  11759853.078272732,
  11768611.953919861,
  11790035.9375,
  11790047.881167987,
  11810012.727178657,
  11837939.86120326,
  11843647.166989485,
  11844701.240905834,
  11847115.946895441,
  11848274.991723655,
  11849527.65492637,
  11849553.125,
  11863568.75,
  11864361.71980095,
  11872835.493446702,
  11874316.33757873,
  11877539.58277209,
  11881538.89014271,
  11882171.063570684,
  11883011.878383249,
  11883020.518726151,
  11885174.61399025,
  11885326.5625,
  11927726.522682112,
  11929181.029014822,
  11953144.276878664,
  11953200.867572479,
  11990762.895079818,
  12007396.399782972,
  12013021.875,
  12014618.332461325,
  12021932.3798308,
  12033897.563695157,
  12046117.653045557,
  12048065.247748006,
  12054915.494098885,
  12056543.117510285,
  12059475.848050015,
  12071532.811497629,
  12076233.748617668,
  12080926.145167224,
  12088715.491758378,
  12098971.114727372,
  12101109.375,
  12105437.273901032,
  12110031.25,
  12116352.202408178,
  12125789.0625,
  12190859.375,
  12197207.177399764,
  12200792.66132302,
  12211179.6875,
  12221448.330772126,
  12228583.993977634,
  12233535.760050531,
  12238300.043272212,
  12238580.25172933,
  12253145.238302752,
  12254297.326422153,
  12258410.678815756,
  12288957.164853279,
  12289310.838695673,
  12289330.253045611,
  12297710.244243065,
  12305519.72271929,
  12310327.398157166,
  12322989.0625,
  12324044.164532494,
  12329820.3125,
  12332577.09308798,
  12337011.81811188,
  12341743.810120482,
  12342501.695288636,
  12342810.497005822,
  12343754.962555766,
  12343804.782191912,
  12352654.674848929,
  12352754.700761702,
  12354356.25,
  12357963.68856273,
  12362640.024070436,
  12369662.946120147,
  12370490.625,
  12371351.5625,
  12372298.410756001,
  12372321.301103534,
  12377263.815551868,
  12381496.984201385,
  12402448.30154804,
  12403637.384970969,
  12409410.643863685,
  12413746.302096846,
  12420857.39869064,
  12428475.024861988,
  12431436.98469129,
  12432332.8125,
  12438610.9375,
  12441056.185473407,
  12449480.577314021,
  12452387.5,
  12513439.0625,
  12521259.519716527,
  12523334.280990949,
  12552088.706217892,
  12565009.510186302,
  12566767.180312986,
  12567655.284750745,
  12569330.369037924,
  12594751.494608602,
  12599437.419326715,
  12605604.046514783,
  12625307.082359781,
  12651582.0887859,
  12664801.392066179,
  12682617.881817695,
  12706485.9375,
  12706814.447803903,
  12707018.344998086,
  12711548.4375,
  12724315.250892384,
  12729414.95502987,
  12730318.531745126,
  12730430.8073495,
  12741551.201999674,
  12749264.0625,
  12761198.4375,
  12766007.8125,
  12784679.6875,
  12790893.003129179,
  12799169.906802809,
  12799207.118685463,
  12799407.8125,
  12802906.25,
  12808428.96446223,
  12818371.875,
  12833260.9375,
  12848030.424400657,
  12853403.784758795,
  12853431.624389589,
  12868784.134136362,
  12877820.080221934,
  12880792.1875,
  12891910.87749689,
  12903998.306881944,
  12923870.334241977,
  12926803.017623859,
  12948012.5,
  12950003.125,
  12990785.986276973,
  12990840.18655162,
  12991378.785256276,
  12995511.925688231,
  13007307.8125,
  13022611.553608151,
  13039106.25,
  13052739.91108321,
  13052741.444961546,
  13067412.415981488,
  13079221.614540862,
  13099546.875,
  13108734.375,
  13109334.091298334,
  13114684.760340426,
  13135956.25,
  13166376.016174827,
  13237218.828458795,
  13238315.625,
  13240207.438232316,
  13284890.243905606,
  13345464.122766785,
  13347806.096652368,
  13358965.42547117,
  13451263.56546446,
  13453270.691184577,
  13461675.0,
  13477257.8125,
  13484950.0,
  13490188.990311217,
  13493904.641638903,
  13495427.723665878,
  13510346.875,
  13514675.331474775,
  13516314.152254121,
  13517257.611157248,
  13517991.908552757,
  13518212.5,
  13520275.979378957,
  13521857.8125,
  13522142.1875,
  13530309.375,
  13532932.037453566,
  13534326.232072355],
 [58.57129563365466,
  23.328665517742504,
  9.36726630655112,
  11.12229462231287,
  14.770376808832328,
  21.59350405221729,
  38.62831621957545,
  77.01988450447634,
  11.616668158650038,
  84.17367176446945,
  30.107573279735462,
  131.6150320839601,
  68.13887929835124,
  62.08080418791789,
  43.74876351245008,
  10.22275849910139,
  5.5312516523736,
  8.249341249518777,
  22.231193594763347,
  15.317049370541051,
  21.54121093407292,
  9.013301064113712,
  20.79998239349255,
  23.694199178673365,
  21.73083945786773,
  9.691076810242192,
  86.82646348004627,
  19.697940112795823,
  23.499044680942433,
  17.169735786988685,
  21.457553358555472,
  107.95347434073373,
  19.828060495144086,
  97.30898216552573,
  6.008093117683541,
  53.78811558150949,
  71.30625411477098,
  17.42578926336175,
  16.525051636394153,
  156.71906533101884,
  8.877242749269037,
  8.72821683525696,
  15.441242654584638,
  8.083370748914245,
  24.32974346333232,
  46.4480821123653,
  23.806962078445217,
  22.46456624981189,
  11.937984629123887,
  72.82971066940117,
  100.12379395068335,
  129.15656932524934,
  70.59139035957244,
  53.024744082246976,
  28.20820245736775,
  68.03861491336768,
  35.411318375397094,
  117.22114238344923,
  6.81480158563829,
  9.915921186435412,
  22.84058066284163,
  29.619554096111013,
  13.209474417145264,
  11.827571077267034,
  16.872197658289267,
  14.587326840690285,
  148.0938748630612,
  23.756885666380132,
  24.989186266702117,
  52.16638492411569,
  38.316952641698954,
  65.48612834509218,
  34.2748543385096,
  15.592457960109797,
  16.091839191889406,
  167.73628832260738,
  48.23302420932012,
  48.58294173998369,
  8.318012451740083,
  60.73228484457848,
  26.762608272901442,
  21.28631381524968,
  21.779050978050748,
  78.63091873995685,
  80.49417947643087,
  17.462517979085703,
  108.31644289537317,
  6.238157040858988,
  32.2703278753571,
  79.75699045791328,
  10.891831667939877,
  32.07951129523607,
  9.101952815995922,
  28.418149886433124,
  22.72810395713356,
  54.14937108586418,
  98.20596981584045,
  18.90434093050252,
  233.18438593839977,
  18.71682281403566,
  9.131622251901165,
  23.91949008104107,
  36.83368791552786,
  25.397582516500652,
  6.194608435214572,
  18.457356755809137,
  12.57417857976828,
  22.582416848868064,
  78.43480847732826,
  17.588622633306084,
  16.710473929955505,
  18.333106236324454,
  86.90799630936115,
  62.77681944800598,
  54.564269372991866,
  56.01944133432377,
  5.374731935814738,
  7.578622078901453,
  42.29528841899147,
  7.24634962864628,
  5.714352045277359,
  7.93963295936251,
  57.15611509588377,
  19.040792997409273,
  32.634211796866566,
  21.47843755477674,
  80.04600620140931,
  8.361298952187017,
  9.920709906400463,
  12.861796582918545,
  69.39299618892508,
  79.38038170434771,
  20.92889574020502,
  137.91747645091036,
  11.610457113291911,
  92.83697778058237,
  39.715756851215275,
  60.51001373835122,
  31.0855005571628,
  23.02139411775415,
  53.0637708230102,
  30.98964325775669,
  15.907995240216321,
  15.857485654583426,
  16.728308241688755,
  5.564303596171225,
  27.191473078258888,
  67.74682316356235,
  31.003666627545083,
  70.2291432428266,
  8.594517134429749,
  46.52990160610279,
  21.222596076104963,
  138.78863745796434,
  61.33617561125287,
  42.07395597781649,
  13.763574370901136,
  48.440778790852136,
  67.25645270511008,
  27.180690309905597,
  61.95360624180139,
  71.36649074128283,
  89.60403856537958,
  6.443113151108977,
  69.26564548931842,
  29.04465649530925,
  51.729288263465754,
  5.816443868104516,
  68.71759332290607,
  18.43620385351919,
  91.60040578478642,
  14.461280941688948,
  26.906858095035457,
  34.99723423917127,
  103.4796809426892,
  5.30784073166951,
  43.386124028612095,
  27.846708942980143,
  23.327517615929494,
  15.598648744887045,
  31.94399856468777,
  15.831728207444906,
  43.260646362737,
  10.139488701670984,
  37.258384513326874,
  42.52023423473864,
  13.062437161671866,
  9.074396730136435,
  80.3124751540522,
  69.775678350818,
  38.86698832746866,
  16.107263380105454,
  14.783698488677725,
  68.43603223327298,
  7.7399636915894,
  7.419129481939966,
  64.54525205567396,
  5.041133182690351,
  5.330204164668274,
  73.0561103338258,
  17.54441027291371,
  76.40916407082887,
  16.566690242784773,
  14.983521269682699,
  49.92082343531145,
  25.051940160910558,
  17.812736683158946,
  42.3203996633591,
  29.95656116168032,
  99.79834681844255,
  29.37950960717422,
  10.226445361361211,
  28.125321407806968,
  65.48826355593854,
  50.88249855363242,
  16.260804840900246,
  59.155422910980285,
  35.31527582744585,
  42.65180576982577,
  27.839812483869153,
  54.703754987110734,
  50.87066539301161,
  32.41372900985752,
  51.81477305323438,
  144.14270219903653,
  14.650642783225805,
  57.021905939918156,
  42.10175670835298,
  26.069872330926994,
  22.42421619960278,
  50.94666920123797,
  66.76695875207633,
  13.263261617830537,
  6.5514907396206565,
  5.359558144839126,
  77.28858932733307,
  81.43332647620434,
  26.867882672116494,
  126.00265525531108,
  8.351013104914326,
  88.19882572811005,
  101.8347967199233,
  72.72251029791693,
  36.4384575571025,
  56.591295044543656,
  40.27546416743414,
  10.823037335478448,
  18.381296649966,
  19.32987228975821,
  60.3936212901804,
  94.72190427933973,
  51.40895065123914,
  5.55431835679427,
  11.855975967339415,
  18.153972827160825,
  18.552050411019373,
  39.7182584843535,
  40.049271077718515,
  22.47539430654753,
  182.69356753659085,
  5.509846518677861,
  28.84135166895968,
  51.12545609546079,
  10.767581495949146,
  79.4841194221117,
  45.75603200191485,
  46.34940941454821,
  73.34973500824889,
  10.769812442546373,
  75.27228548640468,
  16.036874973093568,
  59.802772766331856,
  7.694964239803823,
  12.9407637201959,
  51.572087540175154,
  32.71859597754264,
  13.212258322110696,
  24.38670023020801,
  30.124652957269248,
  86.28483780167355,
  72.9400169310025,
  20.79423998308172,
  23.76429993496005,
  50.29315939663727,
  6.07658556545266,
  18.157991048637694,
  18.07517983009098,
  10.050859097506073,
  17.525070344403257,
  5.377812779467981,
  23.636555277077612,
  7.7057272640520535,
  78.68841629943648,
  20.16373155490395,
  5.146281443664074,
  61.912316977308265,
  12.00699428808731,
  24.67238943527335,
  99.29041396605,
  18.09336445328164,
  48.571636881385416,
  31.45912946661526,
  8.990295061406366,
  103.97809716307002,
  50.16633216897899,
  7.8640030826939835,
  12.149744617196212,
  7.904661609394422,
  13.238449143753693,
  91.2459983296458,
  13.265494400395282,
  42.63067409745607,
  6.566492661525614,
  88.83926493228911,
  42.48907152462587,
  81.66792630073043,
  23.282256627242273,
  7.331131452739423,
  59.09539057205833,
  19.950924659928237,
  77.79831581620144,
  7.722787896633274,
  40.52139076002147,
  25.929206594139185,
  16.69588680799624,
  81.06882832560599,
  21.893964156425387,
  5.312424091433347,
  61.822419910812826,
  17.39840254945882,
  60.46415762647308,
  13.340928216794993,
  141.18143884851554,
  17.018811193889686,
  59.02258412203743,
  116.39418952261948,
  90.94118772227516,
  70.59337289782707,
  132.56462889682385,
  14.626173391873497,
  50.29797505234686,
  28.57552352958628,
  10.519327166887544,
  37.01713326461523,
  30.75026396146229,
  14.032134907702076,
  20.425479256665515,
  50.40822569690259,
  6.328279744220919,
  13.149106538650674,
  6.450728254195457,
  5.49545863678443,
  18.656630670812575,
  40.99519087454638,
  15.200781981343658,
  14.317915559037155,
  24.868616030177666,
  10.157382009763202,
  7.326284635079966,
  13.962481840405337,
  20.012385652509902,
  5.751101387378411,
  26.20806013452035,
  34.59345917226382,
  17.5235764804364,
  20.98968102498617,
  9.056752559902506,
  175.68096518954238,
  11.359136541810196,
  15.860227939307938,
  9.105044676174387,
  15.151655117490192,
  62.32496375039588,
  24.18798496487939,
  86.60167694538129,
  80.68024762671142,
  67.97514803533488,
  7.00404445104064,
  13.492083552928797,
  24.66546823235344,
  5.1110381785127235,
  31.086559576735475,
  8.134193960647435,
  12.261396551254059,
  179.16369037630952,
  38.77790099737566,
  24.933410800817235,
  11.68303703572287,
  9.527111942095727,
  28.800696847361685,
  67.45924510874437,
  112.76853617011083,
  55.254770725096094,
  8.59952766459498,
  60.48536942173719,
  14.086083989183908,
  76.44651225158047,
  73.2761537423953,
  196.1114856797066,
  89.6537487409889,
  88.0923069077297,
  117.32547584660453,
  7.309232664133151,
  28.780179569480687,
  38.25224693036286,
  41.94458145944858,
  11.291360973015614,
  6.461562052628999,
  24.6917061140915,
  7.2242124767958185,
  10.927257236813185,
  47.62362388351568,
  183.85233359153742,
  32.579774661073756,
  64.77866611519207,
  10.838996283655824,
  41.169472353994664,
  10.369449001081755,
  56.92070214626546,
  6.982504920670008,
  15.884393659436185,
  86.06267732080241,
  5.474440189745283,
  23.47877023620519,
  94.55988509164382,
  14.910057425333306,
  7.003133155046004,
  89.6337161186774,
  75.26944018128385,
  19.089546876209987,
  58.44417850116866,
  17.9832987985655,
  54.310640335551035,
  25.530683750967224,
  125.6281123274559,
  20.36003032389581,
  9.512086670118443,
  17.386667160663894,
  135.85215721557338,
  14.791883535850683,
  71.79947226306216,
  25.192155261843396,
  20.283515577844902,
  70.61856220748568,
  6.6722021539051735,
  37.21439795166085,
  39.96894729959123,
  8.990523694002215,
  6.1627703443934285,
  56.81213546390012,
  23.358726498640763,
  7.556569321964341,
  20.265975355084027,
  8.792474473951097,
  26.385622965167485,
  92.85036261086564,
  43.452585673475284,
  12.816397223053084,
  19.23800909366875,
  31.621045576972197,
  63.525278237650426,
  43.28584199422976,
  12.300678089842004,
  42.72075730552848,
  116.84460293206509,
  25.53864554984964,
  9.766984687656507,
  49.37138930022858,
  45.176521371135365,
  65.8531420050303,
  6.245036241017265,
  15.925733943875045,
  8.609329995630876,
  20.695347864925772,
  28.32518865582491,
  6.447290640514619,
  12.939082660556197,
  26.712849786981838,
  31.215084883327602,
  7.95110559306821,
  20.641922700487836,
  28.365839962095663,
  77.63607194104131,
  5.976337210244874,
  38.35073950444072,
  26.382361538335587,
  10.206453831429773,
  31.376326235489408,
  11.3418073405289,
  55.8422275790988,
  133.5400734161776,
  33.09038242349268,
  180.12528554154386,
  14.992298895427966,
  75.67504340030679,
  10.589535888279855,
  111.67574095779477,
  15.22930633244243,
  119.76377374978344,
  60.58671191836012,
  52.538773934390456,
  20.7863115040151,
  42.56229480817096,
  15.374737044093054,
  95.77710169576866,
  127.69122158678363,
  95.94470632849146,
  20.795425774432267,
  16.10995609864871,
  5.941573585537059,
  19.21111454187715,
  39.609573071167645,
  34.07889488678319,
  74.77133389595116,
  7.960237648452039,
  42.22765535669677,
  86.78925805911857,
  5.295578478558962,
  45.4946385886513,
  111.65075026848433,
  15.807871731274638,
  71.10647599244821,
  25.796377837171956,
  5.491799561378873,
  27.24999653394341,
  9.264108361484686,
  10.25461058113253,
  8.781940846114253,
  24.403069070465996,
  25.56000560421904,
  12.419381338979749,
  107.3726571542254,
  70.06755545484957,
  22.49616117626991,
  103.59667879436307,
  5.329263754375793,
  49.81147627110015,
  38.684270314870076,
  10.499787209669593,
  13.725060745174016,
  57.077735427687706,
  30.675260319551306,
  68.57777226656908,
  29.339771056622904,
  19.339551561044335,
  9.010408070409003,
  10.192559568245219,
  86.68882818776349,
  37.52342577134934,
  17.937975413383175,
  68.64176254882415,
  58.402195942291335,
  61.007124080438274,
  20.08474691225381,
  38.12336010072819,
  8.16318340389225,
  10.783065820229,
  59.217371572387265,
  34.402161605715136,
  15.550415113744474,
  78.06052672145165,
  14.547449226693427,
  58.9723803443803,
  42.5418929792828,
  45.59566316994359,
  6.4204890279454245,
  17.16039493547367,
  20.78250383230657,
  46.72032240025519,
  17.957708277714797,
  72.29367086727052,
  62.54387065716355,
  12.707495638153997,
  14.630656905864686,
  24.269312420828086,
  67.79129075254171,
  91.66759331944496,
  131.746227286324,
  50.32337164519672,
  51.37164059874975,
  8.017361132634237,
  75.58998345415668,
  52.676468565988024,
  22.84627108128996,
  13.663940877592156,
  15.17630386222057,
  101.21006935474807,
  6.724315397593212,
  13.3341723286413,
  14.69499562292919,
  60.461644636278635,
  45.253975109568394,
  47.63846298454778,
  53.35919978459143,
  13.479818128479813,
  10.55344042729989,
  77.45686375400152,
  16.375809228978582,
  9.408466211584761,
  66.4613206584992,
  64.94149757647627,
  70.6780238001515,
  82.74955534392683,
  73.85440796874556,
  58.33454333398977,
  44.20099986827395,
  10.09964640222436,
  17.344598528906207])

Energy spectrum

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

([7327426.165155034,
  7799377.79122286,
  7799390.667111297,
  7939804.684034693,
  7956131.510062348,
  8109033.885006448,
  8209116.604035509,
  8216208.976214902,
  8335283.104312408,
  8486787.5,
  8556695.013555547,
  8608280.537665417,
  8789835.9375,
  8830321.875,
  8893587.5,
  8893589.957035514,
  8909754.64975903,
  8921381.17891598,
  8922395.911433162,
  8968225.095157664,
  8968268.254809383,
  8969201.160707876,
  8970016.615093034,
  8979272.431575812,
  8979939.62047793,
  8993147.271797227,
  9007339.0625,
  9009068.049013903,
  9012516.867536742,
  9033400.48309831,
  9055802.4091866,
  9095473.4375,
  9095495.205135133,
  9106713.365973506,
  9131542.94382511,
  9143060.9375,
  9154457.631342614,
  9154467.82935047,
  9158136.056151805,
  9166148.13453587,
  9200683.784165535,
  9217205.438665004,
  9217976.397194602,
  9220791.782393917,
  9230026.246870363,
  9246042.1875,
  9251360.641888585,
  9278027.71716614,
  9285183.01007829,
  9296638.925568322,
  9299884.049148375,
  9345634.234690037,
  9346292.1875,
  9346492.1875,
  9347976.003853073,
  9348831.25,
  9348848.4375,
  9349648.16009287,
  9352925.251617005,
  9353805.306362716,
  9400261.237052396,
  9401669.840147445,
  9422966.644377237,
  9423806.690948995,
  9432250.182095895,
  9432695.216365414,
  9433114.0625,
  9443974.70048084,
  9446769.984196296,
  9450565.625,
  9453778.125,
  9455848.4375,
  9460913.985221656,
  9463617.348407395,
  9466659.647269633,
  9483931.777406503,
  9508966.540260786,
  9514589.957339026,
  9514901.62147496,
  9517695.3125,
  9521261.893517284,
  9522652.604017746,
  9523066.422492987,
  9523903.125,
  9524551.5625,
  9524792.721884431,
  9525302.290748097,
  9528440.391628487,
  9552768.36184259,
  9559687.178782301,
  9559999.8436646,
  9560656.25,
  9561341.903074494,
  9579026.212965334,
  9589400.645778496,
  9595899.504652567,
  9599269.238749076,
  9605294.50687347,
  9609730.7235725,
  9610199.167961268,
  9613146.628884837,
  9614226.275596865,
  9618464.663514046,
  9664529.894618401,
  9664663.465508314,
  9665144.442758262,
  9668211.700415751,
  9700166.790179495,
  9723083.913256293,
  9730400.555775732,
  9736848.83526096,
  9741232.42610989,
  9742377.981699638,
  9742823.4375,
  9765887.035624983,
  9766693.75,
  9766719.67725108,
  9770410.868924772,
  9780678.125,
  9781144.204739803,
  9785720.431165475,
  9788358.655751918,
  9808085.9375,
  9808583.588169497,
  9813212.39343075,
  9815598.193608623,
  9824971.875,
  9826076.701176478,
  9831420.93589777,
  9835396.19625881,
  9835639.0625,
  9835970.3125,
  9836261.35081455,
  9836515.450608224,
  9837967.475851655,
  9838809.375,
  9856243.75,
  9856481.048540115,
  9885837.5,
  9887371.517387621,
  9895981.25,
  9898386.723926332,
  9900897.107390271,
  9901403.346210882,
  9901406.651192874,
  9907061.564476078,
  9912917.60361787,
  9916457.8125,
  9916729.719260572,
  9944321.875,
  9945217.727706514,
  9945235.571232207,
  9950724.472448848,
  9953787.713318588,
  9964235.9375,
  9964544.201671857,
  9965747.700305002,
  9966795.3125,
  9968779.6875,
  9969680.755824512,
  9970300.0,
  9974410.9375,
  9974676.19605752,
  9986223.151097342,
  9990026.5625,
  9990045.3125,
  9999723.843523532,
  10003210.970740348,
  10007898.4375,
  10007951.926766513,
  10012539.194053937,
  10017620.1805305,
  10022336.710041035,
  10025506.25,
  10026787.009956984,
  10027405.22051793,
  10030017.1875,
  10034018.71479335,
  10043229.35362895,
  10064962.6623484,
  10064971.687430203,
  10071823.652972111,
  10078284.375,
  10080052.385105353,
  10081281.503864676,
  10086465.625,
  10089349.386021024,
  10090369.992063241,
  10095773.4375,
  10096840.183617594,
  10097425.0,
  10107094.355760887,
  10126921.832971035,
  10127321.875,
  10142575.362841457,
  10149932.273067508,
  10160179.6875,
  10180140.60061183,
  10180251.216514101,
  10183453.125,
  10190345.236790597,
  10201112.5,
  10223676.932226567,
  10224796.885978809,
  10236476.5625,
  10248871.140190983,
  10252125.59280453,
  10255750.0,
  10257955.108714499,
  10270148.309936251,
  10279271.600184523,
  10292621.210173285,
  10294989.041992733,
  10295976.5625,
  10296107.141611105,
  10296147.18283717,
  10308536.8888846,
  10312792.1875,
  10322008.61804829,
  10325643.714425078,
  10326206.25,
  10358231.25,
  10360787.5,
  10364907.8125,
  10370717.1875,
  10380779.60903569,
  10387392.1875,
  10387644.272717165,
  10388113.212121796,
  10389802.23951119,
  10410451.5625,
  10416700.0,
  10422747.914414411,
  10424831.648875559,
  10425022.380816001,
  10426523.468860503,
  10427165.625,
  10431043.770637976,
  10432017.376520744,
  10440187.61989328,
  10445104.6875,
  10445145.3125,
  10448706.25,
  10454005.239069661,
  10455601.5625,
  10455634.375,
  10458771.670926899,
  10458947.556441935,
  10458978.939359512,
  10461587.5,
  10461907.43368539,
  10462978.125,
  10462995.742463369,
  10465705.108299486,
  10465730.334377792,
  10474415.157231178,
  10477700.0,
  10478309.519962726,
  10479475.971809018,
  10481060.883694548,
  10487465.136796782,
  10495787.160922995,
  10501945.3125,
  10503381.006701568,
  10508401.175898485,
  10508667.1875,
  10509104.6875,
  10513143.679788072,
  10518677.883292628,
  10529090.275011314,
  10529220.572978143,
  10529242.73658636,
  10531850.246684035,
  10531909.937488614,
  10534096.875,
  10538176.38933034,
  10539913.423931781,
  10542124.748207744,
  10544145.3125,
  10550090.476358475,
  10565461.025291791,
  10566777.551967675,
  10567358.185831482,
  10568051.5625,
  10580235.826049905,
  10585429.060469655,
  10591216.259630535,
  10595089.806242248,
  10596014.607546404,
  10598147.886318216,
  10603913.213103363,
  10607828.399749015,
  10620724.530082503,
  10621599.355544081,
  10621700.723934634,
  10625354.53769265,
  10626814.546452133,
  10642616.993597535,
  10642817.068049252,
  10650596.99334149,
  10655201.5625,
  10656549.276483875,
  10662374.567019334,
  10663439.084789084,
  10665101.091205105,
  10667427.230182495,
  10667959.943337955,
  10668053.887993434,
  10672405.717701185,
  10683301.5625,
  10686336.977740964,
  10696645.3125,
  10698866.11348568,
  10708080.354843566,
  10708385.9375,
  10711884.375,
  10731785.761690062,
  10734933.554308563,
  10736203.125,
  10742700.892635493,
  10745604.6875,
  10746926.881369725,
  10748553.125,
  10751181.103338206,
  10751209.941040728,
  10760750.66745447,
  10767169.851767238,
  10767230.363151476,
  10785190.625,
  10787125.523350062,
  10787492.1875,
  10805005.962528912,
  10805239.0625,
  10806127.018236192,
  10807394.87982792,
  10810114.637627069,
  10812429.6875,
  10815804.6875,
  10822979.66390588,
  10857645.227232529,
  10863868.75,
  10918338.336707894,
  10919619.884979524,
  10945764.0625,
  10964353.57477003,
  10994955.989037808,
  10998174.385874284,
  10998653.125,
  11002853.593632428,
  11006450.996597279,
  11016062.26404281,
  11018686.699659476,
  11025283.719425179,
  11032134.266904904,
  11053821.941399861,
  11095335.763671527,
  11125766.250932002,
  11132193.025735555,
  11150030.413336134,
  11155254.404516628,
  11164391.960484443,
  11166353.799916288,
  11176426.433260964,
  11185793.75,
  11232250.489999108,
  11235793.198009068,
  11236318.40628443,
  11237661.520932302,
  11255545.455682302,
  11255893.971820354,
  11263565.344589077,
  11264400.054062905,
  11274832.30043398,
  11278684.95155125,
  11280390.625,
  11285168.75,
  11287679.6875,
  11312735.47120625,
  11324629.247025544,
  11324632.065573448,
  11328610.003886038,
  11329234.145130951,
  11337037.200278776,
  11357041.246543927,
  11375579.13056528,
  11375700.0,
  11403106.99226797,
  11405376.890927367,
  11413892.410228908,
  11415845.069434233,
  11416109.375,
  11442310.9375,
  11473715.625,
  11478972.423351046,
  11484207.64891236,
  11516124.755877858,
  11522331.25,
  11536606.746160336,
  11560892.797188083,
  11583493.392771447,
  11589855.767051289,
  11617951.5625,
  11620610.044228334,
  11620613.713904215,
  11623857.8125,
  11623964.0625,
  11625526.664737193,
  11625709.210187174,
  11656000.560084011,
  11657281.627008902,
  11657748.305358363,
  11668675.143138427,
  11669969.509165164,
  11739986.777641645,
  11745700.0,
  11759813.869591098,
  11759853.078272732,
  11768611.953919861,
  11790035.9375,
  11790047.881167987,
  11810012.727178657,
  11837939.86120326,
  11843647.166989485,
  11844701.240905834,
  11847115.946895441,
  11848274.991723655,
  11849527.65492637,
  11849553.125,
  11863568.75,
  11864361.71980095,
  11872835.493446702,
  11874316.33757873,
  11877539.58277209,
  11881538.89014271,
  11882171.063570684,
  11883011.878383249,
  11883020.518726151,
  11885174.61399025,
  11885326.5625,
  11927726.522682112,
  11929181.029014822,
  11953144.276878664,
  11953200.867572479,
  11990762.895079818,
  12007396.399782972,
  12013021.875,
  12014618.332461325,
  12021932.3798308,
  12033897.563695157,
  12046117.653045557,
  12048065.247748006,
  12054915.494098885,
  12056543.117510285,
  12059475.848050015,
  12071532.811497629,
  12076233.748617668,
  12080926.145167224,
  12088715.491758378,
  12098971.114727372,
  12101109.375,
  12105437.273901032,
  12110031.25,
  12116352.202408178,
  12125789.0625,
  12190859.375,
  12197207.177399764,
  12200792.66132302,
  12211179.6875,
  12221448.330772126,
  12228583.993977634,
  12233535.760050531,
  12238300.043272212,
  12238580.25172933,
  12253145.238302752,
  12254297.326422153,
  12258410.678815756,
  12288957.164853279,
  12289310.838695673,
  12289330.253045611,
  12297710.244243065,
  12305519.72271929,
  12310327.398157166,
  12322989.0625,
  12324044.164532494,
  12329820.3125,
  12332577.09308798,
  12337011.81811188,
  12341743.810120482,
  12342501.695288636,
  12342810.497005822,
  12343754.962555766,
  12343804.782191912,
  12352654.674848929,
  12352754.700761702,
  12354356.25,
  12357963.68856273,
  12362640.024070436,
  12369662.946120147,
  12370490.625,
  12371351.5625,
  12372298.410756001,
  12372321.301103534,
  12377263.815551868,
  12381496.984201385,
  12402448.30154804,
  12403637.384970969,
  12409410.643863685,
  12413746.302096846,
  12420857.39869064,
  12428475.024861988,
  12431436.98469129,
  12432332.8125,
  12438610.9375,
  12441056.185473407,
  12449480.577314021,
  12452387.5,
  12513439.0625,
  12521259.519716527,
  12523334.280990949,
  12552088.706217892,
  12565009.510186302,
  12566767.180312986,
  12567655.284750745,
  12569330.369037924,
  12594751.494608602,
  12599437.419326715,
  12605604.046514783,
  12625307.082359781,
  12651582.0887859,
  12664801.392066179,
  12682617.881817695,
  12706485.9375,
  12706814.447803903,
  12707018.344998086,
  12711548.4375,
  12724315.250892384,
  12729414.95502987,
  12730318.531745126,
  12730430.8073495,
  12741551.201999674,
  12749264.0625,
  12761198.4375,
  12766007.8125,
  12784679.6875,
  12790893.003129179,
  12799169.906802809,
  12799207.118685463,
  12799407.8125,
  12802906.25,
  12808428.96446223,
  12818371.875,
  12833260.9375,
  12848030.424400657,
  12853403.784758795,
  12853431.624389589,
  12868784.134136362,
  12877820.080221934,
  12880792.1875,
  12891910.87749689,
  12903998.306881944,
  12923870.334241977,
  12926803.017623859,
  12948012.5,
  12950003.125,
  12990785.986276973,
  12990840.18655162,
  12991378.785256276,
  12995511.925688231,
  13007307.8125,
  13022611.553608151,
  13039106.25,
  13052739.91108321,
  13052741.444961546,
  13067412.415981488,
  13079221.614540862,
  13099546.875,
  13108734.375,
  13109334.091298334,
  13114684.760340426,
  13135956.25,
  13166376.016174827,
  13237218.828458795,
  13238315.625,
  13240207.438232316,
  13284890.243905606,
  13345464.122766785,
  13347806.096652368,
  13358965.42547117,
  13451263.56546446,
  13453270.691184577,
  13461675.0,
  13477257.8125,
  13484950.0,
  13490188.990311217,
  13493904.641638903,
  13495427.723665878,
  13510346.875,
  13514675.331474775,
  13516314.152254121,
  13517257.611157248,
  13517991.908552757,
  13518212.5,
  13520275.979378957,
  13521857.8125,
  13522142.1875,
  13530309.375,
  13532932.037453566,
  13534326.232072355],
 [58.57129563365466,
  23.328665517742504,
  9.36726630655112,
  11.12229462231287,
  14.770376808832328,
  21.59350405221729,
  38.62831621957545,
  77.01988450447634,
  11.616668158650038,
  84.17367176446945,
  30.107573279735462,
  131.6150320839601,
  68.13887929835124,
  62.08080418791789,
  43.74876351245008,
  10.22275849910139,
  5.5312516523736,
  8.249341249518777,
  22.231193594763347,
  15.317049370541051,
  21.54121093407292,
  9.013301064113712,
  20.79998239349255,
  23.694199178673365,
  21.73083945786773,
  9.691076810242192,
  86.82646348004627,
  19.697940112795823,
  23.499044680942433,
  17.169735786988685,
  21.457553358555472,
  107.95347434073373,
  19.828060495144086,
  97.30898216552573,
  6.008093117683541,
  53.78811558150949,
  71.30625411477098,
  17.42578926336175,
  16.525051636394153,
  156.71906533101884,
  8.877242749269037,
  8.72821683525696,
  15.441242654584638,
  8.083370748914245,
  24.32974346333232,
  46.4480821123653,
  23.806962078445217,
  22.46456624981189,
  11.937984629123887,
  72.82971066940117,
  100.12379395068335,
  129.15656932524934,
  70.59139035957244,
  53.024744082246976,
  28.20820245736775,
  68.03861491336768,
  35.411318375397094,
  117.22114238344923,
  6.81480158563829,
  9.915921186435412,
  22.84058066284163,
  29.619554096111013,
  13.209474417145264,
  11.827571077267034,
  16.872197658289267,
  14.587326840690285,
  148.0938748630612,
  23.756885666380132,
  24.989186266702117,
  52.16638492411569,
  38.316952641698954,
  65.48612834509218,
  34.2748543385096,
  15.592457960109797,
  16.091839191889406,
  167.73628832260738,
  48.23302420932012,
  48.58294173998369,
  8.318012451740083,
  60.73228484457848,
  26.762608272901442,
  21.28631381524968,
  21.779050978050748,
  78.63091873995685,
  80.49417947643087,
  17.462517979085703,
  108.31644289537317,
  6.238157040858988,
  32.2703278753571,
  79.75699045791328,
  10.891831667939877,
  32.07951129523607,
  9.101952815995922,
  28.418149886433124,
  22.72810395713356,
  54.14937108586418,
  98.20596981584045,
  18.90434093050252,
  233.18438593839977,
  18.71682281403566,
  9.131622251901165,
  23.91949008104107,
  36.83368791552786,
  25.397582516500652,
  6.194608435214572,
  18.457356755809137,
  12.57417857976828,
  22.582416848868064,
  78.43480847732826,
  17.588622633306084,
  16.710473929955505,
  18.333106236324454,
  86.90799630936115,
  62.77681944800598,
  54.564269372991866,
  56.01944133432377,
  5.374731935814738,
  7.578622078901453,
  42.29528841899147,
  7.24634962864628,
  5.714352045277359,
  7.93963295936251,
  57.15611509588377,
  19.040792997409273,
  32.634211796866566,
  21.47843755477674,
  80.04600620140931,
  8.361298952187017,
  9.920709906400463,
  12.861796582918545,
  69.39299618892508,
  79.38038170434771,
  20.92889574020502,
  137.91747645091036,
  11.610457113291911,
  92.83697778058237,
  39.715756851215275,
  60.51001373835122,
  31.0855005571628,
  23.02139411775415,
  53.0637708230102,
  30.98964325775669,
  15.907995240216321,
  15.857485654583426,
  16.728308241688755,
  5.564303596171225,
  27.191473078258888,
  67.74682316356235,
  31.003666627545083,
  70.2291432428266,
  8.594517134429749,
  46.52990160610279,
  21.222596076104963,
  138.78863745796434,
  61.33617561125287,
  42.07395597781649,
  13.763574370901136,
  48.440778790852136,
  67.25645270511008,
  27.180690309905597,
  61.95360624180139,
  71.36649074128283,
  89.60403856537958,
  6.443113151108977,
  69.26564548931842,
  29.04465649530925,
  51.729288263465754,
  5.816443868104516,
  68.71759332290607,
  18.43620385351919,
  91.60040578478642,
  14.461280941688948,
  26.906858095035457,
  34.99723423917127,
  103.4796809426892,
  5.30784073166951,
  43.386124028612095,
  27.846708942980143,
  23.327517615929494,
  15.598648744887045,
  31.94399856468777,
  15.831728207444906,
  43.260646362737,
  10.139488701670984,
  37.258384513326874,
  42.52023423473864,
  13.062437161671866,
  9.074396730136435,
  80.3124751540522,
  69.775678350818,
  38.86698832746866,
  16.107263380105454,
  14.783698488677725,
  68.43603223327298,
  7.7399636915894,
  7.419129481939966,
  64.54525205567396,
  5.041133182690351,
  5.330204164668274,
  73.0561103338258,
  17.54441027291371,
  76.40916407082887,
  16.566690242784773,
  14.983521269682699,
  49.92082343531145,
  25.051940160910558,
  17.812736683158946,
  42.3203996633591,
  29.95656116168032,
  99.79834681844255,
  29.37950960717422,
  10.226445361361211,
  28.125321407806968,
  65.48826355593854,
  50.88249855363242,
  16.260804840900246,
  59.155422910980285,
  35.31527582744585,
  42.65180576982577,
  27.839812483869153,
  54.703754987110734,
  50.87066539301161,
  32.41372900985752,
  51.81477305323438,
  144.14270219903653,
  14.650642783225805,
  57.021905939918156,
  42.10175670835298,
  26.069872330926994,
  22.42421619960278,
  50.94666920123797,
  66.76695875207633,
  13.263261617830537,
  6.5514907396206565,
  5.359558144839126,
  77.28858932733307,
  81.43332647620434,
  26.867882672116494,
  126.00265525531108,
  8.351013104914326,
  88.19882572811005,
  101.8347967199233,
  72.72251029791693,
  36.4384575571025,
  56.591295044543656,
  40.27546416743414,
  10.823037335478448,
  18.381296649966,
  19.32987228975821,
  60.3936212901804,
  94.72190427933973,
  51.40895065123914,
  5.55431835679427,
  11.855975967339415,
  18.153972827160825,
  18.552050411019373,
  39.7182584843535,
  40.049271077718515,
  22.47539430654753,
  182.69356753659085,
  5.509846518677861,
  28.84135166895968,
  51.12545609546079,
  10.767581495949146,
  79.4841194221117,
  45.75603200191485,
  46.34940941454821,
  73.34973500824889,
  10.769812442546373,
  75.27228548640468,
  16.036874973093568,
  59.802772766331856,
  7.694964239803823,
  12.9407637201959,
  51.572087540175154,
  32.71859597754264,
  13.212258322110696,
  24.38670023020801,
  30.124652957269248,
  86.28483780167355,
  72.9400169310025,
  20.79423998308172,
  23.76429993496005,
  50.29315939663727,
  6.07658556545266,
  18.157991048637694,
  18.07517983009098,
  10.050859097506073,
  17.525070344403257,
  5.377812779467981,
  23.636555277077612,
  7.7057272640520535,
  78.68841629943648,
  20.16373155490395,
  5.146281443664074,
  61.912316977308265,
  12.00699428808731,
  24.67238943527335,
  99.29041396605,
  18.09336445328164,
  48.571636881385416,
  31.45912946661526,
  8.990295061406366,
  103.97809716307002,
  50.16633216897899,
  7.8640030826939835,
  12.149744617196212,
  7.904661609394422,
  13.238449143753693,
  91.2459983296458,
  13.265494400395282,
  42.63067409745607,
  6.566492661525614,
  88.83926493228911,
  42.48907152462587,
  81.66792630073043,
  23.282256627242273,
  7.331131452739423,
  59.09539057205833,
  19.950924659928237,
  77.79831581620144,
  7.722787896633274,
  40.52139076002147,
  25.929206594139185,
  16.69588680799624,
  81.06882832560599,
  21.893964156425387,
  5.312424091433347,
  61.822419910812826,
  17.39840254945882,
  60.46415762647308,
  13.340928216794993,
  141.18143884851554,
  17.018811193889686,
  59.02258412203743,
  116.39418952261948,
  90.94118772227516,
  70.59337289782707,
  132.56462889682385,
  14.626173391873497,
  50.29797505234686,
  28.57552352958628,
  10.519327166887544,
  37.01713326461523,
  30.75026396146229,
  14.032134907702076,
  20.425479256665515,
  50.40822569690259,
  6.328279744220919,
  13.149106538650674,
  6.450728254195457,
  5.49545863678443,
  18.656630670812575,
  40.99519087454638,
  15.200781981343658,
  14.317915559037155,
  24.868616030177666,
  10.157382009763202,
  7.326284635079966,
  13.962481840405337,
  20.012385652509902,
  5.751101387378411,
  26.20806013452035,
  34.59345917226382,
  17.5235764804364,
  20.98968102498617,
  9.056752559902506,
  175.68096518954238,
  11.359136541810196,
  15.860227939307938,
  9.105044676174387,
  15.151655117490192,
  62.32496375039588,
  24.18798496487939,
  86.60167694538129,
  80.68024762671142,
  67.97514803533488,
  7.00404445104064,
  13.492083552928797,
  24.66546823235344,
  5.1110381785127235,
  31.086559576735475,
  8.134193960647435,
  12.261396551254059,
  179.16369037630952,
  38.77790099737566,
  24.933410800817235,
  11.68303703572287,
  9.527111942095727,
  28.800696847361685,
  67.45924510874437,
  112.76853617011083,
  55.254770725096094,
  8.59952766459498,
  60.48536942173719,
  14.086083989183908,
  76.44651225158047,
  73.2761537423953,
  196.1114856797066,
  89.6537487409889,
  88.0923069077297,
  117.32547584660453,
  7.309232664133151,
  28.780179569480687,
  38.25224693036286,
  41.94458145944858,
  11.291360973015614,
  6.461562052628999,
  24.6917061140915,
  7.2242124767958185,
  10.927257236813185,
  47.62362388351568,
  183.85233359153742,
  32.579774661073756,
  64.77866611519207,
  10.838996283655824,
  41.169472353994664,
  10.369449001081755,
  56.92070214626546,
  6.982504920670008,
  15.884393659436185,
  86.06267732080241,
  5.474440189745283,
  23.47877023620519,
  94.55988509164382,
  14.910057425333306,
  7.003133155046004,
  89.6337161186774,
  75.26944018128385,
  19.089546876209987,
  58.44417850116866,
  17.9832987985655,
  54.310640335551035,
  25.530683750967224,
  125.6281123274559,
  20.36003032389581,
  9.512086670118443,
  17.386667160663894,
  135.85215721557338,
  14.791883535850683,
  71.79947226306216,
  25.192155261843396,
  20.283515577844902,
  70.61856220748568,
  6.6722021539051735,
  37.21439795166085,
  39.96894729959123,
  8.990523694002215,
  6.1627703443934285,
  56.81213546390012,
  23.358726498640763,
  7.556569321964341,
  20.265975355084027,
  8.792474473951097,
  26.385622965167485,
  92.85036261086564,
  43.452585673475284,
  12.816397223053084,
  19.23800909366875,
  31.621045576972197,
  63.525278237650426,
  43.28584199422976,
  12.300678089842004,
  42.72075730552848,
  116.84460293206509,
  25.53864554984964,
  9.766984687656507,
  49.37138930022858,
  45.176521371135365,
  65.8531420050303,
  6.245036241017265,
  15.925733943875045,
  8.609329995630876,
  20.695347864925772,
  28.32518865582491,
  6.447290640514619,
  12.939082660556197,
  26.712849786981838,
  31.215084883327602,
  7.95110559306821,
  20.641922700487836,
  28.365839962095663,
  77.63607194104131,
  5.976337210244874,
  38.35073950444072,
  26.382361538335587,
  10.206453831429773,
  31.376326235489408,
  11.3418073405289,
  55.8422275790988,
  133.5400734161776,
  33.09038242349268,
  180.12528554154386,
  14.992298895427966,
  75.67504340030679,
  10.589535888279855,
  111.67574095779477,
  15.22930633244243,
  119.76377374978344,
  60.58671191836012,
  52.538773934390456,
  20.7863115040151,
  42.56229480817096,
  15.374737044093054,
  95.77710169576866,
  127.69122158678363,
  95.94470632849146,
  20.795425774432267,
  16.10995609864871,
  5.941573585537059,
  19.21111454187715,
  39.609573071167645,
  34.07889488678319,
  74.77133389595116,
  7.960237648452039,
  42.22765535669677,
  86.78925805911857,
  5.295578478558962,
  45.4946385886513,
  111.65075026848433,
  15.807871731274638,
  71.10647599244821,
  25.796377837171956,
  5.491799561378873,
  27.24999653394341,
  9.264108361484686,
  10.25461058113253,
  8.781940846114253,
  24.403069070465996,
  25.56000560421904,
  12.419381338979749,
  107.3726571542254,
  70.06755545484957,
  22.49616117626991,
  103.59667879436307,
  5.329263754375793,
  49.81147627110015,
  38.684270314870076,
  10.499787209669593,
  13.725060745174016,
  57.077735427687706,
  30.675260319551306,
  68.57777226656908,
  29.339771056622904,
  19.339551561044335,
  9.010408070409003,
  10.192559568245219,
  86.68882818776349,
  37.52342577134934,
  17.937975413383175,
  68.64176254882415,
  58.402195942291335,
  61.007124080438274,
  20.08474691225381,
  38.12336010072819,
  8.16318340389225,
  10.783065820229,
  59.217371572387265,
  34.402161605715136,
  15.550415113744474,
  78.06052672145165,
  14.547449226693427,
  58.9723803443803,
  42.5418929792828,
  45.59566316994359,
  6.4204890279454245,
  17.16039493547367,
  20.78250383230657,
  46.72032240025519,
  17.957708277714797,
  72.29367086727052,
  62.54387065716355,
  12.707495638153997,
  14.630656905864686,
  24.269312420828086,
  67.79129075254171,
  91.66759331944496,
  131.746227286324,
  50.32337164519672,
  51.37164059874975,
  8.017361132634237,
  75.58998345415668,
  52.676468565988024,
  22.84627108128996,
  13.663940877592156,
  15.17630386222057,
  101.21006935474807,
  6.724315397593212,
  13.3341723286413,
  14.69499562292919,
  60.461644636278635,
  45.253975109568394,
  47.63846298454778,
  53.35919978459143,
  13.479818128479813,
  10.55344042729989,
  77.45686375400152,
  16.375809228978582,
  9.408466211584761,
  66.4613206584992,
  64.94149757647627,
  70.6780238001515,
  82.74955534392683,
  73.85440796874556,
  58.33454333398977,
  44.20099986827395,
  10.09964640222436,
  17.344598528906207])

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

([7327426.165155034,
  7799377.79122286,
  7799390.667111297,
  7939804.684034693,
  7956131.510062348,
  8109033.885006448,
  8209116.604035509,
  8216208.976214902,
  8335283.104312408,
  8486787.5,
  8556695.013555547,
  8608280.537665417,
  8789835.9375,
  8830321.875,
  8893587.5,
  8893589.957035514,
  8909754.64975903,
  8921381.17891598,
  8922395.911433162,
  8968225.095157664,
  8968268.254809383,
  8969201.160707876,
  8970016.615093034,
  8979272.431575812,
  8979939.62047793,
  8993147.271797227,
  9007339.0625,
  9009068.049013903,
  9012516.867536742,
  9033400.48309831,
  9055802.4091866,
  9095473.4375,
  9095495.205135133,
  9106713.365973506,
  9131542.94382511,
  9143060.9375,
  9154457.631342614,
  9154467.82935047,
  9158136.056151805,
  9166148.13453587,
  9200683.784165535,
  9217205.438665004,
  9217976.397194602,
  9220791.782393917,
  9230026.246870363,
  9246042.1875,
  9251360.641888585,
  9278027.71716614,
  9285183.01007829,
  9296638.925568322,
  9299884.049148375,
  9345634.234690037,
  9346292.1875,
  9346492.1875,
  9347976.003853073,
  9348831.25,
  9348848.4375,
  9349648.16009287,
  9352925.251617005,
  9353805.306362716,
  9400261.237052396,
  9401669.840147445,
  9422966.644377237,
  9423806.690948995,
  9432250.182095895,
  9432695.216365414,
  9433114.0625,
  9443974.70048084,
  9446769.984196296,
  9450565.625,
  9453778.125,
  9455848.4375,
  9460913.985221656,
  9463617.348407395,
  9466659.647269633,
  9483931.777406503,
  9508966.540260786,
  9514589.957339026,
  9514901.62147496,
  9517695.3125,
  9521261.893517284,
  9522652.604017746,
  9523066.422492987,
  9523903.125,
  9524551.5625,
  9524792.721884431,
  9525302.290748097,
  9528440.391628487,
  9552768.36184259,
  9559687.178782301,
  9559999.8436646,
  9560656.25,
  9561341.903074494,
  9579026.212965334,
  9589400.645778496,
  9595899.504652567,
  9599269.238749076,
  9605294.50687347,
  9609730.7235725,
  9610199.167961268,
  9613146.628884837,
  9614226.275596865,
  9618464.663514046,
  9664529.894618401,
  9664663.465508314,
  9665144.442758262,
  9668211.700415751,
  9700166.790179495,
  9723083.913256293,
  9730400.555775732,
  9736848.83526096,
  9741232.42610989,
  9742377.981699638,
  9742823.4375,
  9765887.035624983,
  9766693.75,
  9766719.67725108,
  9770410.868924772,
  9780678.125,
  9781144.204739803,
  9785720.431165475,
  9788358.655751918,
  9808085.9375,
  9808583.588169497,
  9813212.39343075,
  9815598.193608623,
  9824971.875,
  9826076.701176478,
  9831420.93589777,
  9835396.19625881,
  9835639.0625,
  9835970.3125,
  9836261.35081455,
  9836515.450608224,
  9837967.475851655,
  9838809.375,
  9856243.75,
  9856481.048540115,
  9885837.5,
  9887371.517387621,
  9895981.25,
  9898386.723926332,
  9900897.107390271,
  9901403.346210882,
  9901406.651192874,
  9907061.564476078,
  9912917.60361787,
  9916457.8125,
  9916729.719260572,
  9944321.875,
  9945217.727706514,
  9945235.571232207,
  9950724.472448848,
  9953787.713318588,
  9964235.9375,
  9964544.201671857,
  9965747.700305002,
  9966795.3125,
  9968779.6875,
  9969680.755824512,
  9970300.0,
  9974410.9375,
  9974676.19605752,
  9986223.151097342,
  9990026.5625,
  9990045.3125,
  9999723.843523532,
  10003210.970740348,
  10007898.4375,
  10007951.926766513,
  10012539.194053937,
  10017620.1805305,
  10022336.710041035,
  10025506.25,
  10026787.009956984,
  10027405.22051793,
  10030017.1875,
  10034018.71479335,
  10043229.35362895,
  10064962.6623484,
  10064971.687430203,
  10071823.652972111,
  10078284.375,
  10080052.385105353,
  10081281.503864676,
  10086465.625,
  10089349.386021024,
  10090369.992063241,
  10095773.4375,
  10096840.183617594,
  10097425.0,
  10107094.355760887,
  10126921.832971035,
  10127321.875,
  10142575.362841457,
  10149932.273067508,
  10160179.6875,
  10180140.60061183,
  10180251.216514101,
  10183453.125,
  10190345.236790597,
  10201112.5,
  10223676.932226567,
  10224796.885978809,
  10236476.5625,
  10248871.140190983,
  10252125.59280453,
  10255750.0,
  10257955.108714499,
  10270148.309936251,
  10279271.600184523,
  10292621.210173285,
  10294989.041992733,
  10295976.5625,
  10296107.141611105,
  10296147.18283717,
  10308536.8888846,
  10312792.1875,
  10322008.61804829,
  10325643.714425078,
  10326206.25,
  10358231.25,
  10360787.5,
  10364907.8125,
  10370717.1875,
  10380779.60903569,
  10387392.1875,
  10387644.272717165,
  10388113.212121796,
  10389802.23951119,
  10410451.5625,
  10416700.0,
  10422747.914414411,
  10424831.648875559,
  10425022.380816001,
  10426523.468860503,
  10427165.625,
  10431043.770637976,
  10432017.376520744,
  10440187.61989328,
  10445104.6875,
  10445145.3125,
  10448706.25,
  10454005.239069661,
  10455601.5625,
  10455634.375,
  10458771.670926899,
  10458947.556441935,
  10458978.939359512,
  10461587.5,
  10461907.43368539,
  10462978.125,
  10462995.742463369,
  10465705.108299486,
  10465730.334377792,
  10474415.157231178,
  10477700.0,
  10478309.519962726,
  10479475.971809018,
  10481060.883694548,
  10487465.136796782,
  10495787.160922995,
  10501945.3125,
  10503381.006701568,
  10508401.175898485,
  10508667.1875,
  10509104.6875,
  10513143.679788072,
  10518677.883292628,
  10529090.275011314,
  10529220.572978143,
  10529242.73658636,
  10531850.246684035,
  10531909.937488614,
  10534096.875,
  10538176.38933034,
  10539913.423931781,
  10542124.748207744,
  10544145.3125,
  10550090.476358475,
  10565461.025291791,
  10566777.551967675,
  10567358.185831482,
  10568051.5625,
  10580235.826049905,
  10585429.060469655,
  10591216.259630535,
  10595089.806242248,
  10596014.607546404,
  10598147.886318216,
  10603913.213103363,
  10607828.399749015,
  10620724.530082503,
  10621599.355544081,
  10621700.723934634,
  10625354.53769265,
  10626814.546452133,
  10642616.993597535,
  10642817.068049252,
  10650596.99334149,
  10655201.5625,
  10656549.276483875,
  10662374.567019334,
  10663439.084789084,
  10665101.091205105,
  10667427.230182495,
  10667959.943337955,
  10668053.887993434,
  10672405.717701185,
  10683301.5625,
  10686336.977740964,
  10696645.3125,
  10698866.11348568,
  10708080.354843566,
  10708385.9375,
  10711884.375,
  10731785.761690062,
  10734933.554308563,
  10736203.125,
  10742700.892635493,
  10745604.6875,
  10746926.881369725,
  10748553.125,
  10751181.103338206,
  10751209.941040728,
  10760750.66745447,
  10767169.851767238,
  10767230.363151476,
  10785190.625,
  10787125.523350062,
  10787492.1875,
  10805005.962528912,
  10805239.0625,
  10806127.018236192,
  10807394.87982792,
  10810114.637627069,
  10812429.6875,
  10815804.6875,
  10822979.66390588,
  10857645.227232529,
  10863868.75,
  10918338.336707894,
  10919619.884979524,
  10945764.0625,
  10964353.57477003,
  10994955.989037808,
  10998174.385874284,
  10998653.125,
  11002853.593632428,
  11006450.996597279,
  11016062.26404281,
  11018686.699659476,
  11025283.719425179,
  11032134.266904904,
  11053821.941399861,
  11095335.763671527,
  11125766.250932002,
  11132193.025735555,
  11150030.413336134,
  11155254.404516628,
  11164391.960484443,
  11166353.799916288,
  11176426.433260964,
  11185793.75,
  11232250.489999108,
  11235793.198009068,
  11236318.40628443,
  11237661.520932302,
  11255545.455682302,
  11255893.971820354,
  11263565.344589077,
  11264400.054062905,
  11274832.30043398,
  11278684.95155125,
  11280390.625,
  11285168.75,
  11287679.6875,
  11312735.47120625,
  11324629.247025544,
  11324632.065573448,
  11328610.003886038,
  11329234.145130951,
  11337037.200278776,
  11357041.246543927,
  11375579.13056528,
  11375700.0,
  11403106.99226797,
  11405376.890927367,
  11413892.410228908,
  11415845.069434233,
  11416109.375,
  11442310.9375,
  11473715.625,
  11478972.423351046,
  11484207.64891236,
  11516124.755877858,
  11522331.25,
  11536606.746160336,
  11560892.797188083,
  11583493.392771447,
  11589855.767051289,
  11617951.5625,
  11620610.044228334,
  11620613.713904215,
  11623857.8125,
  11623964.0625,
  11625526.664737193,
  11625709.210187174,
  11656000.560084011,
  11657281.627008902,
  11657748.305358363,
  11668675.143138427,
  11669969.509165164,
  11739986.777641645,
  11745700.0,
  11759813.869591098,
  11759853.078272732,
  11768611.953919861,
  11790035.9375,
  11790047.881167987,
  11810012.727178657,
  11837939.86120326,
  11843647.166989485,
  11844701.240905834,
  11847115.946895441,
  11848274.991723655,
  11849527.65492637,
  11849553.125,
  11863568.75,
  11864361.71980095,
  11872835.493446702,
  11874316.33757873,
  11877539.58277209,
  11881538.89014271,
  11882171.063570684,
  11883011.878383249,
  11883020.518726151,
  11885174.61399025,
  11885326.5625,
  11927726.522682112,
  11929181.029014822,
  11953144.276878664,
  11953200.867572479,
  11990762.895079818,
  12007396.399782972,
  12013021.875,
  12014618.332461325,
  12021932.3798308,
  12033897.563695157,
  12046117.653045557,
  12048065.247748006,
  12054915.494098885,
  12056543.117510285,
  12059475.848050015,
  12071532.811497629,
  12076233.748617668,
  12080926.145167224,
  12088715.491758378,
  12098971.114727372,
  12101109.375,
  12105437.273901032,
  12110031.25,
  12116352.202408178,
  12125789.0625,
  12190859.375,
  12197207.177399764,
  12200792.66132302,
  12211179.6875,
  12221448.330772126,
  12228583.993977634,
  12233535.760050531,
  12238300.043272212,
  12238580.25172933,
  12253145.238302752,
  12254297.326422153,
  12258410.678815756,
  12288957.164853279,
  12289310.838695673,
  12289330.253045611,
  12297710.244243065,
  12305519.72271929,
  12310327.398157166,
  12322989.0625,
  12324044.164532494,
  12329820.3125,
  12332577.09308798,
  12337011.81811188,
  12341743.810120482,
  12342501.695288636,
  12342810.497005822,
  12343754.962555766,
  12343804.782191912,
  12352654.674848929,
  12352754.700761702,
  12354356.25,
  12357963.68856273,
  12362640.024070436,
  12369662.946120147,
  12370490.625,
  12371351.5625,
  12372298.410756001,
  12372321.301103534,
  12377263.815551868,
  12381496.984201385,
  12402448.30154804,
  12403637.384970969,
  12409410.643863685,
  12413746.302096846,
  12420857.39869064,
  12428475.024861988,
  12431436.98469129,
  12432332.8125,
  12438610.9375,
  12441056.185473407,
  12449480.577314021,
  12452387.5,
  12513439.0625,
  12521259.519716527,
  12523334.280990949,
  12552088.706217892,
  12565009.510186302,
  12566767.180312986,
  12567655.284750745,
  12569330.369037924,
  12594751.494608602,
  12599437.419326715,
  12605604.046514783,
  12625307.082359781,
  12651582.0887859,
  12664801.392066179,
  12682617.881817695,
  12706485.9375,
  12706814.447803903,
  12707018.344998086,
  12711548.4375,
  12724315.250892384,
  12729414.95502987,
  12730318.531745126,
  12730430.8073495,
  12741551.201999674,
  12749264.0625,
  12761198.4375,
  12766007.8125,
  12784679.6875,
  12790893.003129179,
  12799169.906802809,
  12799207.118685463,
  12799407.8125,
  12802906.25,
  12808428.96446223,
  12818371.875,
  12833260.9375,
  12848030.424400657,
  12853403.784758795,
  12853431.624389589,
  12868784.134136362,
  12877820.080221934,
  12880792.1875,
  12891910.87749689,
  12903998.306881944,
  12923870.334241977,
  12926803.017623859,
  12948012.5,
  12950003.125,
  12990785.986276973,
  12990840.18655162,
  12991378.785256276,
  12995511.925688231,
  13007307.8125,
  13022611.553608151,
  13039106.25,
  13052739.91108321,
  13052741.444961546,
  13067412.415981488,
  13079221.614540862,
  13099546.875,
  13108734.375,
  13109334.091298334,
  13114684.760340426,
  13135956.25,
  13166376.016174827,
  13237218.828458795,
  13238315.625,
  13240207.438232316,
  13284890.243905606,
  13345464.122766785,
  13347806.096652368,
  13358965.42547117,
  13451263.56546446,
  13453270.691184577,
  13461675.0,
  13477257.8125,
  13484950.0,
  13490188.990311217,
  13493904.641638903,
  13495427.723665878,
  13510346.875,
  13514675.331474775,
  13516314.152254121,
  13517257.611157248,
  13517991.908552757,
  13518212.5,
  13520275.979378957,
  13521857.8125,
  13522142.1875,
  13530309.375,
  13532932.037453566,
  13534326.232072355],
 [58.57129563365466,
  23.328665517742504,
  9.36726630655112,
  11.12229462231287,
  14.770376808832328,
  21.59350405221729,
  38.62831621957545,
  77.01988450447634,
  11.616668158650038,
  84.17367176446945,
  30.107573279735462,
  131.6150320839601,
  68.13887929835124,
  62.08080418791789,
  43.74876351245008,
  10.22275849910139,
  5.5312516523736,
  8.249341249518777,
  22.231193594763347,
  15.317049370541051,
  21.54121093407292,
  9.013301064113712,
  20.79998239349255,
  23.694199178673365,
  21.73083945786773,
  9.691076810242192,
  86.82646348004627,
  19.697940112795823,
  23.499044680942433,
  17.169735786988685,
  21.457553358555472,
  107.95347434073373,
  19.828060495144086,
  97.30898216552573,
  6.008093117683541,
  53.78811558150949,
  71.30625411477098,
  17.42578926336175,
  16.525051636394153,
  156.71906533101884,
  8.877242749269037,
  8.72821683525696,
  15.441242654584638,
  8.083370748914245,
  24.32974346333232,
  46.4480821123653,
  23.806962078445217,
  22.46456624981189,
  11.937984629123887,
  72.82971066940117,
  100.12379395068335,
  129.15656932524934,
  70.59139035957244,
  53.024744082246976,
  28.20820245736775,
  68.03861491336768,
  35.411318375397094,
  117.22114238344923,
  6.81480158563829,
  9.915921186435412,
  22.84058066284163,
  29.619554096111013,
  13.209474417145264,
  11.827571077267034,
  16.872197658289267,
  14.587326840690285,
  148.0938748630612,
  23.756885666380132,
  24.989186266702117,
  52.16638492411569,
  38.316952641698954,
  65.48612834509218,
  34.2748543385096,
  15.592457960109797,
  16.091839191889406,
  167.73628832260738,
  48.23302420932012,
  48.58294173998369,
  8.318012451740083,
  60.73228484457848,
  26.762608272901442,
  21.28631381524968,
  21.779050978050748,
  78.63091873995685,
  80.49417947643087,
  17.462517979085703,
  108.31644289537317,
  6.238157040858988,
  32.2703278753571,
  79.75699045791328,
  10.891831667939877,
  32.07951129523607,
  9.101952815995922,
  28.418149886433124,
  22.72810395713356,
  54.14937108586418,
  98.20596981584045,
  18.90434093050252,
  233.18438593839977,
  18.71682281403566,
  9.131622251901165,
  23.91949008104107,
  36.83368791552786,
  25.397582516500652,
  6.194608435214572,
  18.457356755809137,
  12.57417857976828,
  22.582416848868064,
  78.43480847732826,
  17.588622633306084,
  16.710473929955505,
  18.333106236324454,
  86.90799630936115,
  62.77681944800598,
  54.564269372991866,
  56.01944133432377,
  5.374731935814738,
  7.578622078901453,
  42.29528841899147,
  7.24634962864628,
  5.714352045277359,
  7.93963295936251,
  57.15611509588377,
  19.040792997409273,
  32.634211796866566,
  21.47843755477674,
  80.04600620140931,
  8.361298952187017,
  9.920709906400463,
  12.861796582918545,
  69.39299618892508,
  79.38038170434771,
  20.92889574020502,
  137.91747645091036,
  11.610457113291911,
  92.83697778058237,
  39.715756851215275,
  60.51001373835122,
  31.0855005571628,
  23.02139411775415,
  53.0637708230102,
  30.98964325775669,
  15.907995240216321,
  15.857485654583426,
  16.728308241688755,
  5.564303596171225,
  27.191473078258888,
  67.74682316356235,
  31.003666627545083,
  70.2291432428266,
  8.594517134429749,
  46.52990160610279,
  21.222596076104963,
  138.78863745796434,
  61.33617561125287,
  42.07395597781649,
  13.763574370901136,
  48.440778790852136,
  67.25645270511008,
  27.180690309905597,
  61.95360624180139,
  71.36649074128283,
  89.60403856537958,
  6.443113151108977,
  69.26564548931842,
  29.04465649530925,
  51.729288263465754,
  5.816443868104516,
  68.71759332290607,
  18.43620385351919,
  91.60040578478642,
  14.461280941688948,
  26.906858095035457,
  34.99723423917127,
  103.4796809426892,
  5.30784073166951,
  43.386124028612095,
  27.846708942980143,
  23.327517615929494,
  15.598648744887045,
  31.94399856468777,
  15.831728207444906,
  43.260646362737,
  10.139488701670984,
  37.258384513326874,
  42.52023423473864,
  13.062437161671866,
  9.074396730136435,
  80.3124751540522,
  69.775678350818,
  38.86698832746866,
  16.107263380105454,
  14.783698488677725,
  68.43603223327298,
  7.7399636915894,
  7.419129481939966,
  64.54525205567396,
  5.041133182690351,
  5.330204164668274,
  73.0561103338258,
  17.54441027291371,
  76.40916407082887,
  16.566690242784773,
  14.983521269682699,
  49.92082343531145,
  25.051940160910558,
  17.812736683158946,
  42.3203996633591,
  29.95656116168032,
  99.79834681844255,
  29.37950960717422,
  10.226445361361211,
  28.125321407806968,
  65.48826355593854,
  50.88249855363242,
  16.260804840900246,
  59.155422910980285,
  35.31527582744585,
  42.65180576982577,
  27.839812483869153,
  54.703754987110734,
  50.87066539301161,
  32.41372900985752,
  51.81477305323438,
  144.14270219903653,
  14.650642783225805,
  57.021905939918156,
  42.10175670835298,
  26.069872330926994,
  22.42421619960278,
  50.94666920123797,
  66.76695875207633,
  13.263261617830537,
  6.5514907396206565,
  5.359558144839126,
  77.28858932733307,
  81.43332647620434,
  26.867882672116494,
  126.00265525531108,
  8.351013104914326,
  88.19882572811005,
  101.8347967199233,
  72.72251029791693,
  36.4384575571025,
  56.591295044543656,
  40.27546416743414,
  10.823037335478448,
  18.381296649966,
  19.32987228975821,
  60.3936212901804,
  94.72190427933973,
  51.40895065123914,
  5.55431835679427,
  11.855975967339415,
  18.153972827160825,
  18.552050411019373,
  39.7182584843535,
  40.049271077718515,
  22.47539430654753,
  182.69356753659085,
  5.509846518677861,
  28.84135166895968,
  51.12545609546079,
  10.767581495949146,
  79.4841194221117,
  45.75603200191485,
  46.34940941454821,
  73.34973500824889,
  10.769812442546373,
  75.27228548640468,
  16.036874973093568,
  59.802772766331856,
  7.694964239803823,
  12.9407637201959,
  51.572087540175154,
  32.71859597754264,
  13.212258322110696,
  24.38670023020801,
  30.124652957269248,
  86.28483780167355,
  72.9400169310025,
  20.79423998308172,
  23.76429993496005,
  50.29315939663727,
  6.07658556545266,
  18.157991048637694,
  18.07517983009098,
  10.050859097506073,
  17.525070344403257,
  5.377812779467981,
  23.636555277077612,
  7.7057272640520535,
  78.68841629943648,
  20.16373155490395,
  5.146281443664074,
  61.912316977308265,
  12.00699428808731,
  24.67238943527335,
  99.29041396605,
  18.09336445328164,
  48.571636881385416,
  31.45912946661526,
  8.990295061406366,
  103.97809716307002,
  50.16633216897899,
  7.8640030826939835,
  12.149744617196212,
  7.904661609394422,
  13.238449143753693,
  91.2459983296458,
  13.265494400395282,
  42.63067409745607,
  6.566492661525614,
  88.83926493228911,
  42.48907152462587,
  81.66792630073043,
  23.282256627242273,
  7.331131452739423,
  59.09539057205833,
  19.950924659928237,
  77.79831581620144,
  7.722787896633274,
  40.52139076002147,
  25.929206594139185,
  16.69588680799624,
  81.06882832560599,
  21.893964156425387,
  5.312424091433347,
  61.822419910812826,
  17.39840254945882,
  60.46415762647308,
  13.340928216794993,
  141.18143884851554,
  17.018811193889686,
  59.02258412203743,
  116.39418952261948,
  90.94118772227516,
  70.59337289782707,
  132.56462889682385,
  14.626173391873497,
  50.29797505234686,
  28.57552352958628,
  10.519327166887544,
  37.01713326461523,
  30.75026396146229,
  14.032134907702076,
  20.425479256665515,
  50.40822569690259,
  6.328279744220919,
  13.149106538650674,
  6.450728254195457,
  5.49545863678443,
  18.656630670812575,
  40.99519087454638,
  15.200781981343658,
  14.317915559037155,
  24.868616030177666,
  10.157382009763202,
  7.326284635079966,
  13.962481840405337,
  20.012385652509902,
  5.751101387378411,
  26.20806013452035,
  34.59345917226382,
  17.5235764804364,
  20.98968102498617,
  9.056752559902506,
  175.68096518954238,
  11.359136541810196,
  15.860227939307938,
  9.105044676174387,
  15.151655117490192,
  62.32496375039588,
  24.18798496487939,
  86.60167694538129,
  80.68024762671142,
  67.97514803533488,
  7.00404445104064,
  13.492083552928797,
  24.66546823235344,
  5.1110381785127235,
  31.086559576735475,
  8.134193960647435,
  12.261396551254059,
  179.16369037630952,
  38.77790099737566,
  24.933410800817235,
  11.68303703572287,
  9.527111942095727,
  28.800696847361685,
  67.45924510874437,
  112.76853617011083,
  55.254770725096094,
  8.59952766459498,
  60.48536942173719,
  14.086083989183908,
  76.44651225158047,
  73.2761537423953,
  196.1114856797066,
  89.6537487409889,
  88.0923069077297,
  117.32547584660453,
  7.309232664133151,
  28.780179569480687,
  38.25224693036286,
  41.94458145944858,
  11.291360973015614,
  6.461562052628999,
  24.6917061140915,
  7.2242124767958185,
  10.927257236813185,
  47.62362388351568,
  183.85233359153742,
  32.579774661073756,
  64.77866611519207,
  10.838996283655824,
  41.169472353994664,
  10.369449001081755,
  56.92070214626546,
  6.982504920670008,
  15.884393659436185,
  86.06267732080241,
  5.474440189745283,
  23.47877023620519,
  94.55988509164382,
  14.910057425333306,
  7.003133155046004,
  89.6337161186774,
  75.26944018128385,
  19.089546876209987,
  58.44417850116866,
  17.9832987985655,
  54.310640335551035,
  25.530683750967224,
  125.6281123274559,
  20.36003032389581,
  9.512086670118443,
  17.386667160663894,
  135.85215721557338,
  14.791883535850683,
  71.79947226306216,
  25.192155261843396,
  20.283515577844902,
  70.61856220748568,
  6.6722021539051735,
  37.21439795166085,
  39.96894729959123,
  8.990523694002215,
  6.1627703443934285,
  56.81213546390012,
  23.358726498640763,
  7.556569321964341,
  20.265975355084027,
  8.792474473951097,
  26.385622965167485,
  92.85036261086564,
  43.452585673475284,
  12.816397223053084,
  19.23800909366875,
  31.621045576972197,
  63.525278237650426,
  43.28584199422976,
  12.300678089842004,
  42.72075730552848,
  116.84460293206509,
  25.53864554984964,
  9.766984687656507,
  49.37138930022858,
  45.176521371135365,
  65.8531420050303,
  6.245036241017265,
  15.925733943875045,
  8.609329995630876,
  20.695347864925772,
  28.32518865582491,
  6.447290640514619,
  12.939082660556197,
  26.712849786981838,
  31.215084883327602,
  7.95110559306821,
  20.641922700487836,
  28.365839962095663,
  77.63607194104131,
  5.976337210244874,
  38.35073950444072,
  26.382361538335587,
  10.206453831429773,
  31.376326235489408,
  11.3418073405289,
  55.8422275790988,
  133.5400734161776,
  33.09038242349268,
  180.12528554154386,
  14.992298895427966,
  75.67504340030679,
  10.589535888279855,
  111.67574095779477,
  15.22930633244243,
  119.76377374978344,
  60.58671191836012,
  52.538773934390456,
  20.7863115040151,
  42.56229480817096,
  15.374737044093054,
  95.77710169576866,
  127.69122158678363,
  95.94470632849146,
  20.795425774432267,
  16.10995609864871,
  5.941573585537059,
  19.21111454187715,
  39.609573071167645,
  34.07889488678319,
  74.77133389595116,
  7.960237648452039,
  42.22765535669677,
  86.78925805911857,
  5.295578478558962,
  45.4946385886513,
  111.65075026848433,
  15.807871731274638,
  71.10647599244821,
  25.796377837171956,
  5.491799561378873,
  27.24999653394341,
  9.264108361484686,
  10.25461058113253,
  8.781940846114253,
  24.403069070465996,
  25.56000560421904,
  12.419381338979749,
  107.3726571542254,
  70.06755545484957,
  22.49616117626991,
  103.59667879436307,
  5.329263754375793,
  49.81147627110015,
  38.684270314870076,
  10.499787209669593,
  13.725060745174016,
  57.077735427687706,
  30.675260319551306,
  68.57777226656908,
  29.339771056622904,
  19.339551561044335,
  9.010408070409003,
  10.192559568245219,
  86.68882818776349,
  37.52342577134934,
  17.937975413383175,
  68.64176254882415,
  58.402195942291335,
  61.007124080438274,
  20.08474691225381,
  38.12336010072819,
  8.16318340389225,
  10.783065820229,
  59.217371572387265,
  34.402161605715136,
  15.550415113744474,
  78.06052672145165,
  14.547449226693427,
  58.9723803443803,
  42.5418929792828,
  45.59566316994359,
  6.4204890279454245,
  17.16039493547367,
  20.78250383230657,
  46.72032240025519,
  17.957708277714797,
  72.29367086727052,
  62.54387065716355,
  12.707495638153997,
  14.630656905864686,
  24.269312420828086,
  67.79129075254171,
  91.66759331944496,
  131.746227286324,
  50.32337164519672,
  51.37164059874975,
  8.017361132634237,
  75.58998345415668,
  52.676468565988024,
  22.84627108128996,
  13.663940877592156,
  15.17630386222057,
  101.21006935474807,
  6.724315397593212,
  13.3341723286413,
  14.69499562292919,
  60.461644636278635,
  45.253975109568394,
  47.63846298454778,
  53.35919978459143,
  13.479818128479813,
  10.55344042729989,
  77.45686375400152,
  16.375809228978582,
  9.408466211584761,
  66.4613206584992,
  64.94149757647627,
  70.6780238001515,
  82.74955534392683,
  73.85440796874556,
  58.33454333398977,
  44.20099986827395,
  10.09964640222436,
  17.344598528906207])

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)