startcount = 1
r = pd.DataFrame
allmatrix = []
for j in range(len(pathstring)):
for p in sorted(
Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.out')):
print(p)
try:
r_single = mu.getAvgR(
str(p)) #, 'zero', 'cylinder_shell.out', p.parts[-2][2])
inseg = mu.gettubeILTT(p, 'adsorbed') # choose itube/otube/iotube
outseg = mu.gettubeOLTT(
p, 'adsorbed'
) # [[# adsorbed chain, fluc], [# adsorbed seg, fluc]]
except ValueError:
continue
# print(inseg)
# print('#adsorbed chain', inseg[0][0])
# print('#adsorbed chain fluc', inseg[0][1])
# print('#adsorbed seg', inseg[1][0])
# print('#adsorbed seg fluc', inseg[1][1])
# print(outseg)
labelstring = 'ph = {:d}'.format(int(p.parts[-2][2]))
print(p.parts)
print(p.parts[-2])
print(p.parts[-2][2])
seg_io_loop_errup = []
num_io_loop = []
num_io_loop_errdown = []
num_io_loop_errup = []
ph = []
for p in sorted(
Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.out')):
# zip(pathstring, linecolor)
print(p)
#io_loops = mu.gettubeIOLTT(p, 'loops') #choose itube/otube/iotube
#io_tails = mu.gettubeIOLTT(p, 'tails')
#io_trains = mu.gettubeIOLTT(p, 'trains')
i_loops = mu.gettubeILTT(p, 'loops') #choose itube/otube/iotube
o_loops = mu.gettubeOLTT(p, 'loops') #choose itube/otube/iotube
io_loops = i_loops + o_loops
## Input: fin - file in
## restype - type of result. Can be 'loops', 'tails', 'trains' or 'adsorbed'
## Output: 2D array of results on the form
## [number of x] [variation (fluctuation)],
## [length of x] [variation],
## [number of segments in x] [variation]
## in the case of x = 'loops', 'tails' or 'trains', or on the form
## [number of adsorbed chains] [variation],
## [number of adsorbed segments] [variation]
## in the case of 'adsorbed'
## Note: loops, tails and trains are average number per (adsorbed?) chain!
if j == 0:
for p in sorted(
Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.out')):
# zip(pathstring, linecolor)
print(p)
io_loops = mu.gettubeIOLTT(p, 'loops') #choose itube/otube/iotube
io_tails = mu.gettubeIOLTT(p, 'tails')
io_trains = mu.gettubeIOLTT(p, 'trains')
i_loops = mu.gettubeILTT(p, 'loops') #choose itube/otube/iotube
i_tails = mu.gettubeILTT(p, 'tails')
i_trains = mu.gettubeILTT(p, 'trains')
o_loops = mu.gettubeOLTT(p, 'loops') #choose itube/otube/iotube
o_tails = mu.gettubeOLTT(p, 'tails')
o_trains = mu.gettubeOLTT(p, 'trains')
## Input: fin - file in
## restype - type of result. Can be 'loops', 'tails', 'trains' or 'adsorbed'
## Output: 2D array of results on the form
## [number of x] [variation (fluctuation)],
## [length of x] [variation],
## [number of segments in x] [variation]
## in the case of x = 'loops', 'tails' or 'trains', or on the form
## [number of adsorbed chains] [variation],
## [number of adsorbed segments] [variation]
## in the case of 'adsorbed'
## Note: loops, tails and trains are average number per (adsorbed?) chain!
seg_io_loop_errup = []
num_io_loop = []
num_io_loop_errdown = []
num_io_loop_errup = []
ph = []
for p in sorted(Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.out')):
# zip(pathstring, linecolor)
print(p)
#io_loops = mu.gettubeIOLTT(p, 'loops') #choose itube/otube/iotube
#io_tails = mu.gettubeIOLTT(p, 'tails')
#io_trains = mu.gettubeIOLTT(p, 'trains')
i_loops = mu.gettubeILTT(p, 'adsorbed') #choose itube/otube/iotube
o_loops = mu.gettubeOLTT(p, 'adsorbed') #choose itube/otube/iotube
io_loops = i_loops + o_loops
if j == 0:
titlestring = 'neg, out'
elif j == 1:
titlestring = 'pos, out'
elif j == 2:
titlestring = 'neg, in'
elif j == 3:
titlestring = 'pos, in'
#titlestring = '{:s}'.format(str(p.parts[-4]))
labelstring = 'pH = {:d}'.format(int(p.parts[-2][2]))