set termoption dash

set style line 1 lt 2 lc rgb "black" lw 1
set style line 2 lt 8 lc rgb "black" lw 1
set style line 3 lt 4 lc rgb "black" lw 1
set style line 4 lt 7 lc rgb "black" lw 1
set style line 5 lt 7 lc rgb "blue" lw 1
set style line 6 lt 1 lc rgb "black" lw 2
set style line 7 lt 1 lc rgb "blue" lw 2

set key font "Helvetica,16"

set xlabel 'Polarizer Orientation, {/Symbol q} ({/Symbol \260})' font "Helvetica,18"
set ylabel 'Probability' font "Helvetica,18"
show xlabel
show ylabel

set grid
set xtics ("-90{/Symbol \260} (V)" -90, "-45{/Symbol \260}" -45, "0{/Symbol \260} (H)" 0, "+45{/Symbol \260}" 45, "90{/Symbol \260} (V)" 90) font "Helvetica,16"
#set ytics -0.2,0.1,0.4 font "Helvetica,16"
set xrange [0:90]
set yrange [0:0.5]

#set label 1 '{/Symbol a}_1 = +45{/Symbol \260}' at -6.5,210
#set label 2 '{/Symbol a}_1 = 0{/Symbol \260}' at -34,160
#set arrow 1 from -7.5,210 to -22.5,210 filled
#set arrow 2 from -35,160 to -48,160 filled

a(x) = a1*sin(x*2*3.142/(x*a4 + a5)) + a2*cos(x*2*3.142/(x*a4 + a5)) + a3
b(x) = b1*sin(x*2*3.142/(x*b4 + b5)) + b2*cos(x*2*3.142/(x*b4 + b5)) + b3
c(x) = c1*sin(x*2*3.142/(x*c4 + c5)) + c2*cos(x*2*3.142/(x*c4 + c5)) + c3
d(x) = d1*sin(x*2*3.142/(x*d4 + d5)) + d2*cos(x*2*3.142/(x*d4 + d5)) + d3

a1 = -4.14854
a2 = -136.198
a3 = 136.426
a4 = -0.0654965
a5 = 178.079

#a1 = -4.14854         +/- 0.3896       (9.39%)
#a2 = -136.198         +/- 0.5569       (0.4089%)
#a3 = 136.426          +/- 0.5616       (0.4117%)
#a4 = -0.0654965       +/- 0.008755     (13.37%)
#a5 = 178.079          +/- 0.6592       (0.3702%)

b1 = -94.7849
b2 = -91.5479
b3 = 131.855
b4 = 0.00410186
b5 = 179.902

#b1 = -94.7849         +/- 0.2527       (0.2666%)
#b2 = -91.5479         +/- 0.4536       (0.4955%)
#b3 = 131.855          +/- 0.3441       (0.261%)
#b4 = 0.00410186       +/- 0.006881     (167.8%)
#b5 = 179.902          +/- 0.3567       (0.1983%)

c1 = -132.596
c2 = 3.43297
c3 = 132.712
c4 = -0.013522
c5 = 178.246

#c1 = -132.596         +/- 0.402        (0.3032%)
#c2 = 3.43297          +/- 0.5435       (15.83%)
#c3 = 132.712          +/- 0.3969       (0.2991%)
#c4 = -0.013522        +/- 0.006282     (46.46%)
#c5 = 178.246          +/- 0.2525       (0.1417%)

d1 = -96.3166
d2 = 98.5453
d3 = 137.926
d4 = -0.00777243
d5 = 179.982

#d1 = -96.3166         +/- 0.5134       (0.5331%)
#d2 = 98.5453          +/- 0.3717       (0.3771%)
#d3 = 137.926          +/- 0.3134       (0.2272%)
#d4 = -0.00777243      +/- 0.004687     (60.31%)
#d5 = 179.982          +/- 0.1761       (0.09786%)

a_max = (sqrt(a1*a1 + a2*a2) + a3)*2
b_max = (sqrt(b1*b1 + b2*b2) + b3)*2
c_max = (sqrt(c1*c1 + c2*c2) + c3)*2
d_max = (sqrt(d1*d1 + d2*d2) + d3)*2

#plot '0_0_visibility_long_integration_20141126.dat' using 1:($5/10)/a_max ls 4 title '',\
#     '1_229227_visibility_long_integration_20141126.dat' using 1:($5/10)/b_max ls 5 title '',\
#     '0_454005_visibility_long_integration_20141126.dat' using 1:($5/10)/c_max ls 4 title '',\
#     '1_675669_visibility_long_integration_20141126.dat' using 1:($5/10)/d_max ls 5 title '',\
#      (a(x)/a_max) ls 6 title '', (b(x)/b_max) ls 7 title '', (c(x)/c_max) ls 6 title '', (d(x)/d_max) ls 7 title ''


plot (a(x)/a_max) ls 6 title '', (b(x)/b_max) ls 7 title '', (c(x)/c_max) ls 6 title '', (d(x)/d_max) ls 7 title ''
      
set terminal postscript enhanced color
set output 'wedge_error_characterisation_curve.eps'
replot