# -*- coding: utf-8 -*-
"""
Created on Sun Aug 28 20:01:19 2022

@author: Yale
"""

path = r'/Volumes/AMBER1/21Jan2022.sptw/ptu/hdf5'
path = r'/Volumes/AMBER1/27Jan2022.sptw/ptu/hdf5'
path = r"/Volumes/AMBER1/9Feb2022/ptu/DDX3X/HDF5s"
path = r'/Volumes/AMBER1/9Feb2022/ptu/DDX3Y/HDF5s/hdf5'
path = r'/Volumes/AMBER1/17Feb2022.sptw/ptu/8uMDDX3Y/HDF5s/DDX3Y_ATP'
path = r'/Volumes/AMBER1/20May2022_2uMDDX3.sptw/ptu/2uM_DDX3Y/HDF5s'
path = r'/Volumes/AMBER1/05Apr2022.sptw/ptu/2uM_DDX3Y/HDF5s/ATP'
path = r'/Volumes/AMBER1/07Mar2022_1uMX.sptw/ptu/HDF5s/1uMDDX3X_ATP/hdf5'
path = r'/Volumes/AMBER1/30Mar2022.sptw/ptu/HDF5s'
path = r'/Volumes/AMBER1/25May2022.sptw/ptu/8uMX/HDF5s'
path = r'/Volumes/AMBER1/09Mar2022.sptw/ptu/1uMDDX3Y/HDF5s'
path = r'/Volumes/AMBER1/07Apr2022.sptw/ptu/HDF5s/1uMDDX3Y_ATP'
path = r'/Volumes/AMBER1/24May2022.sptw/ptu/1uMY/HDF5s'
path = r'/Volumes/AMBER1/16Mar2022.sptw/ptu/8uMDDX3X/HDF5s'
path = r'/Volumes/AMBER1/27Apr2022.sptw/ptu/8uMY/HDF5s '
path = r'/Volumes/AMBER1/05Apr2022.sptw/ptu/6uM_DDX3Y/HDF5s/ATP/hdf5'
path = r'/Volumes/AMBER1/27Apr2022.sptw/ptu/6uMY/HDF5s'
path = r'/Volumes/AMBER1/04May2022.sptw/ptu/HDF5s/6uMY'
path = r'/Volumes/AMBER1/26Apr2022.sptw/ptu/6uMY_ATP/HDF5s'
path = r'/Volumes/AMBER1/13Oct2021_Mutants.sptw/hdf5/hdf5'
path = r'/Volumes/AMBER1/02Mar2022.sptw/ptu/2uMDDX3X/HDF5s'
path = r'/Volumes/AMBER1/06Apr2022.sptw/ptu/0.5uM_DDX3X/HDF5s'
path = r'/Volumes/AMBER1/02Mar2022.sptw/ptu/2uMDDX3X/HDF5s'
path = r'/Volumes/AMBER1/21Jan2022.sptw/ptu/hdf5/DL_LF'
path = r'/Volumes/AMBER1/30Mar2022.sptw/ptu/HDF5s'
path = r'/Volumes/AMBER1/07Apr2022.sptw/ptu/HDF5s/1uMDDX3Y/FCS_HF'
path = r'/Volumes/AMBER1/25May2022.sptw/ptu/1uMX/HDF5s/FCS_HF'
path = r'/Volumes/AMBER1/27Apr2022.sptw/ptu/8uMY/HDF5s/FCS_HF'
path = r'/Volumes/AMBER1/25May2022.sptw/ptu/8uMX/HDF5s/FCS_HF'
path = r'/Volumes/AMBER1/09Mar2022.sptw/ptu/1uMDDX3Y/HDF5s/FCS_HF'
path = r'/Volumes/AMBER1/17Feb2022.sptw/ptu/8uMDDX3Y/HDF5s/DDX3Y/FCS_DL_HF'
path = r'/Volumes/AMBER1/17Feb2022.sptw/ptu/8uMDDX3Y/HDF5s/DDX3Y_ATP/FCS_DL_LF'
path = r"/Volumes/AMBER1/Hex_102022.sptw/HDF5s/xxx"
path = r"/Volumes/AMBER1/01Nov2022_nucl.sptw/HDF5s"
path = r"/Volumes/AMBER1/10Nov2022_deltaNC.sptw/HDF5s"
path = r"/Volumes/AMBER1/08Nov2022_nucl.sptw/HDF5s"
path = r"/Volumes/AMBER1/DelDel_anneal_24Oct2022.sptw/HDF5s/HDF5s0"
path = r"/Volumes/AMBER1/02Mar2022.sptw/ptu/2uMDDX3Y/HDF5s"
path = r"/Volumes/AMBER1/18Nov2022_deltaNC_WT.sptw/HDF5s/HDF5s0"
path = r"/Volumes/AMBER1/17Nov2022_ss_ssMBP_newRNA.sptw/GroupMeas_3/HDF5s"
path = r"/Volumes/AMBER1/17Nov2022_ss_ssMBP_newRNA.sptw/HDF5s"
path = r"/Volumes/AMBER1/17Nov2022_ss_ssMBP_newRNA.sptw/GroupMeas_4/HDF5s"
path = "/Volumes/AMBER1/18Nov2022_deltaNC_WT.sptw/GroupMeas_4/HDF5s"
path = r"/Volumes/AMBER1/16Dec2022_nucl.sptw/HDF5s/csvs/FCS"
path = r"/Volumes/AMBER1/15Dec2022_deltaCdeltaN.sptw/HDF5s"
path = r"/Volumes/AMBER1/15Dec2022_deltaCdeltaN.sptw/HDF5s/HDF5s_deltaNDDX3Y_ATP"
path = r"/Volumes/AMBER1/29Jun2022.sptw/HDF5s/HDF5s1"
path = r"/Volumes/AMBER1/29Dec2022_domainswap.sptw/HDF5s"
path = r"/Volumes/AMBER1/10Jan2023_delta_MBP.sptw/HDF5s/HDF5s4_2uMdeltadeltaX"
path = r"/Volumes/AMBER4/26Jan2023_mCherryNucl.sptw/HDF5s"
path = r"/Volumes/AMBER4/20April2023_mCherry_RNA.sptw/HDF5s"
path = r"/Volumes/AMBER4/13June2023_Anneal_MCOtrunc.sptw/HDF5s/HDF5s8"
path = r"/Volumes/AMBER1/10Jan2023_delta_MBP.sptw/HDF5s"
path = r"/Volumes/AMBER4/29June2023_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER4/28June2023_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER4/26June2023_Mut532.sptw/HDF5s"
path = r"/Volumes/AMBER4/18July2023_TruncMCO.sptw/GroupMeas_19/HDF5s"
path = r"/Volumes/AMBER4/28June2023_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER4/26June2023_Mut532_R534H_L505V.sptw/HDF5s"
path = r"/Volumes/AMBER4/29June2023_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER4/12July2023_Mutant.sptw/HDF5s"
path = r"/Volumes/AMBER4/13July_2023_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER4/31Aug2023_G325EWT.sptw/HDF5s"
path = r"/Volumes/AMBER4/24Aug2023_R351W_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER4/23Aug2023_Mut_A497V.sptw/HDF5s"
path = r"/Volumes/AMBER4/31Aug2023_G325EWT.sptw/HDF5s"
path = r"/Volumes/AMBER4/22Aug2023_Mut_T275M.sptw/HDF5s"
path = r"/Volumes/AMBER4/04Aug2023_E449G.sptw/HDF5s"
path = r"/Volumes/AMBER4/28June2023_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER4/27July023_R475G_Y.sptw/HDF5s"
path = r"/Volumes/AMBER1/10Mar2022.sptw/ptu/HDF5s"
path = r"/Volumes/AMBER4/20Oct2023_E449_dd.sptw/GroupMeas_3/HDF5s"


