#set term x11 persist
set term cairolatex pdf size 5in,3in
set output "rmds.tex"

ptsize = 0.2

# R ColorBrewer theme (minus yellow)
set linetype 1 lc rgb '#E41A1C' pt 7 ps ptsize
set linetype 2 lc rgb '#377EB8' pt 7 ps ptsize
set linetype 3 lc rgb '#4DAF4A' pt 7 ps ptsize
set linetype 4 lc rgb '#984EA3' pt 7 ps ptsize
set linetype 5 lc rgb '#FF7F00' pt 7 ps ptsize
set linetype 6  lc rgb '#A65628' pt 7 ps ptsize 
set linetype 7  lc rgb '#F781BF' pt 7 ps ptsize 
set linetype 8  lc rgb '#999999' pt 7 ps ptsize 
set linetype 9 lc rgb '#8DD3C7' pt 7 ps ptsize 
set linetype 10 lc rgb '#BEBADA' pt 7 ps ptsize 
set linetype 11 lc rgb '#FB8072' pt 7 ps ptsize 

set format x "$10^{%L}$"
set format y "$10^{%L}$"

set xlabel "Simulation Time [s]"
set ylabel "RMSD [m]"
set offsets graph 0.1, graph 0.1, graph 0.1, graph 0.1
set logscale xy
set key top left
set grid

set arrow from 2.0,0.00001 to 2.0,0.42 nohead lt 1 dt 2
set arrow from 0.10,0.00001 to 0.10, 0.032 nohead lt 2 dt 2
set arrow from 0.010,0.00001 to 0.010, 0.00268 nohead lt 3 dt 2
set arrow from 0.0010,0.00001 to 0.0010, 0.00032 nohead lt 4 dt 2

plot "g05means.out" with linespoints lt 1 smooth csplines title "$\\gamma = 0.5$", \
    "g05means.out" with points lt 1  notitle, \
    "g10means.out" with linespoints lt 2 smooth csplines title "$\\gamma = 1.0$", \
    "g10means.out" with points lt 2  notitle,\
    "g100means.out" with linespoints lt 3 smooth csplines title "$\\gamma = 100.0$", \
    "g100means.out" with points lt 3  notitle, \
    "g1000means.out" with linespoints lt 4 smooth csplines title "$\\gamma = 1000.0$", \
    "g1000means.out" with points lt 4 notitle, \
    0.5*x with lines dt 2 lc rgb "black" notitle