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

([4377943.284777576,
  12255914.0625,
  12256329.440357525,
  12350185.856029617,
  12350435.9375,
  12350455.795000589,
  12354271.41548004,
  12355335.9375,
  12355503.667421496,
  12357477.079338957,
  12357818.614609193,
  12358233.043507664,
  12358295.461459577,
  12358806.25,
  12358926.464205928,
  12359014.0625,
  12359614.0625,
  12359625.573930643,
  12359647.259285161,
  12361462.5,
  12362081.25,
  12363108.426152745,
  12364302.568847135,
  12364984.375,
  12365433.53820789,
  12366037.5,
  12366396.668463694,
  12367611.886960475,
  12367621.875,
  12367638.369155396,
  12367762.343800675,
  12368516.549478827,
  12368541.515387911,
  12370076.5625,
  12370356.204935547,
  12370868.607021285,
  12373801.034123015,
  12373909.375,
  12374256.25,
  12374508.079869175,
  12374638.46673763,
  12374809.375,
  12374815.279637367,
  12375015.215434952,
  12375051.384495022,
  12375280.00012539,
  12375770.359766342,
  12375889.62138288,
  12375995.3125,
  12376019.614683805,
  12376542.779199434,
  12376676.212086761,
  12376758.560438618,
  12377829.6875,
  12377854.6875,
  12378148.085868357,
  12378183.181248715,
  12378411.063585557,
  12378417.581000097,
  12378618.638100153,
  12378670.148967575,
  12378772.644742787,
  12378923.4375,
  12378938.26618486,
  12378976.44432619,
  12379085.246972604,
  12379132.123131301,
  12379310.842022989,
  12379349.519611906,
  12379478.271606395,
  12379509.959710054,
  12379799.350521581,
  12380064.317919713,
  12380101.5625,
  12380143.75,
  12380200.0,
  12380218.75,
  12380274.852569217,
  12380294.331618657,
  12381185.9375,
  12381279.6875,
  12381288.470141595,
  12381313.644901635,
  12381318.75,
  12381426.722140273,
  12381432.021509686,
  12381864.343178418,
  12382220.3125,
  12382343.75,
  12382349.039412165,
  12382441.851744793,
  12382465.867682574,
  12382475.164130617,
  12382577.896651376,
  12382610.805211566,
  12382774.694935815,
  12382960.483284507,
  12383078.439428175,
  12383258.472484248,
  12383308.977514647,
  12383321.652865456,
  12383435.582282286,
  12383533.493772384,
  12383562.704805171,
  12383843.05684846,
  12384185.9375,
  12384298.655672044,
  12384396.331523027,
  12385641.562205112,
  12385953.125,
  12386310.92283924,
  12386329.147376813,
  12386360.819116259,
  12387310.825105557,
  12387451.47191883,
  12387648.4375,
  12388776.5625,
  12389169.561119694,
  12389481.932590725,
  12389542.1875,
  12389628.740562841,
  12389871.064415647,
  12390442.721593175,
  12390457.336101025,
  12390826.210399674,
  12393743.75,
  12393751.355175387,
  12393857.8125,
  12394288.840428794,
  12394389.550966877,
  12394433.437107764,
  12394705.25648869,
  12394729.6875,
  12394862.5,
  12395015.0445802,
  12395199.845917802,
  12395277.452527585,
  12395645.3125,
  12395672.214333445,
  12395982.8125,
  12396072.46789041,
  12396078.125,
  12396109.798390329,
  12396301.380054422,
  12396313.64245891,
  12396635.9375,
  12396658.224157428,
  12396751.944499362,
  12396796.508937938,
  12396881.25,
  12397480.787809664,
  12397582.7021443,
  12397798.4375,
  12397985.9375,
  12398543.49004211,
  12398870.818241073,
  12399457.8125,
  12399552.3246644,
  12399707.8125,
  12401185.927070305,
  12401365.61845149,
  12404167.073746992,
  12404339.319597904,
  12404643.562496297,
  12405013.269530231,
  12405142.1875,
  12406502.663532255,
  12407585.9375,
  12408482.8125,
  12410842.1875,
  12414386.349441666,
  12416169.084054116,
  12417307.776898224,
  12418067.144471522,
  12444225.697179591,
  12444520.810785094,
  12444631.357740153,
  12445208.784452088,
  12445995.3125,
  12449925.584971754,
  12451404.6875,
  12451457.8125,
  12452351.66657799,
  12455740.625,
  12455882.672177993,
  12456330.724715944,
  12457113.081415823,
  12457602.50124031,
  12458931.56982687,
  12459931.140894866,
  12460046.734162364,
  12460122.542642342,
  12460406.920009531,
  12460532.8125,
  12460666.284695176,
  12460679.6875,
  12460686.068091732,
  12460816.307987213,
  12460976.5625,
  12460999.144124873,
  12461115.465272877,
  12461380.304726053,
  12461489.932332514,
  12461561.080045262,
  12461745.70062896,
  12461885.600787194,
  12461991.573433017,
  12462043.193655215,
  12462136.029631585,
  12462443.793387676,
  12462461.386566192,
  12462659.24080902,
  12463060.080401521,
  12463799.766469764,
  12464104.278917091,
  12464865.50202505,
  12465157.656546237,
  12465292.1875,
  12465684.375,
  12465959.748769395,
  12466034.807243405,
  12466044.15024325,
  12466137.5,
  12466221.875,
  12467146.756972058,
  12468006.498682642,
  12470824.582563516,
  12471129.6875,
  12473804.6875,
  12474046.875,
  12474610.9375,
  12474653.125,
  12477406.079069437,
  12477711.960540265,
  12477880.690797407,
  12477902.313297536,
  12478436.948230572,
  12482798.781729758,
  12482814.766896904,
  12482963.427293448,
  12482982.8125,
  12483013.539455382,
  12483075.977168985,
  12483080.980597656,
  12483116.010108437,
  12483134.941611752,
  12483142.1875,
  12483190.03370018,
  12483199.903074529,
  12483263.069600042,
  12483271.292567734,
  12483373.076325474,
  12483407.612917833,
  12483433.410456043,
  12483437.914202208,
  12483440.997874942,
  12483441.835304867,
  12483470.057304375,
  12483475.0,
  12483531.490946762,
  12483534.375,
  12483553.803175142,
  12483613.17324111,
  12483614.0625,
  12483614.252150819,
  12483625.945432298,
  12483626.5625,
  12483699.233020019,
  12483704.6875,
  12483709.375,
  12483716.299569217,
  12483717.1875,
  12483739.0625,
  12483750.0,
  12483802.927420411,
  12483808.967707543,
  12483839.0625,
  12483839.0625,
  12483875.897974633,
  12483889.693701714,
  12483922.533697639,
  12483933.90565027,
  12483949.15722256,
  12484039.670209765,
  12484048.4375,
  12484051.195774453,
  12484084.986598872,
  12484088.907235784,
  12484106.25,
  12484116.640465459,
  12484125.150760682,
  12484142.486732727,
  12484172.874919172,
  12484175.0,
  12484249.44943961,
  12484252.673398979,
  12484276.215081908,
  12484334.375,
  12484351.505218223,
  12484368.296546234,
  12484377.051737294,
  12484396.348905586,
  12484452.359739888,
  12484513.632896487,
  12484523.453744859,
  12484529.368383035,
  12484531.850962084,
  12484535.746577486,
  12484537.442165133,
  12484541.447840609,
  12484543.631709643,
  12484554.372816136,
  12484576.517665122,
  12484578.125,
  12484578.281292558,
  12484585.9375,
  12484589.0625,
  12484602.589748371,
  12484606.198046824,
  12484624.281988807,
  12484650.787377324,
  12484660.9375,
  12484675.752636196,
  12484676.521624032,
  12484678.728347644,
  12484679.423367925,
  12484681.671221549,
  12484690.97312798,
  12484704.6875,
  12484707.203258539,
  12484716.361892644,
  12484718.799264,
  12484733.400892474,
  12484762.5,
  12484763.852532241,
  12484791.783537438,
  12484803.125,
  12484833.165727066,
  12484841.065890003,
  12484842.668689597,
  12484843.75,
  12484864.0625,
  12484864.721153002,
  12484865.625,
  12484865.86928223,
  12484914.0625,
  12484920.3125,
  12484921.31379078,
  12484921.875,
  12484932.8125,
  12484934.25045108,
  12484940.870542074,
  12484946.751185931,
  12484947.693685168,
  12484951.317583993,
  12484962.606082259,
  12484987.477190815,
  12484988.530671118,
  12485000.654375043,
  12485010.133116776,
  12485012.5,
  12485021.773685321,
  12485033.204493854,
  12485040.906354176,
  12485044.034115646,
  12485045.259240624,
  12485072.388641825,
  12485109.375,
  12485132.670030052,
  12485138.156129232,
  12485149.949379055,
  12485153.935033027,
  12485159.3514224,
  12485160.825105151,
  12485162.5,
  12485174.850765664,
  12485181.25,
  12485182.573843576,
  12485187.5,
  12485189.0625,
  12485202.495316336,
  12485208.658545874,
  12485233.728667852,
  12485272.920024453,
  12485277.311232613,
  12485280.789746601,
  12485331.432438534,
  12485334.375,
  12485343.78695632,
  12485401.400514033,
  12485406.064010208,
  12485419.203550613,
  12485443.767725334,
  12485451.35982336,
  12485457.8125,
  12485470.3125,
  12485485.9375,
  12485511.646069713,
  12485521.155412925,
  12485554.929352881,
  12485570.178409105,
  12485573.4375,
  12485587.5,
  12485593.75,
  12485608.543041794,
  12485674.996756652,
  12485682.860724129,
  12485703.635479359,
  12485718.353978677,
  12485740.205124361,
  12485744.129773118,
  12485801.5625,
  12485804.97459279,
  12485810.893955808,
  12485811.014899192,
  12485824.77754257,
  12485830.6829768,
  12485870.914166197,
  12485899.37263224,
  12485906.24892023,
  12485907.907916911,
  12485914.409838758,
  12485917.1875,
  12485953.894877963,
  12485965.741874095,
  12486070.3125,
  12486078.97830789,
  12486079.898505686,
  12486112.332051605,
  12486131.25,
  12486190.159488682,
  12486190.366810668,
  12486201.130727619,
  12486254.6875,
  12486268.318494622,
  12486546.875,
  12486573.4375,
  12486633.783059942,
  12486809.962634478,
  12486896.875,
  12487112.271738444,
  12487461.508081421,
  12487511.743821198,
  12487782.071446456,
  12487783.011805061,
  12487837.744988529,
  12487839.0625,
  12487930.1051326,
  12488016.141207064,
  12488045.239556307,
  12488153.125,
  12488211.537450785,
  12488238.99854379,
  12488263.192916077,
  12488299.45941862,
  12488322.676409803,
  12488359.375,
  12488361.074303394,
  12488380.635380907,
  12488390.545051357,
  12488414.933630584,
  12488435.440400722,
  12488442.1875,
  12488495.273064727,
  12488496.875,
  12488510.59185944,
  12488544.17255787,
  12488567.524606531,
  12488588.327682339,
  12488608.64707329,
  12488635.9375,
  12488638.870789122,
  12488667.68441781,
  12488681.359977726,
  12488685.9375,
  12488721.771628581,
  12488746.071470276,
  12488753.7384552,
  12488770.759002756,
  12488816.571193267,
  12488860.083632132,
  12488957.8125,
  12488963.59457016,
  12488999.754093206,
  12489067.860075442,
  12489071.898554211,
  12489102.014430651,
  12489186.43841051,
  12489190.484312676,
  12489196.35716567,
  12489212.308330331,
  12489238.785092873,
  12489296.875,
  12489395.3125,
  12489441.442679115,
  12489512.796229681,
  12489514.0625,
  12489563.984592019,
  12489575.0,
  12489585.391687103,
  12489596.518509097,
  12489598.261168296,
  12489618.75,
  12489621.008741235,
  12489631.018182522,
  12489644.809926324,
  12489656.47743189,
  12489677.966998324,
  12489706.25,
  12489711.284182452,
  12489749.100217065,
  12489771.798574172,
  12489792.843319522,
  12489812.5,
  12489817.80145353,
  12489839.207011605,
  12489908.81881227,
  12489927.645674584,
  12489927.714997852,
  12489932.8125,
  12489960.9375,
  12490130.883691173,
  12490156.761285847,
  12490193.10247562,
  12490204.907194976,
  12490257.699484117,
  12490289.787799528,
  12490303.941453287,
  12490343.75,
  12490365.3529407,
  12490407.8125,
  12490438.5473968,
  12490521.02339622,
  12490593.814094542,
  12490627.453184415,
  12490670.132433543,
  12490683.75735773,
  12490686.518592007,
  12490737.119451202,
  12490820.022300342,
  12490876.217348158,
  12490919.734678585,
  12490931.25,
  12491022.662480893,
  12491156.25,
  12491162.009963434,
  12491190.625,
  12491257.8125,
  12491296.875,
  12491399.153855111,
  12491449.426244935,
  12491600.223154653,
  12491624.939001372,
  12491734.521849304,
  12491973.87556011,
  12492004.529601667,
  12492169.688616702,
  12492319.02989759,
  12492461.89282896,
  12492468.564718042,
  12492498.198679073,
  12492513.327596009,
  12492564.153597169,
  12492623.288618093,
  12492701.5625,
  12493039.011530682,
  12493847.03076862,
  12493918.350042826,
  12493926.216887426,
  12493959.204417396,
  12494196.44261006,
  12494491.889644405,
  12494646.765198428,
  12494701.420367543,
  12495587.958964828,
  12495743.75,
  12495754.568365885,
  12495832.56004037,
  12495894.115550946,
  12496513.52320837,
  12496560.558882074,
  12496692.476343887,
  12496884.375,
  12496895.3125,
  12496900.0,
  12497196.875,
  12497247.085662022,
  12497724.913023826,
  12497737.5,
  12498024.80860055,
  12498070.606707606,
  12498298.361328416,
  12498300.0,
  12498333.840802947,
  12498648.014628079,
  12498769.815556247,
  12498871.562471412,
  12499148.4375,
  12500088.826117525,
  12501196.50931726,
  12501992.18207627,
  12502064.374650987,
  12503242.597424798,
  12504282.8125,
  12506257.252408797,
  12506575.770622978,
  12506587.5,
  12507656.25,
  12508159.535052627,
  12508375.83123053,
  12509292.051979573,
  12509547.41517504,
  12511390.206322545,
  12516028.126519538,
  373853869.6446979,
  373853870.12083316,
  373853870.8804961,
  373853873.758095,
  373853901.7736579],
 [7.004302851715175,
  54.2576156490845,
  52.7827389588267,
  50.357055805034804,
  84.88691555930873,
  13.451005207967631,
  37.474570952328,
  48.31037786080955,
  70.93875680701878,
  17.36086264441774,
  92.87774448067515,
  7.688982170720358,
  6.793614087306694,
  44.60195710740117,
  7.5680691137891,
  51.26032352770066,
  52.775889813085335,
  21.43108277118567,
  16.640939237347535,
  35.035939337396606,
  36.6024684837862,
  18.209113547467805,
  43.31381393721329,
  88.67862736395894,
  6.363457327753961,
  66.0216006285424,
  8.180296070374004,
  13.026764143938262,
  34.156679043828504,
  19.75929476696141,
  25.77954790495339,
  32.780327284854856,
  12.878156688888136,
  44.877564915847245,
  7.587019555074294,
  10.911341887162749,
  72.88546053844883,
  38.95503365303942,
  48.2157197870433,
  16.068933879874905,
  20.58820386433075,
  54.09716674459274,
  10.627632353280555,
  150.66802807242107,
  15.352020608677297,
  16.283923237368274,
  8.311731724044057,
  12.13023167617171,
  72.92585432215718,
  19.723332358476174,
  9.887409569480281,
  54.62312045103624,
  18.840361418313126,
  53.059243698315264,
  28.984531191684297,
  92.98324569478643,
  87.54781679144352,
  13.311946311933655,
  196.2597482540032,
  143.58960223654444,
  74.24594720791205,
  18.997492205317556,
  39.05979334971462,
  12.581703743166399,
  10.947884422484599,
  47.44537166522119,
  39.65412861729904,
  26.769453105470635,
  58.741996208758316,
  58.039444158532234,
  34.21384982641598,
  5.884772675791998,
  16.01345310123795,
  56.18019973872723,
  41.80077085039523,
  52.7835127493145,
  42.874367104633464,
  75.39338216402459,
  57.381499354199875,
  82.75570976709292,
  30.814960507178643,
  20.618730001447922,
  13.557536093316953,
  84.07985093763213,
  7.1640582108129,
  71.49154059157988,
  7.1972412858931705,
  74.37991624126794,
  50.9426347706555,
  26.96445422676303,
  13.799568108258393,
  29.654766555404105,
  11.394325189897078,
  99.40537091712687,
  72.169510322715,
  13.893637795889227,
  19.981690038911612,
  16.292650661737497,
  19.178135176252784,
  55.903287662091294,
  40.087610994720364,
  6.051922579601536,
  18.380163276680474,
  29.833401939009743,
  64.75096314754072,
  41.69102500136764,
  70.20106682384063,
  10.372661754478065,
  20.328089075063424,
  39.78241579450581,
  7.117603114454929,
  8.912957808242266,
  7.56097423474301,
  6.741273434609248,
  157.5087950932332,
  38.502386929481354,
  47.90003321755987,
  12.686550200455416,
  9.98397618907522,
  73.5716276338428,
  7.293245638372215,
  18.868770062679392,
  17.40085327646871,
  13.391151812794407,
  76.95218041428222,
  143.84704343640388,
  21.019144970879033,
  119.70922790299124,
  24.79997120257438,
  11.987920190635021,
  47.51887547267849,
  54.846988076010426,
  35.26507073396717,
  64.42799532131818,
  87.17265577492589,
  45.72133186257568,
  50.342673504470895,
  54.75249330360912,
  16.414537058338713,
  33.913248319500404,
  12.200874101335492,
  31.701004971059195,
  5.106881428560929,
  25.378212555741918,
  145.50320410519132,
  31.20017081253198,
  16.783007181597725,
  7.22562359797847,
  22.942040776782378,
  91.59286340540658,
  5.381094068944606,
  14.53693168268741,
  61.05627494577508,
  56.30043106973886,
  23.92231310381323,
  5.177309726449087,
  52.91093571436809,
  122.26197683506172,
  39.60483104421091,
  21.435926330469293,
  43.4313252792459,
  17.676044869099616,
  43.51183254410058,
  6.243795761382695,
  31.05185256128948,
  50.585746674960376,
  23.92848960483153,
  62.18069239650571,
  82.51481921657395,
  77.05730428469288,
  94.48178940736011,
  48.11122034051288,
  34.74683324295322,
  47.18362293916658,
  18.487298006451784,
  92.5403351761172,
  7.642675884445015,
  16.543419812152752,
  44.65363405110824,
  78.68769895852688,
  55.11437889409917,
  33.731281491752995,
  9.420985272507867,
  39.54605405314805,
  41.56544720619895,
  13.257328290267752,
  9.876691390291787,
  20.34140535493127,
  7.723159920603721,
  14.541343412468777,
  17.783245530486912,
  18.289647394742786,
  9.970409553460744,
  55.846599662397956,
  9.959335835894477,
  50.64429440312175,
  26.403157272545364,
  12.394371238955232,
  78.96261134634669,
  130.59362756550067,
  13.148077411442976,
  70.05581898248795,
  10.193893022464694,
  96.32368503245443,
  7.749920382279488,
  36.99474452122499,
  13.024135217484517,
  49.613292408978396,
  5.132476107843655,
  7.619006009319628,
  5.9794908351563265,
  14.453569057794926,
  167.66465604011847,
  36.459816197461514,
  21.668572589839002,
  29.199583527466935,
  26.96100048679822,
  44.30584305497363,
  43.36511135554704,
  11.766920959404809,
  16.89203317410678,
  16.723368278439867,
  38.29377460664038,
  79.41469961341296,
  15.880273033591848,
  30.467718351745688,
  22.358247053487638,
  31.71109269081241,
  53.34441290664619,
  59.06453587539118,
  59.63862412631337,
  34.535156214472444,
  14.327709857483987,
  47.17062367592397,
  101.99609410893586,
  18.860771961106515,
  20.99176312790026,
  27.327591941047455,
  18.536447383394865,
  92.14102266984797,
  71.342591476699,
  13.263372264995027,
  50.776525750768364,
  86.13653436808794,
  16.645164892386383,
  64.83821524928753,
  72.52555679904937,
  25.235613953484524,
  14.584511118508702,
  8.612655771167354,
  94.82866937456829,
  75.72663622248855,
  33.43178986384958,
  29.536987563039947,
  29.206621540734183,
  24.177563463963,
  17.612723189326772,
  24.340439120916738,
  94.22795613391362,
  9.543800394052464,
  58.82793264271234,
  94.47109449229434,
  5.220920218424627,
  82.83838358595963,
  8.389480426354208,
  30.34547082245505,
  69.86694020126463,
  53.42303347536786,
  52.550420806932244,
  46.18233416077441,
  6.906090549187321,
  51.157083810099515,
  42.80377581257708,
  34.40652033329609,
  9.169745327021296,
  13.587329868822412,
  55.087905846047995,
  52.88516976655463,
  123.03580849722142,
  9.928969637227564,
  5.750002162543088,
  28.25382798963022,
  26.33959711933607,
  12.232059925358381,
  72.18310174586425,
  6.050775552412595,
  9.908268129065732,
  108.26000720346725,
  54.28341322582772,
  21.03639537707841,
  19.799223952626363,
  77.44624123492947,
  37.59402633192926,
  73.42979569055228,
  69.12461772150091,
  13.460583203859652,
  63.24880513424552,
  45.98950820977657,
  65.55827183899794,
  21.970339581547233,
  11.989449085009378,
  21.236684927961594,
  174.409413670405,
  107.92047303779387,
  141.05378729203525,
  74.63349848735584,
  17.863132382222727,
  25.251780797547397,
  11.03849469217633,
  10.139757204859007,
  26.38701041786543,
  10.669054229905635,
  27.667725661070612,
  55.23129586761987,
  12.556336105944542,
  68.2992782217281,
  36.66744423772254,
  41.17865015946295,
  27.537142937699215,
  119.90196614019536,
  19.12708698036304,
  65.63377074346508,
  5.6570844134998675,
  57.30944142408754,
  17.8779411943219,
  125.64664343840886,
  28.409435250080143,
  28.764762739705642,
  93.55518154446904,
  7.988695433015603,
  153.09994430758638,
  73.61869899666469,
  9.883958342222924,
  89.81762025170687,
  24.97319244343395,
  10.549999342514699,
  67.27213341003649,
  26.6067678600676,
  27.084375626173784,
  9.774357256048104,
  122.67576648572899,
  52.04371866078608,
  11.689104752251565,
  41.16667644221455,
  78.74668938720963,
  100.06660526156661,
  73.68786333515878,
  28.854773850467076,
  34.19099954820316,
  29.28381452393699,
  47.88492508774891,
  14.29782781489353,
  38.619368181779244,
  12.704465547437419,
  78.3327196547419,
  65.67829758080144,
  116.62608577174841,
  8.918852869246702,
  9.953641023560472,
  55.56046798900046,
  67.95776145625702,
  26.67033309422833,
  14.735778190697665,
  44.44708721087042,
  11.60286912011437,
  59.979078412341295,
  5.740246385936458,
  72.16279462989951,
  25.9956460623521,
  8.666961055875104,
  27.56147093067576,
  6.622038386396409,
  28.066581979028847,
  85.8965392046531,
  32.507970542814704,
  53.75308184294728,
  61.99086464273856,
  43.13756800121971,
  33.42272652884848,
  42.491834194569535,
  12.984372213600585,
  109.16699091385988,
  7.972004532452814,
  115.95675369980508,
  31.60747703140596,
  8.970362227862847,
  5.348212292398322,
  170.32447438112524,
  25.39461971029664,
  25.689753420439597,
  128.4108574894461,
  22.923517291466485,
  23.879023214921144,
  32.850704516803845,
  94.241834519035,
  47.674956807162374,
  42.874078537524895,
  267.69045911495505,
  26.979484728330867,
  7.693118407566914,
  26.129684498104048,
  79.94621217989265,
  72.29553750566461,
  50.67407470133528,
  11.158528124498075,
  82.79205761535891,
  19.157737790208685,
  17.32440397627579,
  20.555302571995973,
  10.529247429324672,
  16.617831539622564,
  79.30467464529444,
  9.586556529180305,
  27.691560468355817,
  15.244960173626433,
  14.156378814162549,
  28.312833940843245,
  144.42400385897858,
  30.843951686429108,
  28.503842966598675,
  33.73055398605692,
  11.717792275688392,
  59.730785744413424,
  98.7516147375063,
  27.017210944780526,
  72.69447911598496,
  29.387392340651942,
  55.48065694190197,
  25.598014438464617,
  56.37725560670995,
  10.470496341939722,
  10.746776600343148,
  10.513823457990869,
  48.16220169554187,
  88.87974662269184,
  43.000628041814466,
  45.45810061817287,
  10.313405625864716,
  5.8587784221719765,
  56.65108657998005,
  63.16020871410067,
  5.744059779469008,
  19.26888098624897,
  12.711023549112907,
  33.50586661377136,
  5.171106749425389,
  83.62372887334072,
  5.076570316611948,
  9.80967873342721,
  32.90073091688129,
  71.11619460163762,
  7.83566819401059,
  135.74274524315902,
  97.51342640508696,
  13.21626051679822,
  19.231859023579307,
  86.74391518532043,
  18.378008196723776,
  70.00226757629704,
  27.03744922006843,
  65.00824796111104,
  81.69910871291269,
  150.24029501106264,
  13.5078216913279,
  87.52256073153305,
  41.99929853602795,
  16.840500379771875,
  16.40353747251429,
  12.725871531799617,
  24.556200229844272,
  38.991822233074046,
  9.1743961229883,
  12.004724580251052,
  11.30407325217084,
  28.82335480698434,
  76.9413871671114,
  12.565114759807928,
  12.244217241795173,
  16.968679101994923,
  10.284049698143138,
  5.568490345122137,
  38.709524829487066,
  13.414745278483165,
  24.467206737540007,
  18.32227261747431,
  9.351399319684976,
  16.998101233912426,
  7.259660085879079,
  94.82048750778833,
  26.50402712375147,
  60.45146088997449,
  24.04368896093938,
  34.22015423624592,
  49.02417389527286,
  12.70207873920776,
  6.824123533046316,
  60.22696514383394,
  59.03014922333422,
  37.33383761241198,
  33.529222118726295,
  23.055007665804016,
  25.470352667080007,
  40.60119960686165,
  18.480932845630722,
  24.662584614048853,
  21.467464795378213,
  24.434821051105246,
  25.75149895464371,
  28.708874142998877,
  11.71658061440018,
  8.67115210440233,
  5.701911266958268,
  33.470786470787885,
  103.39594277767675,
  129.0272434588124,
  6.103257190964134,
  53.80152999506884,
  15.18107246135871,
  13.57947767131657,
  83.12847209064773,
  37.91035844579114,
  5.035424514497084,
  90.07608276274134,
  20.137724540153233,
  15.850670565731386,
  104.74351913363358,
  54.972242220134206,
  19.345840830152646,
  52.315840961158194,
  21.827068816322328,
  69.80099229568637,
  13.284871379796368,
  12.456256639028917,
  11.229390706034335,
  12.878889905664694,
  14.296897577943332,
  77.82817288487384,
  6.582211199591967,
  13.665549855630045,
  77.7139704877213,
  9.055668241211315,
  5.772428704689789,
  71.80452160687348,
  17.060298215282053,
  101.57377530970028,
  5.676264744774738,
  62.753506469608396,
  74.6317907709306,
  47.66046472211516,
  12.468472489162199,
  5.772712952145254,
  15.866236108754977,
  5.289206764185928,
  9.45873791405088,
  63.797398144249144,
  110.19268323613245,
  89.99823758314976,
  16.12767166966444,
  20.486722472118874,
  25.33573871581273,
  82.96305957983394,
  34.740123416116816,
  8.335345774516425,
  7.550262830344939,
  63.97387455254626,
  11.048962226756503,
  5.622720151163016,
  22.552717093928834,
  38.902615652346185,
  5.80146551371251,
  13.515108358102598,
  21.545861529608064,
  21.17567992026929,
  19.985889717332828,
  33.8836065728687,
  67.14101480271866,
  16.26312299478976,
  24.377611928969046,
  5.173963365721217,
  13.219708422485201,
  22.795897220886978,
  9.588192821459481,
  60.719493319813644,
  31.887376943856992,
  61.63608914382314,
  72.03500696873377,
  45.43001952559559,
  9.258926999080897,
  38.95530840917322,
  10.841160035606705,
  9.593211240830037,
  16.7112933995029,
  43.773150754746744,
  21.158617058540894,
  69.0223895100502,
  63.2394972887217,
  5.596985157989508,
  87.56914860691914,
  58.911324525149524,
  13.70960487998605,
  14.920876919768775,
  6.121173059664683,
  6.131791039219603,
  91.84862022944479,
  5.43822957912523,
  8.740972912735097,
  68.20349774504608,
  45.2005446764171,
  15.588788948110123,
  90.53629904683305,
  5.927689844767063,
  17.496416889395206,
  34.3152364079604,
  7.1218511547982555,
  354.39322684690893,
  167.89540900358708,
  17.66175664818565,
  18.934687567299672,
  158.36026077864338])