import matplotlib.pyplot as plt
import numpy as np

fN = path + "results.txt"

fH = open(path + '/' + 'results.txt', 'r')

results = []
runNames = []   # list of run Names
tauDs = []   # list of tauDs
Ns = []      # list of ns
nRuns = 0    # number of runs
for line in fH:
#    print (line)
    results.append(line)

fH.close()
    
lenRes = len(results)

for k in range(lenRes):
#    print (results[k])
    if ('%%' in results[k]):
#        print (results[k][3:])
        nRuns += 1
        runNames.append(results[k][3:-6])
    if ('tauD0:' in results[k]):
#        print (results[k])
        tauDs.append(float(results[k].split()[1]))
    if ('N0:' in results[k]):
#        print (results[k])
        Ns.append(float(results[k].split()[1]))
        

print ("runNames")
print (runNames) 
print (' ')       
print ("tauDs")
print (tauDs)
print (' ')       
print ("Ns")        
print (Ns)        

x = np.arange(nRuns) 

fig1 = plt.figure(figsize=(6, 4), dpi=600)
ax1 = plt.subplot()
ax1.set_ylabel("tauD (s)")
ax1.plot(runNames, tauDs, "o", markersize = 5)
ax1.plot(runNames, tauDs, "-", markersize = 5)
plt.setp(ax1.get_xticklabels(), rotation=30, ha='right')
ax1.set_ylim(0, 1.2*max(tauDs))
plt.show()

fig2 = plt.figure(figsize=(6, 4), dpi=600)
ax2 = plt.subplot()
ax2.set_ylabel("N")
plt.plot(runNames, Ns, "o", markersize = 5)
plt.plot(runNames, Ns, "-", markersize = 5)
plt.setp(ax2.get_xticklabels(), rotation=30, ha='right')
ax2.set_ylim(0, 1.2*max(Ns))
plt.show()