def Oct17():
a= loadtxt('/home/chris/Dropbox/DataWeiss/151017/thioldots.csv', delimiter = ',', unpack = True, skiprows=1, usecols = (0,6,3,1,2,4,5))
for i in range(1,7):
r = polyfit(a[0,0:100],a[i,0:100],1)
a[i]-=rs.polyeval(r,a[0])
peak = uv.findpeak(a[0],a[i],(410,420))
#print peak
diameter = -0.000000066521*peak[0]**3+0.00019557*peak[0]**2-0.092352*peak[0]+13.29
#print 'diam', diameter
epsilon = 21536*diameter**2.3
# print 'eps', epsilon
print 'CONC', 5*peak[1]/epsilon
plot(a[0],a[i])
legend(['4.2','6.0',',6.4','7.7','9.7','11.5'])
print 'THIOL EXCHANGE'
a= loadtxt('/home/chris/Dropbox/DataWeiss/151017/thiolcapped doits.csv', delimiter = ',', unpack = True, skiprows=1, usecols = (0,1))
r = polyfit(a[0,0:100],a[1,0:100],1)
a[1]-=rs.polyeval(r,a[0])
peak = uv.findpeak(a[0],a[1],(410,420))
#print peak
diameter = -0.000000066521*peak[0]**3+0.00019557*peak[0]**2-0.092352*peak[0]+13.29
#print 'diam', diameter
epsilon = 21536*diameter**2.3
# print 'eps', epsilon
print 'CONC', 5*peak[1]/epsilon
a= loadtxt('/home/chris/Dropbox/DataWeiss/151017/thioldots.csv', delimiter = ',', unpack = True, skiprows=1, usecols = (0,3,1,2,4,5,6))
for i in range(1,7):
r = polyfit(a[0,0:100],a[i,0:100],1)
a[i]-=rs.polyeval(r,a[0])
peak = uv.findpeak(a[0],a[i],(410,420))
#print peak
diameter = -0.000000066521*peak[0]**3+0.00019557*peak[0]**2-0.092352*peak[0]+13.29
#print 'diam', diameter
epsilon = 21536*diameter**2.3
# print 'eps', epsilon
print 'CONC', 5*peak[1]/epsilon
plot(a[0],a[i])
legend(['200','400',',800','2000','3200','4000'])
return 0
def Oct13():
a= loadtxt('/home/chris/Dropbox/DataWeiss/151013/NMR samples.csv', delimiter = ',', unpack = True, skiprows=1)
for i in range(1,6):
r = polyfit(a[0,0:100],a[i,0:100],1)
a[i]-=rs.polyeval(r,a[0])
peak = uv.findpeak(a[0],a[i],(410,420))
print peak
diameter = -0.000000066521*peak[0]**3+0.00019557*peak[0]**2-0.092352*peak[0]+13.29
print 'diam', diameter
epsilon = 21536*diameter**2.3
print 'eps', epsilon
print 'conc', 5*peak[1]/epsilon
plot(a[0],a[i])
legend(['10.5','9.5','9.3','8.6','6.8'])
return 0
