# -*- 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/HDF5s9"
path = r"/Volumes/AMBER1/29Dec2022_domainswap.sptw/HDF5s"
path = r"/Volumes/AMBER4/27Jan2023_mCherryNuclY.sptw/HDF5s/HDF5s6_50nMDDX3Y"
path = r"/Volumes/AMBER4/16Jan2023_MCO_RNA.sptw/HDF5s"
path = r"/Volumes/AMBER4/10March2023_MCORNA.sptw/HDF5s"
path = r"/Volumes/AMBER4/22March2023_MCOtrunc.sptw/HDF5s"
path = r"/Volumes/AMBER4/10March2023_MCORNA.sptw/HDF5s/HDF5s6_200nMDDX3X_2mMATP3"
path = r"/Volumes/AMBER4/24March2023_MCORNA.sptw/GroupMeas_9/HDF5s"
path = r"/Volumes/AMBER4/15Dec2022_deltaCdeltaN.sptw/HDF5s/HDF5s_deltaNDDX3Y_ATP"
path = r"/Volumes/AMBER1/DelDel_anneal_24Oct2022.sptw/HDF5s/HDF5s2"
path = r"/Volumes/AMBER4/20April2023_mCherry_RNA.sptw/HDF5s"
path = r"/Volumes/AMBER4/25April2023_mCherryRNA.sptw/HDF5s"
path = r"/Volumes/AMBER4/26April2023_mCherryRNA.sptw/HDF5s"
path = r"/Volumes/AMBER4/13June2023_Anneal_MCOtrunc.sptw/HDF5s/HDF5s1"
path = r"/Volumes/AMBER4/28June2023_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER4/28June2023_Mut.sptw/HDF5s"
path = r"/Volumes/AMBER1/10Mar2022.sptw/ptu/HDF5s"
path = r"/Volumes/AMBER4/20Oct2023_E449_dd.sptw/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 ('tauDa:' in results[k]):
#        print (results[k])
        tauDs.append(float(results[k].split()[1]))
    if ('GN:' 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()