#!/usr/bin/gnuplot

reset

#set terminal postscript enhanced
#set output "slide-zavislost-vdrift-E-beta-model.ps"
set terminal postscript eps enhanced 20
set output "slide-zavislost-vdrift-E-beta-model.eps"

set line style 1 linetype 7 linewidth 1
set line style 2 linetype 1 linewidth 3
set line style 3 linetype 4 linewidth 0.1

set logscale xy
set rmargin 2
set tmargin 0
set bmargin 3
set lmargin 9


# fit [1e-4:1e-2] [] f(x) "zavislost-vdrift-E_T-300-normal-model.dat" using 1:2 via a; replot
# a               = 0.667449         +/- 0.002494     (0.3736%)
# 
# fit [1e-4:1e-3] [] f(x) "zavislost-vdrift-E_T-300-normal-model.dat" using 1:2 via a; replot
# a               = 0.653661         +/- 0.002711     (0.4147%)
# fit [1e-3:1e-2] [] f(x) "zavislost-vdrift-E_T-300-normal-model.dat" using 1:2 via a; replot
# a               = 0.668264         +/- 0.002686     (0.4019%)

k = 891.053540967816
a = 0.668264

set multiplot

set key top left Left reverse samplen 4
set noxtics
set noytics

set pointsize 1.2

plot [:10][2e-2:4e3] \
	"zavislost-vdrift-E_T-300-beta-0-model.dat"		title "{/Symbol b} = 0.00" with lp 2, \
	"zavislost-vdrift-E_T-300-beta-0.05-model.dat"	title "{/Symbol b} = 0.05" with lp 3, \
	"zavislost-vdrift-E_T-300-beta-0.1-model.dat"	title "{/Symbol b} = 0.10" with lp 4, \
	"zavislost-vdrift-E_T-300-beta-0.15-model.dat"	title "{/Symbol b} = 0.15" with lp 5

set grid linestyle 3
set arrow from 1e-2,graph 0 to 1e-2,graph 0.5 nohead linestyle 3
set arrow from 2e-3,1e0  to 7e-2,1e0  nohead linestyle 3
set arrow from 2e-3,1e-1 to 7e-2,1e-1 nohead linestyle 3

set key bottom right Right noreverse samplen 4
set noxtics
#set noytics
set xlabel "{/Symbol a}"
#set ylabel "v_{drift} [m.s^{-1}]"
set label "v_{drift} [m.s^{-1}]" at graph -0.15, graph 0.5 center rotate

set xtics ("10^{-4}" 1e-4, "10^{-3}" 1e-3, "10^{-2}" 1e-2, "10^{-1}" 1e-1, "1" 1, "10^1" 10, "10^2" 1e2 )
set ytics ("10^4" 1e4, "10^3" 1e3, "10^{-2}" 1e-2, "10^{-1}" 1e-1, "10^0" 1, "10^1" 10, "10^2" 1e2 )


plot [:10][2e-2:4e3] \
	"zavislost-vdrift-E_T-300-beta-0.2-model.dat"	title "{/Symbol b} = 0.20" with lp 6, \
	"zavislost-vdrift-E_T-300-beta-0.25-model.dat"	title "{/Symbol b} = 0.25" with lp 7, \
	"zavislost-vdrift-E_T-300-beta-0.5-model.dat"	title "{/Symbol b} = 0.50" with lp 8, \
	"zavislost-vdrift-E_T-300-beta-1.0-model.dat"	title "{/Symbol b} = 1.00" with lp 9, \
	k * sqrt(x) notitle  with lines 1, \
	k * x notitle  with lines 1