Energy spectrum

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

([4377943.284777576,
  12255914.0625,
  12256329.440357525,
  12350185.856029617,
  12350435.9375,
  12350455.795000589,
  12354271.41548004,
  12355335.9375,
  12355503.667421496,
  12357477.079338957,
  12357818.614609193,
  12358233.043507664,
  12358295.461459577,
  12358806.25,
  12358926.464205928,
  12359014.0625,
  12359614.0625,
  12359625.573930643,
  12359647.259285161,
  12361462.5,
  12362081.25,
  12363108.426152745,
  12364302.568847135,
  12364984.375,
  12365433.53820789,
  12366037.5,
  12366396.668463694,
  12367611.886960475,
  12367621.875,
  12367638.369155396,
  12367762.343800675,
  12368516.549478827,
  12368541.515387911,
  12370076.5625,
  12370356.204935547,
  12370868.607021285,
  12373801.034123015,
  12373909.375,
  12374256.25,
  12374508.079869175,
  12374638.46673763,
  12374809.375,
  12374815.279637367,
  12375015.215434952,
  12375051.384495022,
  12375280.00012539,
  12375770.359766342,
  12375889.62138288,
  12375995.3125,
  12376019.614683805,
  12376542.779199434,
  12376676.212086761,
  12376758.560438618,
  12377829.6875,
  12377854.6875,
  12378148.085868357,
  12378183.181248715,
  12378411.063585557,
  12378417.581000097,
  12378618.638100153,
  12378670.148967575,
  12378772.644742787,
  12378923.4375,
  12378938.26618486,
  12378976.44432619,
  12379085.246972604,
  12379132.123131301,
  12379310.842022989,
  12379349.519611906,
  12379478.271606395,
  12379509.959710054,
  12379799.350521581,
  12380064.317919713,
  12380101.5625,
  12380143.75,
  12380200.0,
  12380218.75,
  12380274.852569217,
  12380294.331618657,
  12381185.9375,
  12381279.6875,
  12381288.470141595,
  12381313.644901635,
  12381318.75,
  12381426.722140273,
  12381432.021509686,
  12381864.343178418,
  12382220.3125,
  12382343.75,
  12382349.039412165,
  12382441.851744793,
  12382465.867682574,
  12382475.164130617,
  12382577.896651376,
  12382610.805211566,
  12382774.694935815,
  12382960.483284507,
  12383078.439428175,
  12383258.472484248,
  12383308.977514647,
  12383321.652865456,
  12383435.582282286,
  12383533.493772384,
  12383562.704805171,
  12383843.05684846,
  12384185.9375,
  12384298.655672044,
  12384396.331523027,
  12385641.562205112,
  12385953.125,
  12386310.92283924,
  12386329.147376813,
  12386360.819116259,
  12387310.825105557,
  12387451.47191883,
  12387648.4375,
  12388776.5625,
  12389169.561119694,
  12389481.932590725,
  12389542.1875,
  12389628.740562841,
  12389871.064415647,
  12390442.721593175,
  12390457.336101025,
  12390826.210399674,
  12393743.75,
  12393751.355175387,
  12393857.8125,
  12394288.840428794,
  12394389.550966877,
  12394433.437107764,
  12394705.25648869,
  12394729.6875,
  12394862.5,
  12395015.0445802,
  12395199.845917802,
  12395277.452527585,
  12395645.3125,
  12395672.214333445,
  12395982.8125,
  12396072.46789041,
  12396078.125,
  12396109.798390329,
  12396301.380054422,
  12396313.64245891,
  12396635.9375,
  12396658.224157428,
  12396751.944499362,
  12396796.508937938,
  12396881.25,
  12397480.787809664,
  12397582.7021443,
  12397798.4375,
  12397985.9375,
  12398543.49004211,
  12398870.818241073,
  12399457.8125,
  12399552.3246644,
  12399707.8125,
  12401185.927070305,
  12401365.61845149,
  12404167.073746992,
  12404339.319597904,
  12404643.562496297,
  12405013.269530231,
  12405142.1875,
  12406502.663532255,
  12407585.9375,
  12408482.8125,
  12410842.1875,
  12414386.349441666,
  12416169.084054116,
  12417307.776898224,
  12418067.144471522,
  12444225.697179591,
  12444520.810785094,
  12444631.357740153,
  12445208.784452088,
  12445995.3125,
  12449925.584971754,
  12451404.6875,
  12451457.8125,
  12452351.66657799,
  12455740.625,
  12455882.672177993,
  12456330.724715944,
  12457113.081415823,
  12457602.50124031,
  12458931.56982687,
  12459931.140894866,
  12460046.734162364,
  12460122.542642342,
  12460406.920009531,
  12460532.8125,
  12460666.284695176,
  12460679.6875,
  12460686.068091732,
  12460816.307987213,
  12460976.5625,
  12460999.144124873,
  12461115.465272877,
  12461380.304726053,
  12461489.932332514,
  12461561.080045262,
  12461745.70062896,
  12461885.600787194,
  12461991.573433017,
  12462043.193655215,
  12462136.029631585,
  12462443.793387676,
  12462461.386566192,
  12462659.24080902,
  12463060.080401521,
  12463799.766469764,
  12464104.278917091,
  12464865.50202505,
  12465157.656546237,
  12465292.1875,
  12465684.375,
  12465959.748769395,
  12466034.807243405,
  12466044.15024325,
  12466137.5,
  12466221.875,
  12467146.756972058,
  12468006.498682642,
  12470824.582563516,
  12471129.6875,
  12473804.6875,
  12474046.875,
  12474610.9375,
  12474653.125,
  12477406.079069437,
  12477711.960540265,
  12477880.690797407,
  12477902.313297536,
  12478436.948230572,
  12482798.781729758,
  12482814.766896904,
  12482963.427293448,
  12482982.8125,
  12483013.539455382,
  12483075.977168985,
  12483080.980597656,
  12483116.010108437,
  12483134.941611752,
  12483142.1875,
  12483190.03370018,
  12483199.903074529,
  12483263.069600042,
  12483271.292567734,
  12483373.076325474,
  12483407.612917833,
  12483433.410456043,
  12483437.914202208,
  12483440.997874942,
  12483441.835304867,
  12483470.057304375,
  12483475.0,
  12483531.490946762,
  12483534.375,
  12483553.803175142,
  12483613.17324111,
  12483614.0625,
  12483614.252150819,
  12483625.945432298,
  12483626.5625,
  12483699.233020019,
  12483704.6875,
  12483709.375,
  12483716.299569217,
  12483717.1875,
  12483739.0625,
  12483750.0,
  12483802.927420411,
  12483808.967707543,
  12483839.0625,
  12483839.0625,
  12483875.897974633,
  12483889.693701714,
  12483922.533697639,
  12483933.90565027,
  12483949.15722256,
  12484039.670209765,
  12484048.4375,
  12484051.195774453,
  12484084.986598872,
  12484088.907235784,
  12484106.25,
  12484116.640465459,
  12484125.150760682,
  12484142.486732727,
  12484172.874919172,
  12484175.0,
  12484249.44943961,
  12484252.673398979,
  12484276.215081908,
  12484334.375,
  12484351.505218223,
  12484368.296546234,
  12484377.051737294,
  12484396.348905586,
  12484452.359739888,
  12484513.632896487,
  12484523.453744859,
  12484529.368383035,
  12484531.850962084,
  12484535.746577486,
  12484537.442165133,
  12484541.447840609,
  12484543.631709643,
  12484554.372816136,
  12484576.517665122,
  12484578.125,
  12484578.281292558,
  12484585.9375,
  12484589.0625,
  12484602.589748371,
  12484606.198046824,
  12484624.281988807,
  12484650.787377324,
  12484660.9375,
  12484675.752636196,
  12484676.521624032,
  12484678.728347644,
  12484679.423367925,
  12484681.671221549,
  12484690.97312798,
  12484704.6875,
  12484707.203258539,
  12484716.361892644,
  12484718.799264,
  12484733.400892474,
  12484762.5,
  12484763.852532241,
  12484791.783537438,
  12484803.125,
  12484833.165727066,
  12484841.065890003,
  12484842.668689597,
  12484843.75,
  12484864.0625,
  12484864.721153002,
  12484865.625,
  12484865.86928223,
  12484914.0625,
  12484920.3125,
  12484921.31379078,
  12484921.875,
  12484932.8125,
  12484934.25045108,
  12484940.870542074,
  12484946.751185931,
  12484947.693685168,
  12484951.317583993,
  12484962.606082259,
  12484987.477190815,
  12484988.530671118,
  12485000.654375043,
  12485010.133116776,
  12485012.5,
  12485021.773685321,
  12485033.204493854,
  12485040.906354176,
  12485044.034115646,
  12485045.259240624,
  12485072.388641825,
  12485109.375,
  12485132.670030052,
  12485138.156129232,
  12485149.949379055,
  12485153.935033027,
  12485159.3514224,
  12485160.825105151,
  12485162.5,
  12485174.850765664,
  12485181.25,
  12485182.573843576,
  12485187.5,
  12485189.0625,
  12485202.495316336,
  12485208.658545874,
  12485233.728667852,
  12485272.920024453,
  12485277.311232613,
  12485280.789746601,
  12485331.432438534,
  12485334.375,
  12485343.78695632,
  12485401.400514033,
  12485406.064010208,
  12485419.203550613,
  12485443.767725334,
  12485451.35982336,
  12485457.8125,
  12485470.3125,
  12485485.9375,
  12485511.646069713,
  12485521.155412925,
  12485554.929352881,
  12485570.178409105,
  12485573.4375,
  12485587.5,
  12485593.75,
  12485608.543041794,
  12485674.996756652,
  12485682.860724129,
  12485703.635479359,
  12485718.353978677,
  12485740.205124361,
  12485744.129773118,
  12485801.5625,
  12485804.97459279,
  12485810.893955808,
  12485811.014899192,
  12485824.77754257,
  12485830.6829768,
  12485870.914166197,
  12485899.37263224,
  12485906.24892023,
  12485907.907916911,
  12485914.409838758,
  12485917.1875,
  12485953.894877963,
  12485965.741874095,
  12486070.3125,
  12486078.97830789,
  12486079.898505686,
  12486112.332051605,
  12486131.25,
  12486190.159488682,
  12486190.366810668,
  12486201.130727619,
  12486254.6875,
  12486268.318494622,
  12486546.875,
  12486573.4375,
  12486633.783059942,
  12486809.962634478,
  12486896.875,
  12487112.271738444,
  12487461.508081421,
  12487511.743821198,
  12487782.071446456,
  12487783.011805061,
  12487837.744988529,
  12487839.0625,
  12487930.1051326,
  12488016.141207064,
  12488045.239556307,
  12488153.125,
  12488211.537450785,
  12488238.99854379,
  12488263.192916077,
  12488299.45941862,
  12488322.676409803,
  12488359.375,
  12488361.074303394,
  12488380.635380907,
  12488390.545051357,
  12488414.933630584,
  12488435.440400722,
  12488442.1875,
  12488495.273064727,
  12488496.875,
  12488510.59185944,
  12488544.17255787,
  12488567.524606531,
  12488588.327682339,
  12488608.64707329,
  12488635.9375,
  12488638.870789122,
  12488667.68441781,
  12488681.359977726,
  12488685.9375,
  12488721.771628581,
  12488746.071470276,
  12488753.7384552,
  12488770.759002756,
  12488816.571193267,
  12488860.083632132,
  12488957.8125,
  12488963.59457016,
  12488999.754093206,
  12489067.860075442,
  12489071.898554211,
  12489102.014430651,
  12489186.43841051,
  12489190.484312676,
  12489196.35716567,
  12489212.308330331,
  12489238.785092873,
  12489296.875,
  12489395.3125,
  12489441.442679115,
  12489512.796229681,
  12489514.0625,
  12489563.984592019,
  12489575.0,
  12489585.391687103,
  12489596.518509097,
  12489598.261168296,
  12489618.75,
  12489621.008741235,
  12489631.018182522,
  12489644.809926324,
  12489656.47743189,
  12489677.966998324,
  12489706.25,
  12489711.284182452,
  12489749.100217065,
  12489771.798574172,
  12489792.843319522,
  12489812.5,
  12489817.80145353,
  12489839.207011605,
  12489908.81881227,
  12489927.645674584,
  12489927.714997852,
  12489932.8125,
  12489960.9375,
  12490130.883691173,
  12490156.761285847,
  12490193.10247562,
  12490204.907194976,
  12490257.699484117,
  12490289.787799528,
  12490303.941453287,
  12490343.75,
  12490365.3529407,
  12490407.8125,
  12490438.5473968,
  12490521.02339622,
  12490593.814094542,
  12490627.453184415,
  12490670.132433543,
  12490683.75735773,
  12490686.518592007,
  12490737.119451202,
  12490820.022300342,
  12490876.217348158,
  12490919.734678585,
  12490931.25,
  12491022.662480893,
  12491156.25,
  12491162.009963434,
  12491190.625,
  12491257.8125,
  12491296.875,
  12491399.153855111,
  12491449.426244935,
  12491600.223154653,
  12491624.939001372,
  12491734.521849304,
  12491973.87556011,
  12492004.529601667,
  12492169.688616702,
  12492319.02989759,
  12492461.89282896,
  12492468.564718042,
  12492498.198679073,
  12492513.327596009,
  12492564.153597169,
  12492623.288618093,
  12492701.5625,
  12493039.011530682,
  12493847.03076862,
  12493918.350042826,
  12493926.216887426,
  12493959.204417396,
  12494196.44261006,
  12494491.889644405,
  12494646.765198428,
  12494701.420367543,
  12495587.958964828,
  12495743.75,
  12495754.568365885,
  12495832.56004037,
  12495894.115550946,
  12496513.52320837,
  12496560.558882074,
  12496692.476343887,
  12496884.375,
  12496895.3125,
  12496900.0,
  12497196.875,
  12497247.085662022,
  12497724.913023826,
  12497737.5,
  12498024.80860055,
  12498070.606707606,
  12498298.361328416,
  12498300.0,
  12498333.840802947,
  12498648.014628079,
  12498769.815556247,
  12498871.562471412,
  12499148.4375,
  12500088.826117525,
  12501196.50931726,
  12501992.18207627,
  12502064.374650987,
  12503242.597424798,
  12504282.8125,
  12506257.252408797,
  12506575.770622978,
  12506587.5,
  12507656.25,
  12508159.535052627,
  12508375.83123053,
  12509292.051979573,
  12509547.41517504,
  12511390.206322545,
  12516028.126519538,
  373853869.6446979,
  373853870.12083316,
  373853870.8804961,
  373853873.758095,
  373853901.7736579],
 [7.004302851715175,
  54.2576156490845,
  52.7827389588267,
  50.357055805034804,
  84.88691555930873,
  13.451005207967631,
  37.474570952328,
  48.31037786080955,
  70.93875680701878,
  17.36086264441774,
  92.87774448067515,
  7.688982170720358,
  6.793614087306694,
  44.60195710740117,
  7.5680691137891,
  51.26032352770066,
  52.775889813085335,
  21.43108277118567,
  16.640939237347535,
  35.035939337396606,
  36.6024684837862,
  18.209113547467805,
  43.31381393721329,
  88.67862736395894,
  6.363457327753961,
  66.0216006285424,
  8.180296070374004,
  13.026764143938262,
  34.156679043828504,
  19.75929476696141,
  25.77954790495339,
  32.780327284854856,
  12.878156688888136,
  44.877564915847245,
  7.587019555074294,
  10.911341887162749,
  72.88546053844883,
  38.95503365303942,
  48.2157197870433,
  16.068933879874905,
  20.58820386433075,
  54.09716674459274,
  10.627632353280555,
  150.66802807242107,
  15.352020608677297,
  16.283923237368274,
  8.311731724044057,
  12.13023167617171,
  72.92585432215718,
  19.723332358476174,
  9.887409569480281,
  54.62312045103624,
  18.840361418313126,
  53.059243698315264,
  28.984531191684297,
  92.98324569478643,
  87.54781679144352,
  13.311946311933655,
  196.2597482540032,
  143.58960223654444,
  74.24594720791205,
  18.997492205317556,
  39.05979334971462,
  12.581703743166399,
  10.947884422484599,
  47.44537166522119,
  39.65412861729904,
  26.769453105470635,
  58.741996208758316,
  58.039444158532234,
  34.21384982641598,
  5.884772675791998,
  16.01345310123795,
  56.18019973872723,
  41.80077085039523,
  52.7835127493145,
  42.874367104633464,
  75.39338216402459,
  57.381499354199875,
  82.75570976709292,
  30.814960507178643,
  20.618730001447922,
  13.557536093316953,
  84.07985093763213,
  7.1640582108129,
  71.49154059157988,
  7.1972412858931705,
  74.37991624126794,
  50.9426347706555,
  26.96445422676303,
  13.799568108258393,
  29.654766555404105,
  11.394325189897078,
  99.40537091712687,
  72.169510322715,
  13.893637795889227,
  19.981690038911612,
  16.292650661737497,
  19.178135176252784,
  55.903287662091294,
  40.087610994720364,
  6.051922579601536,
  18.380163276680474,
  29.833401939009743,
  64.75096314754072,
  41.69102500136764,
  70.20106682384063,
  10.372661754478065,
  20.328089075063424,
  39.78241579450581,
  7.117603114454929,
  8.912957808242266,
  7.56097423474301,
  6.741273434609248,
  157.5087950932332,
  38.502386929481354,
  47.90003321755987,
  12.686550200455416,
  9.98397618907522,
  73.5716276338428,
  7.293245638372215,
  18.868770062679392,
  17.40085327646871,
  13.391151812794407,
  76.95218041428222,
  143.84704343640388,
  21.019144970879033,
  119.70922790299124,
  24.79997120257438,
  11.987920190635021,
  47.51887547267849,
  54.846988076010426,
  35.26507073396717,
  64.42799532131818,
  87.17265577492589,
  45.72133186257568,
  50.342673504470895,
  54.75249330360912,
  16.414537058338713,
  33.913248319500404,
  12.200874101335492,
  31.701004971059195,
  5.106881428560929,
  25.378212555741918,
  145.50320410519132,
  31.20017081253198,
  16.783007181597725,
  7.22562359797847,
  22.942040776782378,
  91.59286340540658,
  5.381094068944606,
  14.53693168268741,
  61.05627494577508,
  56.30043106973886,
  23.92231310381323,
  5.177309726449087,
  52.91093571436809,
  122.26197683506172,
  39.60483104421091,
  21.435926330469293,
  43.4313252792459,
  17.676044869099616,
  43.51183254410058,
  6.243795761382695,
  31.05185256128948,
  50.585746674960376,
  23.92848960483153,
  62.18069239650571,
  82.51481921657395,
  77.05730428469288,
  94.48178940736011,
  48.11122034051288,
  34.74683324295322,
  47.18362293916658,
  18.487298006451784,
  92.5403351761172,
  7.642675884445015,
  16.543419812152752,
  44.65363405110824,
  78.68769895852688,
  55.11437889409917,
  33.731281491752995,
  9.420985272507867,
  39.54605405314805,
  41.56544720619895,
  13.257328290267752,
  9.876691390291787,
  20.34140535493127,
  7.723159920603721,
  14.541343412468777,
  17.783245530486912,
  18.289647394742786,
  9.970409553460744,
  55.846599662397956,
  9.959335835894477,
  50.64429440312175,
  26.403157272545364,
  12.394371238955232,
  78.96261134634669,
  130.59362756550067,
  13.148077411442976,
  70.05581898248795,
  10.193893022464694,
  96.32368503245443,
  7.749920382279488,
  36.99474452122499,
  13.024135217484517,
  49.613292408978396,
  5.132476107843655,
  7.619006009319628,
  5.9794908351563265,
  14.453569057794926,
  167.66465604011847,
  36.459816197461514,
  21.668572589839002,
  29.199583527466935,
  26.96100048679822,
  44.30584305497363,
  43.36511135554704,
  11.766920959404809,
  16.89203317410678,
  16.723368278439867,
  38.29377460664038,
  79.41469961341296,
  15.880273033591848,
  30.467718351745688,
  22.358247053487638,
  31.71109269081241,
  53.34441290664619,
  59.06453587539118,
  59.63862412631337,
  34.535156214472444,
  14.327709857483987,
  47.17062367592397,
  101.99609410893586,
  18.860771961106515,
  20.99176312790026,
  27.327591941047455,
  18.536447383394865,
  92.14102266984797,
  71.342591476699,
  13.263372264995027,
  50.776525750768364,
  86.13653436808794,
  16.645164892386383,
  64.83821524928753,
  72.52555679904937,
  25.235613953484524,
  14.584511118508702,
  8.612655771167354,
  94.82866937456829,
  75.72663622248855,
  33.43178986384958,
  29.536987563039947,
  29.206621540734183,
  24.177563463963,
  17.612723189326772,
  24.340439120916738,
  94.22795613391362,
  9.543800394052464,
  58.82793264271234,
  94.47109449229434,
  5.220920218424627,
  82.83838358595963,
  8.389480426354208,
  30.34547082245505,
  69.86694020126463,
  53.42303347536786,
  52.550420806932244,
  46.18233416077441,
  6.906090549187321,
  51.157083810099515,
  42.80377581257708,
  34.40652033329609,
  9.169745327021296,
  13.587329868822412,
  55.087905846047995,
  52.88516976655463,
  123.03580849722142,
  9.928969637227564,
  5.750002162543088,
  28.25382798963022,
  26.33959711933607,
  12.232059925358381,
  72.18310174586425,
  6.050775552412595,
  9.908268129065732,
  108.26000720346725,
  54.28341322582772,
  21.03639537707841,
  19.799223952626363,
  77.44624123492947,
  37.59402633192926,
  73.42979569055228,
  69.12461772150091,
  13.460583203859652,
  63.24880513424552,
  45.98950820977657,
  65.55827183899794,
  21.970339581547233,
  11.989449085009378,
  21.236684927961594,
  174.409413670405,
  107.92047303779387,
  141.05378729203525,
  74.63349848735584,
  17.863132382222727,
  25.251780797547397,
  11.03849469217633,
  10.139757204859007,
  26.38701041786543,
  10.669054229905635,
  27.667725661070612,
  55.23129586761987,
  12.556336105944542,
  68.2992782217281,
  36.66744423772254,
  41.17865015946295,
  27.537142937699215,
  119.90196614019536,
  19.12708698036304,
  65.63377074346508,
  5.6570844134998675,
  57.30944142408754,
  17.8779411943219,
  125.64664343840886,
  28.409435250080143,
  28.764762739705642,
  93.55518154446904,
  7.988695433015603,
  153.09994430758638,
  73.61869899666469,
  9.883958342222924,
  89.81762025170687,
  24.97319244343395,
  10.549999342514699,
  67.27213341003649,
  26.6067678600676,
  27.084375626173784,
  9.774357256048104,
  122.67576648572899,
  52.04371866078608,
  11.689104752251565,
  41.16667644221455,
  78.74668938720963,
  100.06660526156661,
  73.68786333515878,
  28.854773850467076,
  34.19099954820316,
  29.28381452393699,
  47.88492508774891,
  14.29782781489353,
  38.619368181779244,
  12.704465547437419,
  78.3327196547419,
  65.67829758080144,
  116.62608577174841,
  8.918852869246702,
  9.953641023560472,
  55.56046798900046,
  67.95776145625702,
  26.67033309422833,
  14.735778190697665,
  44.44708721087042,
  11.60286912011437,
  59.979078412341295,
  5.740246385936458,
  72.16279462989951,
  25.9956460623521,
  8.666961055875104,
  27.56147093067576,
  6.622038386396409,
  28.066581979028847,
  85.8965392046531,
  32.507970542814704,
  53.75308184294728,
  61.99086464273856,
  43.13756800121971,
  33.42272652884848,
  42.491834194569535,
  12.984372213600585,
  109.16699091385988,
  7.972004532452814,
  115.95675369980508,
  31.60747703140596,
  8.970362227862847,
  5.348212292398322,
  170.32447438112524,
  25.39461971029664,
  25.689753420439597,
  128.4108574894461,
  22.923517291466485,
  23.879023214921144,
  32.850704516803845,
  94.241834519035,
  47.674956807162374,
  42.874078537524895,
  267.69045911495505,
  26.979484728330867,
  7.693118407566914,
  26.129684498104048,
  79.94621217989265,
  72.29553750566461,
  50.67407470133528,
  11.158528124498075,
  82.79205761535891,
  19.157737790208685,
  17.32440397627579,
  20.555302571995973,
  10.529247429324672,
  16.617831539622564,
  79.30467464529444,
  9.586556529180305,
  27.691560468355817,
  15.244960173626433,
  14.156378814162549,
  28.312833940843245,
  144.42400385897858,
  30.843951686429108,
  28.503842966598675,
  33.73055398605692,
  11.717792275688392,
  59.730785744413424,
  98.7516147375063,
  27.017210944780526,
  72.69447911598496,
  29.387392340651942,
  55.48065694190197,
  25.598014438464617,
  56.37725560670995,
  10.470496341939722,
  10.746776600343148,
  10.513823457990869,
  48.16220169554187,
  88.87974662269184,
  43.000628041814466,
  45.45810061817287,
  10.313405625864716,
  5.8587784221719765,
  56.65108657998005,
  63.16020871410067,
  5.744059779469008,
  19.26888098624897,
  12.711023549112907,
  33.50586661377136,
  5.171106749425389,
  83.62372887334072,
  5.076570316611948,
  9.80967873342721,
  32.90073091688129,
  71.11619460163762,
  7.83566819401059,
  135.74274524315902,
  97.51342640508696,
  13.21626051679822,
  19.231859023579307,
  86.74391518532043,
  18.378008196723776,
  70.00226757629704,
  27.03744922006843,
  65.00824796111104,
  81.69910871291269,
  150.24029501106264,
  13.5078216913279,
  87.52256073153305,
  41.99929853602795,
  16.840500379771875,
  16.40353747251429,
  12.725871531799617,
  24.556200229844272,
  38.991822233074046,
  9.1743961229883,
  12.004724580251052,
  11.30407325217084,
  28.82335480698434,
  76.9413871671114,
  12.565114759807928,
  12.244217241795173,
  16.968679101994923,
  10.284049698143138,
  5.568490345122137,
  38.709524829487066,
  13.414745278483165,
  24.467206737540007,
  18.32227261747431,
  9.351399319684976,
  16.998101233912426,
  7.259660085879079,
  94.82048750778833,
  26.50402712375147,
  60.45146088997449,
  24.04368896093938,
  34.22015423624592,
  49.02417389527286,
  12.70207873920776,
  6.824123533046316,
  60.22696514383394,
  59.03014922333422,
  37.33383761241198,
  33.529222118726295,
  23.055007665804016,
  25.470352667080007,
  40.60119960686165,
  18.480932845630722,
  24.662584614048853,
  21.467464795378213,
  24.434821051105246,
  25.75149895464371,
  28.708874142998877,
  11.71658061440018,
  8.67115210440233,
  5.701911266958268,
  33.470786470787885,
  103.39594277767675,
  129.0272434588124,
  6.103257190964134,
  53.80152999506884,
  15.18107246135871,
  13.57947767131657,
  83.12847209064773,
  37.91035844579114,
  5.035424514497084,
  90.07608276274134,
  20.137724540153233,
  15.850670565731386,
  104.74351913363358,
  54.972242220134206,
  19.345840830152646,
  52.315840961158194,
  21.827068816322328,
  69.80099229568637,
  13.284871379796368,
  12.456256639028917,
  11.229390706034335,
  12.878889905664694,
  14.296897577943332,
  77.82817288487384,
  6.582211199591967,
  13.665549855630045,
  77.7139704877213,
  9.055668241211315,
  5.772428704689789,
  71.80452160687348,
  17.060298215282053,
  101.57377530970028,
  5.676264744774738,
  62.753506469608396,
  74.6317907709306,
  47.66046472211516,
  12.468472489162199,
  5.772712952145254,
  15.866236108754977,
  5.289206764185928,
  9.45873791405088,
  63.797398144249144,
  110.19268323613245,
  89.99823758314976,
  16.12767166966444,
  20.486722472118874,
  25.33573871581273,
  82.96305957983394,
  34.740123416116816,
  8.335345774516425,
  7.550262830344939,
  63.97387455254626,
  11.048962226756503,
  5.622720151163016,
  22.552717093928834,
  38.902615652346185,
  5.80146551371251,
  13.515108358102598,
  21.545861529608064,
  21.17567992026929,
  19.985889717332828,
  33.8836065728687,
  67.14101480271866,
  16.26312299478976,
  24.377611928969046,
  5.173963365721217,
  13.219708422485201,
  22.795897220886978,
  9.588192821459481,
  60.719493319813644,
  31.887376943856992,
  61.63608914382314,
  72.03500696873377,
  45.43001952559559,
  9.258926999080897,
  38.95530840917322,
  10.841160035606705,
  9.593211240830037,
  16.7112933995029,
  43.773150754746744,
  21.158617058540894,
  69.0223895100502,
  63.2394972887217,
  5.596985157989508,
  87.56914860691914,
  58.911324525149524,
  13.70960487998605,
  14.920876919768775,
  6.121173059664683,
  6.131791039219603,
  91.84862022944479,
  5.43822957912523,
  8.740972912735097,
  68.20349774504608,
  45.2005446764171,
  15.588788948110123,
  90.53629904683305,
  5.927689844767063,
  17.496416889395206,
  34.3152364079604,
  7.1218511547982555,
  354.39322684690893,
  167.89540900358708,
  17.66175664818565,
  18.934687567299672,
  158.36026077864338])

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

