#!/usr/bin/python

#Deposition for the FIRST sample

from process import *
import odict
import time

tcsdir = "../../research/TCS/"
imagedir = "au_deposition/"

smooth = [4,7]

os.system("mkdir -p " + imagedir.strip("/"))

toplot = odict.odict([
	("Si (initial)" , [tcsdir + "2012-08-09/074447.dat",1]),
	("+Au (3min)", [tcsdir + "2012-08-09/121616.dat", 1]),
	("+Au (8min)", [tcsdir + "2012-08-09/161028.dat", 1]),
	("+Au (8min) a 7hrs", [tcsdir + "2012-08-09/233758.dat", 1]),
	("+Au (8min) a 15hrs", [tcsdir + "2012-08-10/065228.dat", 1]),
	("+Au (18min)", [tcsdir + "2012-08-10/104218.dat", 1]),
	("+Au (18min) a 4hrs", [tcsdir + "2012-08-10/142746.dat", 1]),
	("+Au (18min) a 8hrs", [tcsdir + "2012-08-10/181247.dat", 1]),
	("+Au (18min) a 12hrs", [tcsdir + "2012-08-10/215747.dat", 1]),
	("+Au (18min) a 16hrs", [tcsdir + "2012-08-11/014200.dat", 1]),
	("+Au (18min) a 20hrs", [tcsdir + "2012-08-11/052314.dat", 1]),
	("+Au (18min) a 24hrs", [tcsdir + "2012-08-11/090345.dat", 1]),
	("+Au (18min) a 28hrs", [tcsdir + "2012-08-11/124357.dat", 1]),
	("+Au (18min) a 32hrs", [tcsdir + "2012-08-11/162419.dat", 1]),
	("+Au (18min) a 36hrs", [tcsdir + "2012-08-11/234511.dat", 1]),
	("+Au (18min) a 40hrs", [tcsdir + "2012-08-12/032517.dat", 1]),
	("+Au (18min) a 44hrs", [tcsdir + "2012-08-12/070514.dat", 1]),
	("+Au (18min) a 48hrs", [tcsdir + "2012-08-12/104457.dat", 1]),
	("+Au (18min) a 52hrs", [tcsdir + "2012-08-12/142446.dat", 1]),
	("+Au (18min) a 56hrs", [tcsdir + "2012-08-12/180443.dat", 1]),
	("+Au (18min) a 60hrs", [tcsdir + "2012-08-12/214440.dat", 1]),
	("+Au (18min) a 96hrs", [tcsdir + "2012-08-14/094851.dat", 1]),
	("+Au (33min)", [tcsdir + "2012-08-16/124728.dat", 1]),
	#("+Au (33min) a 4hrs", [tcsdir + "2012-08-16/163328.dat", 1]),	
	#("+Au (33min) a 5hrs", [tcsdir + "2012-08-16/170616.dat", 1]),
	#("+Au (33min) a 18hrs", [tcsdir + "2012-08-17/062341.dat", 1]),
	("+Au (50min)", [tcsdir + "2012-08-17/074929.dat", 1]),
	("+Au (50min) a 1hr", [tcsdir + "2012-08-17/084544.dat", 1]),
	("+Au (83min) a 57hr", [tcsdir + "2012-08-19/172754.dat", 1]),	
	("+Au (100min)", [tcsdir + "2012-08-20/au_on_si(50min)/103405.dat", 1]),
	("+Au (100min) a 1min", [tcsdir + "2012-08-20/103534.dat", 1]),
	("+Au (100min) a 1hr", [tcsdir + "2012-08-20/113127.dat", 1]),
	("+Au (100min) a 2hrs", [tcsdir + "2012-08-20/124236.dat", 1]),
	#("+Au (100min) a 4hrs", [tcsdir + "2012-08-20/143917.dat", 1]),
	("+Au (125min)", [tcsdir + "2012-08-20/153656.dat", 1]),
	("+Au (125min) a 20hrs", [tcsdir + "2012-08-21/191221.dat", 1]),
	("+Au (125min) a 7days", [tcsdir + "2012-08-27/124551.dat", 1]),
	("+Au (145min)", [tcsdir + "2012-08-27/141116.dat", 1])	
])

# I(E) (normalised) data
gnuplot.reset()
gnuplot("set term postscript colour")
gnuplot("set output \""+imagedir+"normalised_data.eps\"")
gnuplot("set title \"Deposition of Au on Si I(E) curves\"")
gnuplot("set key top left")

for f in toplot.items():
	d = ShowData(f[1][0], title=f[0], plot=gnuplot.plot, calibrate=True, normalise=True)



time.sleep(1)

# TCS (normalised) data
gnuplot.reset()
gnuplot("set term postscript colour")
gnuplot("set output \""+imagedir+"tcs_normalised.eps\"")
gnuplot("set title \"Deposition of Au on Si S(E) curves\"")
gnuplot("set key top right")

for f in toplot.items():
	d = ShowTCS(f[1][0], title=f[0], plot=gnuplot.plot, calibrate=True, normalise=True, with_="lp", step=f[1][1], smooth=smooth)


time.sleep(1)

#os.system("cd "+imagedir+"; for i in *.eps; do epstopdf --autorotate=All $i; done; rm *.eps")