([4377943.284777576,
  12255914.0625,
  12256329.440357525,
  12350185.856029617,
  12350435.9375,
  12350455.795000589,
  12354271.41548004,
  12355335.9375,
  12355503.667421496,
  12357477.079338957,
  12357818.614609193,
  12358233.043507664,
  12358295.461459577,
  12358806.25,
  12358926.464205928,
  12359014.0625,
  12359614.0625,
  12359625.573930643,
  12359647.259285161,
  12361462.5,
  12362081.25,
  12363108.426152745,
  12364302.568847135,
  12364984.375,
  12365433.53820789,
  12366037.5,
  12366396.668463694,
  12367611.886960475,
  12367621.875,
  12367638.369155396,
  12367762.343800675,
  12368516.549478827,
  12368541.515387911,
  12370076.5625,
  12370356.204935547,
  12370868.607021285,
  12373801.034123015,
  12373909.375,
  12374256.25,
  12374508.079869175,
  12374638.46673763,
  12374809.375,
  12374815.279637367,
  12375015.215434952,
  12375051.384495022,
  12375280.00012539,
  12375770.359766342,
  12375889.62138288,
  12375995.3125,
  12376019.614683805,
  12376542.779199434,
  12376676.212086761,
  12376758.560438618,
  12377829.6875,
  12377854.6875,
  12378148.085868357,
  12378183.181248715,
  12378411.063585557,
  12378417.581000097,
  12378618.638100153,
  12378670.148967575,
  12378772.644742787,
  12378923.4375,
  12378938.26618486,
  12378976.44432619,
  12379085.246972604,
  12379132.123131301,
  12379310.842022989,
  12379349.519611906,
  12379478.271606395,
  12379509.959710054,
  12379799.350521581,
  12380064.317919713,
  12380101.5625,
  12380143.75,
  12380200.0,
  12380218.75,
  12380274.852569217,
  12380294.331618657,
  12381185.9375,
  12381279.6875,
  12381288.470141595,
  12381313.644901635,
  12381318.75,
  12381426.722140273,
  12381432.021509686,
  12381864.343178418,
  12382220.3125,
  12382343.75,
  12382349.039412165,
  12382441.851744793,
  12382465.867682574,
  12382475.164130617,
  12382577.896651376,
  12382610.805211566,
  12382774.694935815,
  12382960.483284507,
  12383078.439428175,
  12383258.472484248,
  12383308.977514647,
  12383321.652865456,
  12383435.582282286,
  12383533.493772384,
  12383562.704805171,
  12383843.05684846,
  12384185.9375,
  12384298.655672044,
  12384396.331523027,
  12385641.562205112,
  12385953.125,
  12386310.92283924,
  12386329.147376813,
  12386360.819116259,
  12387310.825105557,
  12387451.47191883,
  12387648.4375,
  12388776.5625,
  12389169.561119694,
  12389481.932590725,
  12389542.1875,
  12389628.740562841,
  12389871.064415647,
  12390442.721593175,
  12390457.336101025,
  12390826.210399674,
  12393743.75,
  12393751.355175387,
  12393857.8125,
  12394288.840428794,
  12394389.550966877,
  12394433.437107764,
  12394705.25648869,
  12394729.6875,
  12394862.5,
  12395015.0445802,
  12395199.845917802,
  12395277.452527585,
  12395645.3125,
  12395672.214333445,
  12395982.8125,
  12396072.46789041,
  12396078.125,
  12396109.798390329,
  12396301.380054422,
  12396313.64245891,
  12396635.9375,
  12396658.224157428,
  12396751.944499362,
  12396796.508937938,
  12396881.25,
  12397480.787809664,
  12397582.7021443,
  12397798.4375,
  12397985.9375,
  12398543.49004211,
  12398870.818241073,
  12399457.8125,
  12399552.3246644,
  12399707.8125,
  12401185.927070305,
  12401365.61845149,
  12404167.073746992,
  12404339.319597904,
  12404643.562496297,
  12405013.269530231,
  12405142.1875,
  12406502.663532255,
  12407585.9375,
  12408482.8125,
  12410842.1875,
  12414386.349441666,
  12416169.084054116,
  12417307.776898224,
  12418067.144471522,
  12444225.697179591,
  12444520.810785094,
  12444631.357740153,
  12445208.784452088,
  12445995.3125,
  12449925.584971754,
  12451404.6875,
  12451457.8125,
  12452351.66657799,
  12455740.625,
  12455882.672177993,
  12456330.724715944,
  12457113.081415823,
  12457602.50124031,
  12458931.56982687,
  12459931.140894866,
  12460046.734162364,
  12460122.542642342,
  12460406.920009531,
  12460532.8125,
  12460666.284695176,
  12460679.6875,
  12460686.068091732,
  12460816.307987213,
  12460976.5625,
  12460999.144124873,
  12461115.465272877,
  12461380.304726053,
  12461489.932332514,
  12461561.080045262,
  12461745.70062896,
  12461885.600787194,
  12461991.573433017,
  12462043.193655215,
  12462136.029631585,
  12462443.793387676,
  12462461.386566192,
  12462659.24080902,
  12463060.080401521,
  12463799.766469764,
  12464104.278917091,
  12464865.50202505,
  12465157.656546237,
  12465292.1875,
  12465684.375,
  12465959.748769395,
  12466034.807243405,
  12466044.15024325,
  12466137.5,
  12466221.875,
  12467146.756972058,
  12468006.498682642,
  12470824.582563516,
  12471129.6875,
  12473804.6875,
  12474046.875,
  12474610.9375,
  12474653.125,
  12477406.079069437,
  12477711.960540265,
  12477880.690797407,
  12477902.313297536,
  12478436.948230572,
  12482798.781729758,
  12482814.766896904,
  12482963.427293448,
  12482982.8125,
  12483013.539455382,
  12483075.977168985,
  12483080.980597656,
  12483116.010108437,
  12483134.941611752,
  12483142.1875,
  12483190.03370018,
  12483199.903074529,
  12483263.069600042,
  12483271.292567734,
  12483373.076325474,
  12483407.612917833,
  12483433.410456043,
  12483437.914202208,
  12483440.997874942,
  12483441.835304867,
  12483470.057304375,
  12483475.0,
  12483531.490946762,
  12483534.375,
  12483553.803175142,
  12483613.17324111,
  12483614.0625,
  12483614.252150819,
  12483625.945432298,
  12483626.5625,
  12483699.233020019,
  12483704.6875,
  12483709.375,
  12483716.299569217,
  12483717.1875,
  12483739.0625,
  12483750.0,
  12483802.927420411,
  12483808.967707543,
  12483839.0625,
  12483839.0625,
  12483875.897974633,
  12483889.693701714,
  12483922.533697639,
  12483933.90565027,
  12483949.15722256,
  12484039.670209765,
  12484048.4375,
  12484051.195774453,
  12484084.986598872,
  12484088.907235784,
  12484106.25,
  12484116.640465459,
  12484125.150760682,
  12484142.486732727,
  12484172.874919172,
  12484175.0,
  12484249.44943961,
  12484252.673398979,
  12484276.215081908,
  12484334.375,
  12484351.505218223,
  12484368.296546234,
  12484377.051737294,
  12484396.348905586,
  12484452.359739888,
  12484513.632896487,
  12484523.453744859,
  12484529.368383035,
  12484531.850962084,
  12484535.746577486,
  12484537.442165133,
  12484541.447840609,
  12484543.631709643,
  12484554.372816136,
  12484576.517665122,
  12484578.125,
  12484578.281292558,
  12484585.9375,
  12484589.0625,
  12484602.589748371,
  12484606.198046824,
  12484624.281988807,
  12484650.787377324,
  12484660.9375,
  12484675.752636196,
  12484676.521624032,
  12484678.728347644,
  12484679.423367925,
  12484681.671221549,
  12484690.97312798,
  12484704.6875,
  12484707.203258539,
  12484716.361892644,
  12484718.799264,
  12484733.400892474,
  12484762.5,
  12484763.852532241,
  12484791.783537438,
  12484803.125,
  12484833.165727066,
  12484841.065890003,
  12484842.668689597,
  12484843.75,
  12484864.0625,
  12484864.721153002,
  12484865.625,
  12484865.86928223,
  12484914.0625,
  12484920.3125,
  12484921.31379078,
  12484921.875,
  12484932.8125,
  12484934.25045108,
  12484940.870542074,
  12484946.751185931,
  12484947.693685168,
  12484951.317583993,
  12484962.606082259,
  12484987.477190815,
  12484988.530671118,
  12485000.654375043,
  12485010.133116776,
  12485012.5,
  12485021.773685321,
  12485033.204493854,
  12485040.906354176,
  12485044.034115646,
  12485045.259240624,
  12485072.388641825,
  12485109.375,
  12485132.670030052,
  12485138.156129232,
  12485149.949379055,
  12485153.935033027,
  12485159.3514224,
  12485160.825105151,
  12485162.5,
  12485174.850765664,
  12485181.25,
  12485182.573843576,
  12485187.5,
  12485189.0625,
  12485202.495316336,
  12485208.658545874,
  12485233.728667852,
  12485272.920024453,
  12485277.311232613,
  12485280.789746601,
  12485331.432438534,
  12485334.375,
  12485343.78695632,
  12485401.400514033,
  12485406.064010208,
  12485419.203550613,
  12485443.767725334,
  12485451.35982336,
  12485457.8125,
  12485470.3125,
  12485485.9375,
  12485511.646069713,
  12485521.155412925,
  12485554.929352881,
  12485570.178409105,
  12485573.4375,
  12485587.5,
  12485593.75,
  12485608.543041794,
  12485674.996756652,
  12485682.860724129,
  12485703.635479359,
  12485718.353978677,
  12485740.205124361,
  12485744.129773118,
  12485801.5625,
  12485804.97459279,
  12485810.893955808,
  12485811.014899192,
  12485824.77754257,
  12485830.6829768,
  12485870.914166197,
  12485899.37263224,
  12485906.24892023,
  12485907.907916911,
  12485914.409838758,
  12485917.1875,
  12485953.894877963,
  12485965.741874095,
  12486070.3125,
  12486078.97830789,
  12486079.898505686,
  12486112.332051605,
  12486131.25,
  12486190.159488682,
  12486190.366810668,
  12486201.130727619,
  12486254.6875,
  12486268.318494622,
  12486546.875,
  12486573.4375,
  12486633.783059942,
  12486809.962634478,
  12486896.875,
  12487112.271738444,
  12487461.508081421,
  12487511.743821198,
  12487782.071446456,
  12487783.011805061,
  12487837.744988529,
  12487839.0625,
  12487930.1051326,
  12488016.141207064,
  12488045.239556307,
  12488153.125,
  12488211.537450785,
  12488238.99854379,
  12488263.192916077,
  12488299.45941862,
  12488322.676409803,
  12488359.375,
  12488361.074303394,
  12488380.635380907,
  12488390.545051357,
  12488414.933630584,
  12488435.440400722,
  12488442.1875,
  12488495.273064727,
  12488496.875,
  12488510.59185944,
  12488544.17255787,
  12488567.524606531,
  12488588.327682339,
  12488608.64707329,
  12488635.9375,
  12488638.870789122,
  12488667.68441781,
  12488681.359977726,
  12488685.9375,
  12488721.771628581,
  12488746.071470276,
  12488753.7384552,
  12488770.759002756,
  12488816.571193267,
  12488860.083632132,
  12488957.8125,
  12488963.59457016,
  12488999.754093206,
  12489067.860075442,
  12489071.898554211,
  12489102.014430651,
  12489186.43841051,
  12489190.484312676,
  12489196.35716567,
  12489212.308330331,
  12489238.785092873,
  12489296.875,
  12489395.3125,
  12489441.442679115,
  12489512.796229681,
  12489514.0625,
  12489563.984592019,
  12489575.0,
  12489585.391687103,
  12489596.518509097,
  12489598.261168296,
  12489618.75,
  12489621.008741235,
  12489631.018182522,
  12489644.809926324,
  12489656.47743189,
  12489677.966998324,
  12489706.25,
  12489711.284182452,
  12489749.100217065,
  12489771.798574172,
  12489792.843319522,
  12489812.5,
  12489817.80145353,
  12489839.207011605,
  12489908.81881227,
  12489927.645674584,
  12489927.714997852,
  12489932.8125,
  12489960.9375,
  12490130.883691173,
  12490156.761285847,
  12490193.10247562,
  12490204.907194976,
  12490257.699484117,
  12490289.787799528,
  12490303.941453287,
  12490343.75,
  12490365.3529407,
  12490407.8125,
  12490438.5473968,
  12490521.02339622,
  12490593.814094542,
  12490627.453184415,
  12490670.132433543,
  12490683.75735773,
  12490686.518592007,
  12490737.119451202,
  12490820.022300342,
  12490876.217348158,
  12490919.734678585,
  12490931.25,
  12491022.662480893,
  12491156.25,
  12491162.009963434,
  12491190.625,
  12491257.8125,
  12491296.875,
  12491399.153855111,
  12491449.426244935,
  12491600.223154653,
  12491624.939001372,
  12491734.521849304,
  12491973.87556011,
  12492004.529601667,
  12492169.688616702,
  12492319.02989759,
  12492461.89282896,
  12492468.564718042,
  12492498.198679073,
  12492513.327596009,
  12492564.153597169,
  12492623.288618093,
  12492701.5625,
  12493039.011530682,
  12493847.03076862,
  12493918.350042826,
  12493926.216887426,
  12493959.204417396,
  12494196.44261006,
  12494491.889644405,
  12494646.765198428,
  12494701.420367543,
  12495587.958964828,
  12495743.75,
  12495754.568365885,
  12495832.56004037,
  12495894.115550946,
  12496513.52320837,
  12496560.558882074,
  12496692.476343887,
  12496884.375,
  12496895.3125,
  12496900.0,
  12497196.875,
  12497247.085662022,
  12497724.913023826,
  12497737.5,
  12498024.80860055,
  12498070.606707606,
  12498298.361328416,
  12498300.0,
  12498333.840802947,
  12498648.014628079,
  12498769.815556247,
  12498871.562471412,
  12499148.4375,
  12500088.826117525,
  12501196.50931726,
  12501992.18207627,
  12502064.374650987,
  12503242.597424798,
  12504282.8125,
  12506257.252408797,
  12506575.770622978,
  12506587.5,
  12507656.25,
  12508159.535052627,
  12508375.83123053,
  12509292.051979573,
  12509547.41517504,
  12511390.206322545,
  12516028.126519538,
  373853869.6446979,
  373853870.12083316,
  373853870.8804961,
  373853873.758095,
  373853901.7736579],
 [7.004302851715175,
  54.2576156490845,
  52.7827389588267,
  50.357055805034804,
  84.88691555930873,
  13.451005207967631,
  37.474570952328,
  48.31037786080955,
  70.93875680701878,
  17.36086264441774,
  92.87774448067515,
  7.688982170720358,
  6.793614087306694,
  44.60195710740117,
  7.5680691137891,
  51.26032352770066,
  52.775889813085335,
  21.43108277118567,
  16.640939237347535,
  35.035939337396606,
  36.6024684837862,
  18.209113547467805,
  43.31381393721329,
  88.67862736395894,
  6.363457327753961,
  66.0216006285424,
  8.180296070374004,
  13.026764143938262,
  34.156679043828504,
  19.75929476696141,
  25.77954790495339,
  32.780327284854856,
  12.878156688888136,
  44.877564915847245,
  7.587019555074294,
  10.911341887162749,
  72.88546053844883,
  38.95503365303942,
  48.2157197870433,
  16.068933879874905,
  20.58820386433075,
  54.09716674459274,
  10.627632353280555,
  150.66802807242107,
  15.352020608677297,
  16.283923237368274,
  8.311731724044057,
  12.13023167617171,
  72.92585432215718,
  19.723332358476174,
  9.887409569480281,
  54.62312045103624,
  18.840361418313126,
  53.059243698315264,
  28.984531191684297,
  92.98324569478643,
  87.54781679144352,
  13.311946311933655,
  196.2597482540032,
  143.58960223654444,
  74.24594720791205,
  18.997492205317556,
  39.05979334971462,
  12.581703743166399,
  10.947884422484599,
  47.44537166522119,
  39.65412861729904,
  26.769453105470635,
  58.741996208758316,
  58.039444158532234,
  34.21384982641598,
  5.884772675791998,
  16.01345310123795,
  56.18019973872723,
  41.80077085039523,
  52.7835127493145,
  42.874367104633464,
  75.39338216402459,
  57.381499354199875,
  82.75570976709292,
  30.814960507178643,
  20.618730001447922,
  13.557536093316953,
  84.07985093763213,
  7.1640582108129,
  71.49154059157988,
  7.1972412858931705,
  74.37991624126794,
  50.9426347706555,
  26.96445422676303,
  13.799568108258393,
  29.654766555404105,
  11.394325189897078,
  99.40537091712687,
  72.169510322715,
  13.893637795889227,
  19.981690038911612,
  16.292650661737497,
  19.178135176252784,
  55.903287662091294,
  40.087610994720364,
  6.051922579601536,
  18.380163276680474,
  29.833401939009743,
  64.75096314754072,
  41.69102500136764,
  70.20106682384063,
  10.372661754478065,
  20.328089075063424,
  39.78241579450581,
  7.117603114454929,
  8.912957808242266,
  7.56097423474301,
  6.741273434609248,
  157.5087950932332,
  38.502386929481354,
  47.90003321755987,
  12.686550200455416,
  9.98397618907522,
  73.5716276338428,
  7.293245638372215,
  18.868770062679392,
  17.40085327646871,
  13.391151812794407,
  76.95218041428222,
  143.84704343640388,
  21.019144970879033,
  119.70922790299124,
  24.79997120257438,
  11.987920190635021,
  47.51887547267849,
  54.846988076010426,
  35.26507073396717,
  64.42799532131818,
  87.17265577492589,
  45.72133186257568,
  50.342673504470895,
  54.75249330360912,
  16.414537058338713,
  33.913248319500404,
  12.200874101335492,
  31.701004971059195,
  5.106881428560929,
  25.378212555741918,
  145.50320410519132,
  31.20017081253198,
  16.783007181597725,
  7.22562359797847,
  22.942040776782378,
  91.59286340540658,
  5.381094068944606,
  14.53693168268741,
  61.05627494577508,
  56.30043106973886,
  23.92231310381323,
  5.177309726449087,
  52.91093571436809,
  122.26197683506172,
  39.60483104421091,
  21.435926330469293,
  43.4313252792459,
  17.676044869099616,
  43.51183254410058,
  6.243795761382695,
  31.05185256128948,
  50.585746674960376,
  23.92848960483153,
  62.18069239650571,
  82.51481921657395,
  77.05730428469288,
  94.48178940736011,
  48.11122034051288,
  34.74683324295322,
  47.18362293916658,
  18.487298006451784,
  92.5403351761172,
  7.642675884445015,
  16.543419812152752,
  44.65363405110824,
  78.68769895852688,
  55.11437889409917,
  33.731281491752995,
  9.420985272507867,
  39.54605405314805,
  41.56544720619895,
  13.257328290267752,
  9.876691390291787,
  20.34140535493127,
  7.723159920603721,
  14.541343412468777,
  17.783245530486912,
  18.289647394742786,
  9.970409553460744,
  55.846599662397956,
  9.959335835894477,
  50.64429440312175,
  26.403157272545364,
  12.394371238955232,
  78.96261134634669,
  130.59362756550067,
  13.148077411442976,
  70.05581898248795,
  10.193893022464694,
  96.32368503245443,
  7.749920382279488,
  36.99474452122499,
  13.024135217484517,
  49.613292408978396,
  5.132476107843655,
  7.619006009319628,
  5.9794908351563265,
  14.453569057794926,
  167.66465604011847,
  36.459816197461514,
  21.668572589839002,
  29.199583527466935,
  26.96100048679822,
  44.30584305497363,
  43.36511135554704,
  11.766920959404809,
  16.89203317410678,
  16.723368278439867,
  38.29377460664038,
  79.41469961341296,
  15.880273033591848,
  30.467718351745688,
  22.358247053487638,
  31.71109269081241,
  53.34441290664619,
  59.06453587539118,
  59.63862412631337,
  34.535156214472444,
  14.327709857483987,
  47.17062367592397,
  101.99609410893586,
  18.860771961106515,
  20.99176312790026,
  27.327591941047455,
  18.536447383394865,
  92.14102266984797,
  71.342591476699,
  13.263372264995027,
  50.776525750768364,
  86.13653436808794,
  16.645164892386383,
  64.83821524928753,
  72.52555679904937,
  25.235613953484524,
  14.584511118508702,
  8.612655771167354,
  94.82866937456829,
  75.72663622248855,
  33.43178986384958,
  29.536987563039947,
  29.206621540734183,
  24.177563463963,
  17.612723189326772,
  24.340439120916738,
  94.22795613391362,
  9.543800394052464,
  58.82793264271234,
  94.47109449229434,
  5.220920218424627,
  82.83838358595963,
  8.389480426354208,
  30.34547082245505,
  69.86694020126463,
  53.42303347536786,
  52.550420806932244,
  46.18233416077441,
  6.906090549187321,
  51.157083810099515,
  42.80377581257708,
  34.40652033329609,
  9.169745327021296,
  13.587329868822412,
  55.087905846047995,
  52.88516976655463,
  123.03580849722142,
  9.928969637227564,
  5.750002162543088,
  28.25382798963022,
  26.33959711933607,
  12.232059925358381,
  72.18310174586425,
  6.050775552412595,
  9.908268129065732,
  108.26000720346725,
  54.28341322582772,
  21.03639537707841,
  19.799223952626363,
  77.44624123492947,
  37.59402633192926,
  73.42979569055228,
  69.12461772150091,
  13.460583203859652,
  63.24880513424552,
  45.98950820977657,
  65.55827183899794,
  21.970339581547233,
  11.989449085009378,
  21.236684927961594,
  174.409413670405,
  107.92047303779387,
  141.05378729203525,
  74.63349848735584,
  17.863132382222727,
  25.251780797547397,
  11.03849469217633,
  10.139757204859007,
  26.38701041786543,
  10.669054229905635,
  27.667725661070612,
  55.23129586761987,
  12.556336105944542,
  68.2992782217281,
  36.66744423772254,
  41.17865015946295,
  27.537142937699215,
  119.90196614019536,
  19.12708698036304,
  65.63377074346508,
  5.6570844134998675,
  57.30944142408754,
  17.8779411943219,
  125.64664343840886,
  28.409435250080143,
  28.764762739705642,
  93.55518154446904,
  7.988695433015603,
  153.09994430758638,
  73.61869899666469,
  9.883958342222924,
  89.81762025170687,
  24.97319244343395,
  10.549999342514699,
  67.27213341003649,
  26.6067678600676,
  27.084375626173784,
  9.774357256048104,
  122.67576648572899,
  52.04371866078608,
  11.689104752251565,
  41.16667644221455,
  78.74668938720963,
  100.06660526156661,
  73.68786333515878,
  28.854773850467076,
  34.19099954820316,
  29.28381452393699,
  47.88492508774891,
  14.29782781489353,
  38.619368181779244,
  12.704465547437419,
  78.3327196547419,
  65.67829758080144,
  116.62608577174841,
  8.918852869246702,
  9.953641023560472,
  55.56046798900046,
  67.95776145625702,
  26.67033309422833,
  14.735778190697665,
  44.44708721087042,
  11.60286912011437,
  59.979078412341295,
  5.740246385936458,
  72.16279462989951,
  25.9956460623521,
  8.666961055875104,
  27.56147093067576,
  6.622038386396409,
  28.066581979028847,
  85.8965392046531,
  32.507970542814704,
  53.75308184294728,
  61.99086464273856,
  43.13756800121971,
  33.42272652884848,
  42.491834194569535,
  12.984372213600585,
  109.16699091385988,
  7.972004532452814,
  115.95675369980508,
  31.60747703140596,
  8.970362227862847,
  5.348212292398322,
  170.32447438112524,
  25.39461971029664,
  25.689753420439597,
  128.4108574894461,
  22.923517291466485,
  23.879023214921144,
  32.850704516803845,
  94.241834519035,
  47.674956807162374,
  42.874078537524895,
  267.69045911495505,
  26.979484728330867,
  7.693118407566914,
  26.129684498104048,
  79.94621217989265,
  72.29553750566461,
  50.67407470133528,
  11.158528124498075,
  82.79205761535891,
  19.157737790208685,
  17.32440397627579,
  20.555302571995973,
  10.529247429324672,
  16.617831539622564,
  79.30467464529444,
  9.586556529180305,
  27.691560468355817,
  15.244960173626433,
  14.156378814162549,
  28.312833940843245,
  144.42400385897858,
  30.843951686429108,
  28.503842966598675,
  33.73055398605692,
  11.717792275688392,
  59.730785744413424,
  98.7516147375063,
  27.017210944780526,
  72.69447911598496,
  29.387392340651942,
  55.48065694190197,
  25.598014438464617,
  56.37725560670995,
  10.470496341939722,
  10.746776600343148,
  10.513823457990869,
  48.16220169554187,
  88.87974662269184,
  43.000628041814466,
  45.45810061817287,
  10.313405625864716,
  5.8587784221719765,
  56.65108657998005,
  63.16020871410067,
  5.744059779469008,
  19.26888098624897,
  12.711023549112907,
  33.50586661377136,
  5.171106749425389,
  83.62372887334072,
  5.076570316611948,
  9.80967873342721,
  32.90073091688129,
  71.11619460163762,
  7.83566819401059,
  135.74274524315902,
  97.51342640508696,
  13.21626051679822,
  19.231859023579307,
  86.74391518532043,
  18.378008196723776,
  70.00226757629704,
  27.03744922006843,
  65.00824796111104,
  81.69910871291269,
  150.24029501106264,
  13.5078216913279,
  87.52256073153305,
  41.99929853602795,
  16.840500379771875,
  16.40353747251429,
  12.725871531799617,
  24.556200229844272,
  38.991822233074046,
  9.1743961229883,
  12.004724580251052,
  11.30407325217084,
  28.82335480698434,
  76.9413871671114,
  12.565114759807928,
  12.244217241795173,
  16.968679101994923,
  10.284049698143138,
  5.568490345122137,
  38.709524829487066,
  13.414745278483165,
  24.467206737540007,
  18.32227261747431,
  9.351399319684976,
  16.998101233912426,
  7.259660085879079,
  94.82048750778833,
  26.50402712375147,
  60.45146088997449,
  24.04368896093938,
  34.22015423624592,
  49.02417389527286,
  12.70207873920776,
  6.824123533046316,
  60.22696514383394,
  59.03014922333422,
  37.33383761241198,
  33.529222118726295,
  23.055007665804016,
  25.470352667080007,
  40.60119960686165,
  18.480932845630722,
  24.662584614048853,
  21.467464795378213,
  24.434821051105246,
  25.75149895464371,
  28.708874142998877,
  11.71658061440018,
  8.67115210440233,
  5.701911266958268,
  33.470786470787885,
  103.39594277767675,
  129.0272434588124,
  6.103257190964134,
  53.80152999506884,
  15.18107246135871,
  13.57947767131657,
  83.12847209064773,
  37.91035844579114,
  5.035424514497084,
  90.07608276274134,
  20.137724540153233,
  15.850670565731386,
  104.74351913363358,
  54.972242220134206,
  19.345840830152646,
  52.315840961158194,
  21.827068816322328,
  69.80099229568637,
  13.284871379796368,
  12.456256639028917,
  11.229390706034335,
  12.878889905664694,
  14.296897577943332,
  77.82817288487384,
  6.582211199591967,
  13.665549855630045,
  77.7139704877213,
  9.055668241211315,
  5.772428704689789,
  71.80452160687348,
  17.060298215282053,
  101.57377530970028,
  5.676264744774738,
  62.753506469608396,
  74.6317907709306,
  47.66046472211516,
  12.468472489162199,
  5.772712952145254,
  15.866236108754977,
  5.289206764185928,
  9.45873791405088,
  63.797398144249144,
  110.19268323613245,
  89.99823758314976,
  16.12767166966444,
  20.486722472118874,
  25.33573871581273,
  82.96305957983394,
  34.740123416116816,
  8.335345774516425,
  7.550262830344939,
  63.97387455254626,
  11.048962226756503,
  5.622720151163016,
  22.552717093928834,
  38.902615652346185,
  5.80146551371251,
  13.515108358102598,
  21.545861529608064,
  21.17567992026929,
  19.985889717332828,
  33.8836065728687,
  67.14101480271866,
  16.26312299478976,
  24.377611928969046,
  5.173963365721217,
  13.219708422485201,
  22.795897220886978,
  9.588192821459481,
  60.719493319813644,
  31.887376943856992,
  61.63608914382314,
  72.03500696873377,
  45.43001952559559,
  9.258926999080897,
  38.95530840917322,
  10.841160035606705,
  9.593211240830037,
  16.7112933995029,
  43.773150754746744,
  21.158617058540894,
  69.0223895100502,
  63.2394972887217,
  5.596985157989508,
  87.56914860691914,
  58.911324525149524,
  13.70960487998605,
  14.920876919768775,
  6.121173059664683,
  6.131791039219603,
  91.84862022944479,
  5.43822957912523,
  8.740972912735097,
  68.20349774504608,
  45.2005446764171,
  15.588788948110123,
  90.53629904683305,
  5.927689844767063,
  17.496416889395206,
  34.3152364079604,
  7.1218511547982555,
  354.39322684690893,
  167.89540900358708,
  17.66175664818565,
  18.934687567299672,
  158.36026077864338])

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)